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XIII. Hxtra-Galactic Radio-Frequency Radiation 


By R. Hansury Brown and C, Hazarp 
Jodrell Bank Experimental Station, University of Manchester* 


[Received November 4, 1951] 


SUMMARY 


It has been shown previously that radio-frequency radiation is being 
received on the earth from the extra-galactic nebula N.G.C. 224 (M. 31). 
This paper reports some initial attempts to extend the measurements of 
extra-galactic radio radiation on a wavelength of 1:89m. A _ radio 
survey has been made of three bright nebulae (in addition to N.G.C. 224) 
and three radio sources have been associated with them. The radio 
intensity to be expected from eight of the major clusters of nebulae has 
been calculated and two of the clusters have been surveyed and identified 
with radio sources. For a selected region of sky an attempt has been 
made to predict the irregularities in the radio isophotes corresponding to 
large scale irregularities in the distribution of extra-galactic nebulae. 
These predictions have been compared with radio observations of the 
region, and it is shown that there is a correlation between the radio 
isophotes and the distribution of extra-galactic nebulae. A preliminary 
examination has been made of the relation between the apparent 
photographic magnitude and the apparent radio magnitude of nebulae. 
The results suggest that there is a definite relationship and that the 
apparent photographic magnitude increases more rapidly than the 
apparent radio magnitude. 


§1, INTRODUCTION 


Since the original experiments of Jansky (1931) it has been recognized 
that radio-frequency radiation is being received on the earth from 
extraterrestrial sources. More recent measurements have shown that the 
isophotes of this radiation correspond well with the general outlines of the 
Galaxy as derived from visual observations. The existence of discrete 
sources of radiation, or radio stars, has been established since 1948 
(Bolton and Stanley 1948, Ryle and Smith 1948), and a recent survey of 
the northern hemisphere (Ryle, Smith and Elsmore 1950) gives the 
position of 50 radio stars. Sal , ; 
The results available at present suggest that the majority of the radio 
radiation received has its origin in the radio stars, and that the radio 
stars themselves are widespread throughout the Galaxy with a 
distribution, and even perhaps a density, similar to that of the more 
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common type of star. Attempts to identify these radio stars with 
known visual objects have so far failed and it appears that the radio 
star may represent a new type of astronomical body. 

In a recent work (Hanbury Brown and Hazard 1950) the authors have 
reported the detection of a radio source which has been identified with the 
great spiral nebula in Andromeda (N.G.C. 244, M. 31). This result 
suggests that, whatever the mechanism may be, radio-frequency 
radiation is to be expected from other extra-galactic nebulae. 

The present paper reports some initial attempts to extend the 
measurements of extra-galactic radio-frequency radiation. The experi- 
ments. reported here, have been concentrated on three main lines of 
attack. Firstly, a survey has been made of the radio emission from four 
bright nebulae which lie within the field of the apparatus. Secondly, 
the results from this experiment have been used to predict the 
radio intensity to be expected from the major clusters of nebulae, and 
these predictions have been tested experimentally on two clusters which 
are conveniently within the view of the apparatus. Thirdly, an attempt 
has been made to detect irregularities in the radio isophotes corresponding 
to the large scale irregularities in the distribution of extra-galactic 
nebulae. 


§2. DESCRIPTION OF APPARATUS 


The apparatus used in these experiments has been described in detail 
elsewhere (Hanbury Brown and Hazard 1951), and only a brief description 
will be given here. The aerial system is a paraboloid of 218 ft. diameter. 
and 126 ft. focal length and at a wavelength of 1-89 m the beam width 
is 2 degrees between half power points. 'The beam can be directed to 
declinations between N 38° and N 68° by tilting the mast which supports 
the primary feed. 

The design of the receiving equipment is based on that described by 
Ryle and Vonberg (1948) in which the power received by the aerial is 
balanced continuously against that from a standard generator. The 
receiver is tuned to 1-89 metres and has a .pre-detector bandwidth of 
2 me/s, a post-detector time-constant of 10 seconds and a noise factor 
of 3-5. The calibration of the whole equipment for the measurement of 
power depends upon the aerial gain and upon the standard generator 
which is a tungsten filament diode type CV. 172. 

_ During the daytime the minimum flux which can be detected by this 

equipment is limited by various forms of interference ; however, at night 
it has been found possible to approach the theoretical limits set by the 
fluctuations in the receiver. On favourable nights it is possible to detect 
on a single record a source whose intensity is about 10-* watts/sq.m/c.p.s. 
and by averaging the results of at least four good observations this limit 
has been extended to 5 10-®6 watts/sq.m/c.p.s. Most of the results 
given in this paper were obtained from records made between midnight 
and sunrise. 
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§3. METHOD oF OBSERVATIONS 


The technique of observation was determined by the mechanical 
difficulty of moving the aerial beam rapidly. For the examination of 
a particular region of the sky the aerial beam was fixed at a constant 
elevation, corresponding to the required declination, and at a constant 
azimuth (approximately due south or due north). With the beam at 
this fixed elevation the power received was recorded for the required 
periods of right ascension. The procedure was repeated until an adequate 
number of good records was thade. The measurements were then repeated 
with the beam at different elevations until the particular region of 
the sky had been surveyed. 


§4. THe Derecrion or [NprvipvaL Exrra-Gatactic NEBULAE 


Following the work on N.G.C. 224 (M. 31) a survey was made of the 
region around three of the bright nebulae in the field of view of the aerial. 
These are N.G.C. 5457 (M. 101), N.G.C. 5194 (M. 51) and N.G.C. 4258. In 
each case a radio source was found close to the coordinates of the nebulae. 
The results are shown in table | together with the previously published 
data on N.G.C. 224 (Hanbury Brown and Hazard 1951). 

The presence of a source, probably a radio star, separated from the 
position of N.G.C. 5457 by about 2 degrees made the measurements of 
the weak source at the position of the nebula unreliable, and it was only 
_possible to establish an upper limit to its intensity. 

The identification of the two radio sources with the two nebulae 
N.G.C. 5194 and N.G.C. 4258 can only be made tentatively. The only 
two methods at present available for testing the identification are the 
comparison of the celestial coordinates of the sources with those of the 
nebulae, and the comparison of their intensities with those predicted 
from the other data. It is not possible with the present apparatus to 
measure the apparent diameter of these sources and to correlate them 
with the visual diameters of the nebulae, as was done in the case of 
N.G.C. 224. A comparison of the celestial coordinates shows that the 
positions of the two sources agree with those of the two nebulae within 
the limits of error of the measurements. It is possible, however, that the 
two sources are not associated with the nebulae and are ‘ foreground ’ 
radio stars which are fortuitously coincident in position with the nebulae. 
Ryle, Smith and Elsmore (1950) have published a survey of the radio 
stars in the northern hemisphere, and _ their data has: been used to 
estimate the probability that either of these sources is a radio star. It 
can be shown from their results that the probability of finding a radio 
star whose intensity is comparable (+2 magnitudes) with that of one of 
the sources, and which lies within the limits of error of its position (about 
+1 degree) is about 1 in 50. The probability that both sources represent 
accidental coincidences with radio stars is therefore very small. 
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The radio intensity to be expected from the two nebulae can be 
predicted from the measurements made on N.G.C. 224 by assuming that 
the ratio of the light flux to the radio flux is constant for all nebulae. 
Inspection of table 1 shows that the two radio sources are of about the 
same intensity as would be-expected from the almost equal photographic 
magnitudes of the two nebulae ; however, the observed value of their 
radio intensity is 2 magnitudes more intense in relation to their 
photographic magnitude than would be expected from the measured 
properties of N.G.C. 224. This discrepancy cannot be explained by known 
errors in the measurements, and is discussed in § 7. 


§5. CLUSTERS oF NEBULAE 


The intensity of the radio-frequency radiation received from the 
nebulae given in table 1 suggests that it should be possible to detect 
the radiation from the large clusters and groups of nebulae. The radio 
intensity to be expected from eight large clusters in the northern 
hemisphere has therefore been estimated by assuming that the ratio 
of radio flux to light flux found for N.G.C. 5194 applies to all nebulae. 
The results are expressed on the scale of radio magnitude defined in 
table 1. The integrated photographic magnitude of each cluster has been 
calculated from the number and the most frequent magnitude of the 
component nebulae as given by Hubble and Humason (1931). The 
frequency distribution of the magnitudes of the nebulae in each cluster 
has been assumed to be symmetrical about the most frequent magnitude, 
and the form of the distribution has been obtained from the data given 
for the Coma cluster. The estimated radio magnitudes of the eight 
clusters are given in table 2. 

Zwicky (1937, 1951) has more recently suggested that the extensions 
of the clusters may be greater than estimated by Hubble and Humason. 
If these larger extensions are taken into account no maximum is observed 
in the luminosity function. The maximum found by Hubble and 
Humason is a selective effect due to a concentration of the brighter 
nebulae towards the centre of the cluster. Zwicky’s results do not, 
- however, seriously affect the figures given in table 2. An analysis of 
the Coma cluster shows that the concentration of nebulae, both in density 
and brightness, towards the centre of the cluster is such that radio 
radiation from the outer regions of the cluster would not be detectable, 
and of the nebulae in the central regions only the brighter magnitude 
groups are important. In the case of the Coma cluster, therefore, the 
results given by Hubble and Humason are satisfactory for the present 
estimate, and it will be assumed that this is also true of the other clusters. 

Only four of the clusters given in table 2 lie within the working field 
of the present aerial system and, so far, only two of them have been 
investigated, namely, the Perseus cluster and the cluster centred on 
N.G.C. 911. In each case a radio-frequency source has been detected 
close to the position of the cluster, and the results of the measurements 
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are shown in table 3. Table 3 also gives the coordinates and the 
predicted radio magnitudes of the two clusters. Inspection of the results 
shows that the celestial coordinates of the two sources agree closely with 
those of the corresponding clusters. | 

The observed radio magnitude of the source close to the Perseus 
cluster differs from the predicted radio magnitude of the cluster by 
2 magnitudes, while the radio magnitude of the second source differs by 
less than 2} magnitudes from the predicted value. As the two clusters 
lie in low galactic latitudes, the apparent photographic magnitudes of 
the component nebulae may have been affected by visual obscuration 
in the Galaxy. If the obscuring matter reduces the intensity of the 
light received from the cluster without attenuating the radio flux, then 
the predicted radio magnitudes in table 3 will be too high. The visual 
obscuration for the Perseus cluster has been estimated by Hubble and 
Humason to be 0-4 magnitude on the basis of the abnormal colour excess 
of its component nebulae, and the predicted radio magnitude of this 
cluster may, therefore, be too great by this amount. The second cluster 
lies in the same region and its predicted radio magnitude may be subject 
to the same error. The discrepancy between the observed and the 
predicted radio intensity from the Perseus cluster suggests either that the 
radiation of a significant number of nebulae has been neglected or that 
the radio magnitudes of the component nebulae have been overestimated. 
For example, the results can be explained by assuming that a nebula 
of apparent photographic magnitude +16 has a radio magnitude of 
+14 (on the scale defined in table 1). It follows that all the predicted 
radio magnitudes in table 2 may be in error by the same amount as 
that predicted for the Perseus cluster, namely 2 magnitudes. 

An attempt has been made to measure the apparent diameter of the 
two radio sources. In the case of the Perseus cluster the component 
nebulae lie in a group of diameter about 2 degrees, and it was hoped 
to measure the diameter of the associated radio source by observing the 
apparent broadening of the aerial beam as it swept across the cluster. 
The degree of broadening to be expected has been calculated for the case 
of the beam passing over the Coma cluster. This particular cluster was 
used in the calculation because the apparent angular distribution of its 
component nebulae has been published in detail (Zwicky 1937) and because 
it is of comparable angular diameter to that in Perseus. The calculations 
showed that the maximum degree of broadening of the beam, measured 
between half-power points, would only be about 25% even if the 
concentration of the brighter nebulae to the centre of the cluster is 
neglected. It is to be expected that the broadening would be much less 
in practice due to this concentration of the brighter nebulae. | The 
observed results show that the apparent broadening of the beam is less . 
than 25%, but the records are not sufficiently good to permit a more 
accurate measurement. Hence the measurement of the apparent angular 
diameter of these sources must await more refined measurements, either 
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with a narrower beam or with an interferometer. It should also be noted 
that the value of intensity given for the sources is not a true integrated 
value since it has been assumed that they are point sources. 

It is possible that the sources are, in fact, radio stars in the Galaxy 
and that their positions happen to coincide by chance with those of the 
clusters. From the data given by Ryle, Smith and Elsmore (1950), it 
may be shown that the chance that a radio star falls within the limits 
of error of the position given for the source near the Perseus cluster is 
about 1 in 500. For the weaker source (near N.G.C. 911) this chance 
is about 1 in 100. The probability that both sources represent accidental 
coincidences is therefore very low. From the close agreement of the 
celestial coordinates given in table 3, and from the reasonable agreement 


between the predicted and the observed intensities, the two radio sources _ 


may be tentatively identified with the clusters. 


§6. Large ScaLE [RREGULARITIES IN THE DISTRIBUTION OF 
Extra-GaLactic NEBULAE 


An attempt has also been made to detect by radio-frequency radiation 
the large scale irregularities in the distribution of extra-galactic nebulae. 
An examination of the distribution of the bright nebulae, as given by 
Shapley and Ames (1932) shows conspicuous grouping in high latitudes 
in the northern galactic hemisphere. Since this grouping is observed 
in a region where the effects of visual obscuration by the Galaxy are 
small, it appears likely that it must correspond to a genuine large scale 
clustering of the nebulae. It is therefore to be expected that these 
irregularities in the distribution of the nebulae will give rise to 
irregularities in the radio isophotes of the corresponding regions of the sky. 


(a) Analysis of the astronomical data on the distribution of nebulae 

For the purpose of a preliminary investigation a region of sky was 
_ chosen between declinations N 44° and N 55° and right ascensions 
10" and 184. This region contains marked grouping in the visual dis- 
tribution of nebulae, and it also lies largely in high galactic latitudes 
where the disturbing effects on the observations due to the radiation 
from the Galaxy are at a minimum. The main features of the distribution 
of nebulae in this region have been found by considering only the 
nebulae fainter than magnitude +13-0. In this preliminary analysis 
the effect of the bright nebulae has been omitted since they are so 
scattered that they must be considered individually for each position 
of the aerial beam. Their effect has, however, been taken into account 
when making the detailed investigation of the radiation received at 
declination N 47° 30’ as given in § 6(c). The analysis has been based 
on a detailed survey made by Shapley and Jones (1940) which gives 
the counts of nebulae on a series of photographic plates taken in a strip 
5 degrees wide in declination and centred on N 43°. These counts cover 
the nebulae in a magnitude range down to +17-5 and are divided into 
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groups at intervals of 0-1 magnitude. The results have been used to 
estimate the variation with right ascension of the radio intensity which 
would be received from these nebulae by a 2 degree beam as it sweeps 
round a line of constant declination centred in this strip. In making this 
estimate the intensity from the nebulae in the beam at every value of 
right ascension was found by weighting the contribution from the nebulae 
according to their photographic magnitudes, and by expressing the sum 
of their radiation in terms of an equivalent number of +13-0 magnitude 
nebulae. The number of nebulae in the beam at each position was found 
from the effective coverage of the beam and from the average density 
of nebulae on each photographic plate. The results of the analysis are 
given in fig. 1 and show clearly the irregularities in nebular distribution 
which can also be seen on the maps of bright nebulae. There are 
concentrations of nebulae to be seen at approximately R.A. 12h, R.A. 132 
and R.A. 16%, and the maps of the distribution of bright nebulae suggest 
that these concentrations may extend over a considerable range of 
declination at about the same values of right ascension. 


Fig. | 
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Calculated variation with right ascension of the radio intensity received from 
nebulae in a magnitude range +13-0 to +17-5. Declination N 43°. 
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_ The intensity of the corresponding irregularities to be expected in 

the radio isophotes was estimated by using: the same relation between 
the radio magnitude and the photographic magnitude as in § 5. 
It can be seen from fig. 1 that the maximum irregularity in the 
distribution of nebulae in the beam corresponds to the intensity from 
about four +-13-0 magnitude nebulae, and hence to a radio magnitude 
of +11-5. The limiting sensitivity of the present equipment corresponds 
to a radio magnitude of about +10, and. hence the maximum irregularity 
shown in fig. 1 will be 1-5 magnitudes below the limits of detection. mince, 
however, the survey extends only to nebulae of magnitude +1 v5, it 
is not possible to predict accurately the magnitude of the corresponding 
irregularities in the radio isophotes as the tendency to cluster clearly 
extends to fainter nebulae. However, the curve in fig. 1 suggests that, 
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if the total ‘contribution from the nebulae outside the range of magnitude 
-+13 to +17-5 is appreciable, then it might be possible to detect the 
irregularities in the radio isophotes with the present equipment. . 


(b) General comparison of the radio observations with the astronomical data 


A preliminary survey has been made of the radio radiation from the 
selected region. The survey consists of a series of curves showing the 
radio intensity versus right ascension observed at a number of different 
declinations. These observed curves all show small irregularities 
superimposed on a gradient of the ‘foreground’ radiation from the 
talaxy, which even at high galactic latitudes is considerable. In order 
to show these irregularities on a suitable scale, it has been necessary to 
subtract from the observed curves an.arbitrary smooth gradient. The 
exact shape of this gradient is not important since it has been taken as 
very smooth compared with the irregularities. The resultant curves are 
shown in fig. 2, and they show the presence of two broad maxima at 


Intensity 
(one division represents fifty 
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Irregularities in the background radiation observed at different declinations. 


R.A. 125 and R.A. 15" to 164, (The two sharp maxima at about R.A. 142 
and R.A. 13" are probably radio stars.) These broad maxima are of 
very low intensity, and to show them clearly above the random noise 
it has been necessary to average at least four records at each value of 
declination. It can be seen from fig. 2 that the maxima persists over 
a range of declination, and this suggests the presence of two extended 
bands of radiation. The correlation between the position of these bands 
and those predicted from the survey of nebulae suggests that the 
irregularities in the radio intensity are associated with the clouds of 
nebulae and that they do not represent effects due to the radiation from 
the Galaxy. 

The intensity of the irregularities is much greater than predicted in 


fig. I, and the discrepancy between the predicted and observed values 
varies from 3 to 4 magnitudes. 
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(c) Detailed comparison of the radio observations with the astronomical 
data for declination N 47° 30’ . 


A more detailed effort has been made to correlate the distribution of 
nebulae with the shape of the radio irregularities observed when the 
beam is directed to declination N 47° 30’. In this analysis the additional 
effect of the bright nebulae was taken into account by considering them 
individually. The radiation from each nebula’ was weighted according 
to its photographic magnitude and its position relative to the beam 
when centred on declination N 47° 30’. The sum of the radiation from the 
bright nebulae in the beam was then found for all values of right ascension, 
and was combined with the results given in fig. 1 to give the contribution 
from all nebulae down to magnitude +17-5. ‘The resultant curve is 
shown in fig. 3(a). This procedure assumes that the irregularities in the 
distribution of faint nebulae (magnitude +13-0 to +17-5) found for 
the strip between declinations N 41° to N 46° can be applied to a strip 
2 degrees wide centred on declination N 47° 30’. 

Fig. 3(6) shows the irregularities observed in the radio isophotes over 
a short period of right ascension in a strip centred on declination 
N 47° 30’. The curves in fig. 3 (a) and fig. 3 (6) both show a broad maximum 
at R.A. 12"; however, the effects of individual: bright nebulae, which 
are so marked in fig. 3 (a), are not clearly shown in the observed curve 
for fig. 3(b). It can also be seen that there is a discrepancy between the 
predicted and the observed intensity of about 3 magnitudes. This 
disagreement is not unexpected since the contribution from nebulae 
fainter than +17-5 has been omitted. It is, however, interesting to find 
that if the contribution from the fainter nebulae is increased by a factor 
of fifteen and is then added to the contribution from the bright nebulae, 
a much closer agreement in shape can be obtained. This is shown in 
fig. 3 (c). 

The sharp peak in the observed curve of fig. 3 (6) at R.A. 122 50™ is 
probably due to a radio star, but in view of the results presented in this 
paper the possibility that it is a cluster of nebulae cannot be excluded. 


§7. Discussion oF RESULTS 


The results given in §§4 and 5 have been examined for any relationship — 
between the apparent photographic and the apparent radio magnitudes 
of the nebulae observed. The available data is scanty and is summarized 
in table 4. In fig. 4 the difference between the two scales of apparent 
magnitude is plotted against the apparent photographic magnitude. 
The value given for a +16-0 magnitude nebula, derived from the Perseus 
cluster must be regarded as uncertain ; however, the results suggest that 
there is a relationship between the two magnitude scales, and that the 
apparent photographic magnitude increases more rapidly than the 
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Comparison of the calculated and observed irregularities in the background 
radiation for a strip centred on declination N 47° 30’. 


(a) Calculated curve for all nebulae down to magnitude +17-5. 
(6) Observed curve. 


(c) Calculated curve after increasing the contribution from nebulae 
fainter than magnitude +13-0 by a factor of fifteen. 
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A possible explanation of this discrepancy between the two scales 
might be found in the difference between the types of nebulae observed ; 
however, this seems to be unlikely. Both N.G.C. 4258 and N.G.C. 224 
are Sb type nebulae and yet they show the discrepancy in the two scales. 
Furthermore, the majority of the nebulae in the Perseus cluster are 
type E nebulae, and the comparatively high intensity observed from this. 
cluster would have to be explained by assuming that the ratio of radio 
flux to light flux is greater for type E than for types Sb or Sc. If this 
assumption were correct then the radio radiation from the bright elliptical 
companion (N.G.C. 205) of N.G.C. 224 would be comparable to, or greater - 
than, the radio radiation from N.G.C. 224 itself, and the measurements 
on N.G.C, 224 indicate that this is not so. It is clear that at present there 
are insufficient experimental results to establish the relationship between 
the scales of radio and photographic magnitude, but the present 
preliminary measurements suggest a divergence between the scales 
which requires further investigation. 

If the relationship between the photographic and the radio magnitudes 
suggested by fig. 4 holds for all nebulae, even approximately, then it is 
unlikely that any known individual extra-galactic nebula will be observed 
as a radio source more intense than N.G.C. 224. Many of the known radio 
sources are in fact much brighter than N.G.C. 224, and the absence of 
correlation between the positions of these sources and those of visual 
objects, together with their small apparent angular diameter (less than 
a few minutes of arc) supports the view that they cannot be identified 
with individual extra-galactic nebulae. This evidence, together with the 
observed concentration of the bright sources towards the galactic plane 
(Westerhout and Oort 1951), suggests that the brighter of the known 
sources are in the Galaxy. The results of a survey of the northern 
hemisphere (Ryle, Smith and Elsmore 1950) indicate that there are at 
least 50 radio sources brighter than radio magnitude +10 (as defined 
in table 1), and a survey of the bright external galaxies shows that in 
this region there are only 10 nebulae with an apparent photographic 
magnitude, and hence an apparent radio magnitude, less than +10. It 
therefore seems probable that less than 20% of the sources known at 
present are individual extra-galactic nebulae. 

The extension of the work to fainter individual extra-galactic nebulae 
presents considerable difficulties, as not only will it be difficult to detect 
such weak sources, especially in the presence of irregularities in the 
radiation from the Galaxy, but it will also be increasingly difficult to 
distinguish nebulae from radio stars. On the other hand, the detection 
of clusters of nebulae promises a method. of extending the radio 
observations to much fainter nebulae. It appears from the results in 
§ 5 that there are several clusters of nebulae whose radio radiation should 
be detectable as weak sources, and it may prove that these clusters can 
be distinguished from radio stars by their apparent angular diameter. 
A study of these clusters should yield information on the radio radiation 
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from faint nebulae, and it may also give data on the relative radio intensity 
from different types of nebula. 

The observations of the large scale irregularities in the distribution 
of extra-galactic nebulae are not yet, sufficiently advanced to warrant 
a detailed discussion. However, the preliminary results suggest that the 
analysis of extra-galactic radio-frequency radiation may provide a new 
tool for examining the metagalaxy. 

In a previous publication (Hanbury Brown and Hazard 1951) the 
authors estimated that the contribution of all extra-galactic nebulae to 
the radio-frequency radiation incident on the earth is about 1 % of the total 
radiation observed. This estimate was made by applying to all nebulae 
the ratio of radio flux to light flux measured for N.G.C. 224. The results 


? Fig. 4 


A 
3 
7 Apparent photographic 
magnitude (M,). 


13 15 


14 16517, 


! 
No becl 


photographic magnitude M,—M,, 
' 
W 


Apparent radio magnitude—apparent 


-4 
Comparison of apparent radio and photographic magnitudes of nebulae. 
The full line represents only the slope of the relationship between M, and 
M,—M,, and to give the same numerical values it should be displaced to 
pass through the point M)=9-7. 


shown in fig. 4 suggest that this ratio cannot be applied to all nebulae and 
that the magnitude of the extra-galactic radiation has been under- 
estimated. If the calculation is corrected by means of the data in fig. 4 
then the value of the extra-galactic contribution rises to approximately 
10%. This corresponds to an effective aerial temperature, averaged. over 
the whole sky, of about 100° K at 100 mc/s. Ina recent paper Westerhout 
and Oort (1951) find that the contours of radio-frequency radiation, 
observed over the whole sky, can be explained in terms of ete Tae from 
the Galaxy superimposed on an isotropic contribution of 600 xK at 
100 me/s. They suggest that this temperature may SO the een 
from extra-galactic nebulae and it is interesting to note that the results 
presented in this paper also suggest that there should be a significant 


extra-galactic contribution. 
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SUMMARY 


A theoretical investigation is made of the distribution of electrons 
round a positively charged impurity dissolved in a monovalent metal. 
Applications are made to dissolved hydrogen, where the impurity is a 
proton, and to atoms such as zinc, gallium, etc., which are usually 
considered to contribute their electrons to the conduction electrons. In 
all cases the positive charge must be screened; and in many cases this is 
shown to occur through the formation of bound states below the level of 
the Fermi distribution. The relation of these results to the Hume— 
Rothery rule is discussed. The ideas introduced are used to calculate the 
heats of solution of hydrogen and of polyvalent metals in the noble metals, 
and to discuss the magnetic properties of these alloys. A detailed dis- 
cussion is given of X-ray emission and absorption, the vacancy left in the 
x-ray shell being here treated as the positive impurity. In certain cases 
quantitative predictions are made about the energies of x-ray absorption 
edges. 

Discussions along the same lines are given of the optical absorption of 
the noble metals, and of the x-ray emission spectra of certain alloys. 


§1. [INTRODUCTION 


THE primary purpose of this paper is to study the distribution of electrons 
in metals such as copper, silver and gold in the neighbourhood of 
substitutional impurities of different valency such as Ni, Zn, Ga, Sn, As. 
The Hume-Rothery rule suggests that the last four of these metals lose 
respectively one, two, three and four electrons to the conduction band ; 
on the other hand the electrical resistance (Mott 1936) and other evidence 
“suggest that there must be a concentration of negative charge round the 
dissolved ion, so that a volume surrounding it is neutral. Similarly 
around Ni or around a vacant site there will be a concentration of positive 
charge. In this paper, we study the nature of this polarization, or con- 
centration of charge, and its effect on the binding energies and magnetic 
properties of alloys, and on the x-ray emission and absorption spectra of 


metals and alloys. 


* Communicated by Professor N. F. Mott, F.R.S. 
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The problem is treated throughout by the Hartree-Fock approximation, 
in which its own wave function is ascribed to each electron. The main 
problem of interest is the following : in the case of attractive potentials 
(e.g. those of Zn, Ga, etc. in Cu, or of a proton in a Fermi gas) do there or do 
there not exist occupied bound levels below the lowest level of the Fermi 
distribution ? It should be pointed out that this is not merely a question 
of finding the most accurate approximation for calculating energies, but is 
capable of direct experimental investigation by X-ray emission or magnetic 
measurements. 

So far, most of the calculations dealing with dilute alloys have been 
concerned with two atoms of the same valency. Thus Huang (1948) has 
dealt with Au in Ag and. Arafa (1949) with Cu in Ag and Ag in Cu. 
Huang’s calculations, none the less, show a polarization round the 
impurity. Huntington’s (1942) calculations for the energy of a vacant 
lattice point in copper deals with an impurity having negative charge (like 
nickel), so that the question of bound states does not arise. Mott (1936, 
Mott and Jones 1936 a) applies the Thomas Fermi method to the case of 
positively charged impurities ; but this method implicitly assumes the 
existence of bound states, and so begs the question. 

Our general conclusion is as follows. For an impurity in a monovalent 
metal carrying an excess charge of 2e, 4e etc. (Ga or As in Cu) there will 
always exist doubly occupied bound states below the level of the Fermi 
distribution. These are subtracted from the first Brillouin zone, so that 
the total number of electronic states remains the same as before ; Jones’s 
(1934) explanation of Hume—Rothery’s rule is thus not affected. The 
screening of the positive charge is entirely due to these electrons. The 
Fermi electrons are affected as follows : there will be an absence of charge 
at the centre of the impurity, compensated by a heaping up of charge in a 
region round it; the total displacement of charge must be zero, so no 
contribution is made to the screening. 

In the case of an excess charge unity, as for example Zn in Cu, an 
occupied level may or may not exist. Three configurations are possible, | 
depending on the relative energies of this level and of the conduction band. 
If no occupied level exists, screening is produced by a shift of the conduc- 
tion electrons, in the same way as for a negative charge (this is probably so 
for Zn or H in Cu etc.). If the bound level is occupied by one electron, this 
screens the charge, and the conduction electrons are polarized in such a way 
as to give no resultant shift of charge. Finally two electrons in the bound 
level would provide with the positive charge of the ion a net negative charge, 
which must be screened by pushing away the Fermi electrons. Configura- 
tions of the last two types appear in X-ray phenomena and probably when 
the alkali or alkaline earths metals are solvents of other metals. With an 
odd excess charge (e.g. Sn in Ag) similar results apply to the outermost 
electron of the dissolved monovalent ion (e.g. Sn*+). In other cases the 
additional electrons are always in bound states. 
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We use in this paper two different approaches. We first study the 
disturbances of a simple Fermi gas by a point charge (§§ 2to 4). This 
illustrates well the principles involved. The study of bound states in 
real metals, our main purpose, is best made by estimating their energy 
relative to that of the Fermi electrons (§ 6). The approximations made 
are justified by the analysis of a simple case, hydrogen in copper (§ 7). 
Applications are made to energies of solution (§ 9), magnetism (§ 10) 
and X-ray spectra (§ 11). 


§2. BEHAVIOUR OF AN ELECTRON IN A GIVEN FreLp or ForcE 


Before studying the conduction band of a metal, we must obtain some 
formulae about the motion of an electron in a spherical field of force : 
conditions for the existence of bound states and the displacement of 
charge towards the centre will be investigated. 

Let then an electron be enclosed in a sphere of large radius R and 
move in a spherically symmetrical field of potential V(r). If V(r) tends 
to zero more rapidly than r~? at infinity,* the possible electronic wave 
functions have radial parts with asymptotic form (see for instance Mott 
and Massey 1949 a) : 


*(r) ~r—? sin (k)'r-+-n,— dlr) (=0.1; 2 Re) 


If we make the wave function vanish at the boundaries of the (large) 
sphere, we find the following quantized values of the momentum, valid 
for large R : 

ky ~ (A+31)n/R—n|R 


where A is an integer. If k, is the value of k,;) when V vanishes : 
hy ~ (ALM) n|R 
and ki —hy ~ —n(k;')/R. . . . . . . ° (4) 


There is thus a one to one correspondence between quantized momenta 
k, and k,, and the difference between them depends only on the phase. 

We discuss now the dependence of 7 on k. 7 tends to zero for large k ; 
it will tend to one of the values Az as k tends to zero. A is here an integer 
or zero: it is zero for repulsive or small attractive potentials, and 
increases with the strength of attractive potential V. This behaviour is. 
shown in fig. 1. For limiting values of V just on the borderline between 
two values of A, 7) will tend to one of the values }7, $7, .... as shown by 
the dotted line in fig. 1. The other phases 7, do not behave in this way 
(cf. Mott and Massey 1949 c). 
ee i a 

* i r) varying more slowly than r-?, and as is well known for 
GRE Mcliy (Mote rat eee 1949 b), m» cannot be defined for zero 
momentum. As seen later, this means for attractive potentials that there is an 
infinity of bound states which the Fermi electrons could occupy. 


MZ 
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We study now the origin of the bound states. Let us assume for 
instance that there is one bound state for the potential V, hence that 
n(0)=7 (curve 1, fig. 1). Then, according to eqn. (1), the states 
kyp=7/R, 27/R .. . . etc., correspond to the various states ko’ in the presence 
of the field, k)=7/R corresponding to the lowest one k)’=O (fig. 2). 
Only the state k,=0 cannot be related to any real value of ky’, but is 
replaced by the bound state, which may thus be considered as subtracted 
from the bottom of the Fermi distribution. It is easily verified that this 
result is general: the bound states for a particle in the potential V are — 
always subtracted in that way from the lowest states in the Fermi 
distribution. 


Fig. 1 


Xm 


4 ; 
z at. units 


Phase shifts 7 (x) for attractive potentials (after Morse). a denotes 2kq-1. 
1:B=1-0,1=0; 2: 8=0-7, 1=0; 3: B=1-0, 1=1; 4:8=0-7, l=1. 


We now suppose that there are N electrons in the sphere of radius R 
and assume that V(r) is the self-consistent potential energy of each 
electron in the field of some impurity and of the Fermi gas thus formed. 
Let ky be the maximum momentum before the impurity is introduced. 
We compute the electronic charge which must be added, when the 
impurity is present, to obtain the same maximum ky. 

According to eqn. (1), /y, corresponds in each component | to a different 
value ky,,’ of the momentum : 


key — ky ~ —R-(ky’). PRE ena gaa (| RS 
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ae ae therefore add electrons to occupy states in the interval between 
m and ky (fig. 2). The number of such states, multiplying (1 a) by a 


density factor, is 
Rr-}(214- 1) (yp ky.) ~ — (20-4 1) *y)(kyp) 


where P Aahtty eh 
oe aN replaced n(hyy’) by (ky) because 7, does not vary rapidly 


Summing over / and over both spins o, the total charge introduced is 
therefore : 


Dome ETL Mp). Siar Moms eae Oh 2) 


Fig. 2 


Displacement of states with s symmetry in a box of radius R, when an attractive 
potential V is introduced (a) without V ; (b) with V. 


The quantity n can now easily be shown to be equal to the total excess 
charge locally attracted by the impurity. For let a be the radius where 
the potential V vanishes (fig. 2). It is a well known result of any 
consideration of normal modes of a vibrating medium, that, if ky(R—a) 
is large, the electronic density in the spherical shell a<r<R is independent 
of the boundary conditions at r=a and r=R, and the same as in an 
unperturbed Fermi gas (von Laue 1914). The total charge below ky in 
the spherical shell is thus not altered when the impurity is introduced. 
The charge added, n, must therfore be localized in the sphere of radius @ 
where the potential V is different from zero. 

In eqn. (2), the bound states are counted as occupied. _ A differentiation 
gives 

—n An (k)/ok] dk 
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for the charge displaced by potential V in the /th component, for values 
of k in range dk. In other words, the attraction of charge is obtained 
by a compression of the states in momentum space. Appendix (1) gives 
another proof of this result. 

The locally displaced charge n must compensate exactly the charge 
introduced by the impurity, so as to give no Coulombian term in V(r). 
Usually, the impurities introduced in the physical problem are electrically 
neutral (H, Zn. . . atoms; vacancies . . .): if the centre is charged, electrons 
are introduced too. Thus the charge introduced in the gas is just equal 
to n, and the Fermi gas conserves the same maximum momentum ky. 

In other words, a neutral impurity introduced into a Fermi gas modifies 
the density of states in the k space in such a way as to keep constant the 
top of the Fermi band, at least for small concentrations where impurities 
do not interact. This point is of importance in computing the energies of 
solution (cf. § 7, 9). 

Consider now as an example the screened Coulomb fields* 


V=—Zr— exp, (—¢r) i 22. 2 


which will be used in the next section. Z and q are two constants, and the 
corresponding phase shifts 7, were computed numerically by Morse (1932) 
for positive values of Z. We reproduce in fig. 1 three examples of varia- 
tions of 7(k). Morse’s notation is 


Ca kee SG WA] prea 


Contributions by p and higher components are negligible up to large 
momenta k, and the charge attracted is practically of pure s character. 
For small attractive potentials (e.g. 8B = 0-7), n) presents a maximum. Thus 
the density of attracted charge decreases with increasing momentum and 
some charge is even repelled for large momenta if the energy of the Fermi 
gas has a sufficient width. 

For stronger attractive potentials (e.g. B=1-0), a bound state appears. 
For all values of & and especially for small values, there is then a deficit of 
charge near the origin. The same occurs in.the limiting case where 
(0) = 37 (B ~ 0-9). mH has then a negative slope not only for potentials of 
form (3) but also for square potential walls of any dimensions,+ hence 
probably for any kind of purely attractive potential V. This repulsion 
will be used in the next section. 

In conclusion, we studied in this ‘section the displacement of the 
quantized levels in an enclosed Fermi gas, when a potential V(r) is intro- 
duced. The bound states appearing in an attractive potential are 
subtracted from the bottom of the Fermi band, and the charge locally 
attracted in the Fermi gas is simply related to the phase shifts of the 
quantized levels. 
ee ES oe ee 

* Atomic units used throughout. 


{ For a square well of radius a, the slope of at the origin is —la 
independently of the depth of the well. th "a 
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1 
§3. SELF Consistent SCREENING OF ADDITIONAL CHARGES BY FERMI 
ELECTRONS. 


The discussion so far applies to a gas of electrons moving in an arbitrary 
field of potential V(r). We shall study now the form of this field when 
it is due to an impurity carrying a point charge Z (e.g. Z=-+1 for a proton). 
To the Coulomb term Zr—! we must add a term due to the displacement of 
the electrons, obtained by a self consistent method. The displaced 
charge must exactly screen the charge Z ; it may be provided either by 
bound electrons having negative energy below the bottom of the Fermi 
distribution, or by electrons of positive energy and thus by the Fermi 
electrons themselves. We restrict ourselves here to the study of screening 
by Fermi electrons and show that for the simple Fermi gas, it occurs 
round negative charges Z only. 

Using the Thomas—Fermi method, Mott (1936) obtained for the self 
consistent potential V round a charge Z in a Fermi gas the screened 
Coulomb field given above by formula (3). ¢ is a constant nearly equal to 
1 atomic unit. For copper for instance, Mott’s computation gives 
q=0-98. Isenberg’s Thomas—Fermi method (1950) corrected for changes 
in exchange and correlation energies gives a formula from which we have 
calculated that q ~ 1:3. 

The phase shifts y, have been shown by Morse (loc. cit.) to depend very 
little on the precise form of the potential V ; and it is certainly sufficient 
for our purpose to assume for V a form of type (3) where q¢ is determined by 
the condition that the screening must be complete, and thus by the 
equation 

Qar-1D'(21-+- 1) (Ky) =Z. 
1 


For a positive charge Z (e.g. a proton) fig. 1 shows that contributions to 
the displaced charge by p, d... components are negligible in the range of 
momenta occupied by the Fermi gas. For copper for instance, ky=0-725, 
hence %,—=1 to 1-5. We must therefore have 


no(ky) = 372. 


The corresponding points in fig. 1 (A for Z=1, B for Z=2 etc.) are all na 
region where at least one bound state exists. Hence a simple Fermi gas 
cannot provide stable screening round a positive charge through unbound 
electrons only. This point will be confirmed in § 6, by the study of bound 
electrons. We must point out however that point A is not very far from 
the limiting curve corresponding to 7(0)=37 ; for actual metals (§ 6) or 
for finite concentrations of impurities, stable screening by Fermi electrons 
round impurities with charge +1 may be possible. 

For negative charges Z, no bound states are of course possible and a 
value of g giving screening by the Fermi electrons is always possible. 

We shall now extend these results on screening by Fermi electrons to 
actual metals. Configurations with bound states will be treated later on by 


a different method. 
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§4, SCREENING BY Fermi ELectrons. EXTENSION TO ACTUAL METALS 


As was shown by several authors (Wannier 1937, Tibbs 1939, Peckar 
1946, Mott and Gurney 1948, Slater 1949), the Wigner-Seitz approxl- 
mation can be applied approximately to the case where an impurity is im- 
bedded in a regular lattice. Thus in a monovalent metal with an impurity 
atom the conduction electrons may be considered as a Fermi gas of 
electrons with a certain effective mass (the same as in the pure metal), 
scattered by the perturbing potential introduced by the impurity, each 
wave function being modulated by the wave function ¢,, corresponding 
to the bottom of the band in the pure metal.* 

But these results are only approximate. For an electron of momentum k, 
there is a corrective potential in the Schrodinger equation, of the form 
(Vd, . Vbo)(b.¢2) 4, if dy exp (ik .r) is the corresponding wave function 
in the perturbing potential only. This is small for the screened Coulomb 
fields acting on the conduction electrons; but it is too large for 
unscreened fields acting on the bound states to give much more than a 
qualitative description. . 

Thus, owing to the modulation by the function ¢,, there appear below 
the 4s band of copper for instance and in an attractive potential of 
increasing strength, first a 1s, then a 2s, a 2p,.... states in an interstitial 
position ; or a 4s,.... state in a substitutional one.t The displacement 
of charge still follows eqn. (2) near the bottom of the conduction band. 
The general results on screening obtained for a Fermi gas must therefore 
still be true for monovalent metals, except perhaps in the limiting case 
Z=+1. This will be confirmed in the more complete picture of the 
bound states which will be obtained in §5 by another method, based . 
on energetic considerations. 


§5. Impurities with Bounp ELECTRONS IN METALS 


We consider here general characteristics of bound electrons round 
impurities in metals. Simple applications will follow in the next section, 
and in §7 a more exact study of a simple case (hydrogen in copper) will 
justify the approximations made. 

The method used so far in this paper has been to study the phase 
shifts , of the conduction electrons and deduce the number of bound 
states from the behaviour of 7, at the bottom of the conduction band. 
For actual metals, due to the periodic field and other corrections, this is 
difficult and we shall now obtain more reliable results by considering the 
energies of the bound states and comparing them with the range of 
energies of the conduction band. 

Sn ne 

* cf. Appendix (2) for a simple proof. 


} The next state is a 4p state, not a 5s one as one should expect if the 


Wigner—Seitz approximation was still valid for those strong perturbing 
potentials. 
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We must first define an energy of reference. We take for our zero 
as is usual, the energy of an electron at rest outside the metal. Now, 
in contrast with the Fermi gas studied in §§2 and 3, the bottom of the 
conduction band in actual metals has a large negative energy —Z). 
Ey may be deduced from the energy of sublimation E,, the first ponent 
of ionization #, and the Fermi energy H,(=0-6 H,,) of the metal 
(Mott and Jones 1936 b) by the equation : 


Ky > H+ HE, +E. 


Table 1 gives values of #y calculated from this formula. Values of E 
are taken from Seitz (1940 a) and Kubaschewski (1951), and H, is alesn 
from the Handbook of Chemistry and Physics (1947). . Corrective terms 
for exchange and correlation energies may be neglected if E, is calculated 
by this formula, and #, may be taken as in a Fermi gas. H, is about 15 ev 
for noble metals and varies from 10 to 6 ev for the Ia group (Li to Cs). 


Table 1. Energy Limit £, of the Conduction Band in Monovalent Metals 
(in ev) 


The main difference between the bound states of an impurity in 
insulators and in metals is therefore the following: in insulators the 
empty conduction band has a small electron affinity so that bound states 
have small negative energies. The radius of their orbits is large and 
they overlap a number of neighbouring atoms. The attractive lattice 
potential acting on them lowers the energy of their series limit down to 
the bottom of the conduction band, in accord with Wannier and Slater’s 
theories (ef. § 4); on the other hand the polarization of the overlapped 
neighbouring atoms reduces the ionization energy of the bound states 
by a factor K-? (K=dielectric constant, cf. Mott and Gurney loc. cit.). 
In metals on the contrary, impurity electrons have energies below —L,, 
hence their radii are small and in general do not appreciably overlap the 
neighbouring atoms. The energy of a bound state must therefore be 
related to the absolute zero of energy, not the bottom of the conduction 
band. The influence of the lattice may be treated asa small perturbation 
and the polarization of neighbouring ions neglected (cf. §§ 7 and 11). 
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In general the charge displaced in the conduction band does not much 
modify the form of the bound states, as will be confirmed. in §7 in the 
special case of hydrogen in copper. Ascribing arbitrarily the interaction 
energy between bound and Fermi electrons to the latter, we may therefore 
consider the bound states of an impurity as those of an electron moving 
in the field of the impurity only; then consider the conduction electrons 
as perturbe. by the resulting ion or atom. We may thus use for these 
bound states the same ionization energies H, as in vacuo. If H; > LZ), 
the bound state will certainly exist. If H,<) it is very likely to be 
replaced by some sort of screening by Fermi electrons. 

We must therefore compare H, with the ionization potentials H,; of 
the various elements. 


§6. EXAMPLES OF IMPURITIES WITH BouND ELECTRONS IN MONOVALENT 
METALS 


We study first in some detail the simplest possible case, an interstitial 
proton. 
6.1. Interstitial Proton 

A proton in vacuo may trap one or even two Is electrons, forming 
respectively a H atom and a negative H~ ion, the corresponding electronic 
energies being —13-5 and —0-7 ev (Massey 1950). Hence in the 
metal three configurations are a priori possible, depending on the value of 
E, for the solvent metal. The limits given are of course only a rough 
order of magnitude, owing to the approximations made; energies are in ev: 

(a) EH)>13-5: there will then be screening by Fermi electrons, giving 
an extra charge in the conduction band. 

(b) 0-7<H)<13-5: there will then be one bound electron, forming 
with the proton a hydrogen atom. The field of this atom polarizes the 
conduction band so as to give no total displacement of charge in the band. 


(c) Hj)<0-7: two bound electrons. The hydrogen negative ion thus 


formed repels one electronic charge in the conduction band. 

The simple Fermi gas considered in §§ 2 and 3 corresponds to the case 
H,=0. We see that in this case, in accord with our previous conclusions, 
there should be two bound electrons in a 1s state. The negative exchange 
energy between these bound electrons and the conduction band may 
increase the limit 0-7 ev of this case by several electron volts. But for 
actual monovalent metals the values of H, listed in table 1 are much 
larger than this and therefore configurations (b) or (c) are the only ones 
possible. Configurations of type (b) are probable in Cs, Rb, K, Na which 
thus absorb hydrogen, if at all, in an atomic form, at least for small 
concentrations.* Noble metals form a limiting case between configurations 


(a) and (b). In §7 a detailed discussion is given of this problem for the 


case of copper. 


_* The ionic hydrides stable at higher concentrations are probably configura- 
tions of type (c) stabilized by the Madelung energy. 


, 
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6.2. Substitutional Impurities 


Table 2 gives the successive ionization energies H; of a number of 
elements, H_ being the ionization energy of the negative ion. We see 
from these values and the values shown in table 1 that dissolved atoms 
in general lose their first electron to the conduction band of the matrix. 
All other electrons are retained in bound states, except perhaps for the 
second electron of some di- and tetravalent metals when the matrix is a 
noble metal. Thus the second ionization energy of Ga is 20-5 ev and 
the energy of the electron is well below —H, (=—15-4 ev) for Cu. 
But for Zn, the difference is smaller, and magnetic properties seem to 
indicate that the ion Zn*++ is screened by Fermi electrons (cf. § 10). 

Some metalloid elements dissolved in metals to the left of them in the 
periodic table probably take from the matrix one or more electrons to fill 
up their np shell, these electrons being in bound states. We have thus a 
configuration of the type (c) (H—-) studied above. The range of the primary 
solution is however very limited. At higher concentrations the diamond 
structure with homopolar binding (type Na, P) appearing in compounds of 
alkali metals with P, As, Sb and Bi may be considered as a configuration 
of this type, where the valency electrons of the metallic atoms just fill the 
zone derived from the s and p electrons of the metalloid element (Brauer 
and Zintl 1937). The same is true for ionic crystals (NaCl). In both 
cases, the configuration of type (c) is stabilized by the Madelung energy. 


Table 2. Some Examples of Ionization Potentials H,. H_isfor Negative lons 


Be Mg Zn Cd Al Ga Si Sn NP] As 45.0 


_ | (0-2) |(—0-9) (0-2) (1-7) (1-0) 23 | 38 
, | 9:28 | 7-61 | 9-36 | 8-96 | 5-96 | 5-97 | 8-12 | 7:30 [14-48 | 10 [13-55 | 12-95 
18-12 | 14-96 |17-86 | 16-84 | 18-74 | 20-43 |16-27 |14-5 | 29-47 | 20-1] 34-93 | 23-67 


§7. SOLUTION OF HyDROGEN IN COPPER 

The purpose of this section is to study whether the charge of a proton 
dissolved in copper is screened by a displacement of the Fermi electrons, or 
by an'electron ina bound state. A self consistent solution is computed for 
the bound screening in the case of a Fermi gas. Corrective terms are 
then added to the energy of solution —£,, which is compared with 
experiment and with the maximum value of HZ, one would expect for 
screening by Fermi electrons. 

Bound screening gives a value for H, in satisfactory agreement with 
experiment. But Fermi screening perhaps occurs, for the corresponding 
maximum value for H, seems a little larger than the experimental value. 
We are clearly in a limiting case. Results obtained for bound screening 
are however certainly valid for other phenomena, e.g. round a proton or a 
divalent substitutional impurity in alkali metals and in X-ray phenomena 
(cf. §§ 8, 10, 11). 
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7.1 Proton with a Bound Electron in a Simple Ferma Gas 

The proton and the bound electron form a hydrogen atom which 
polarizes the gas. The process of successive approximations was therefore 
started with a normal hydrogen atom and plane waves for the gas. 
Finally we obtained the wave function of the bound electron in the field of 
the proton and of the charge displaced in the Fermi gas ; and the wave 
function of the Fermi electrons in the field of the proton, of the bound 
electron and of the displaced charge. The solution was judged sufficiently 
self consistent when the total charge displaced in the Fermi gas was reduced 
to 0-1 electronic charge. 

We had to compute a number of rapidly oscillating wave functions with 
s symmetry, for which the usual methods are not very appropriate. We 
devised therefore the following method. 

Decomposition of Schrédinger’s radial equation. The Schrédinger 
equation 
(d2bjdr*)--4 K (r) b= O02 Wat, ee 
may be split into two differential equations of the first order. By analogy 
with the case where K is constant, we write 


p=a(r) sin u(r), with = (1/s) (0b/0r)=K cot u. 
It is easily shown that (4) is then equivalent to the system 
Cuno Ke 1 O@awee OA 


2K or 2D 2u+KkK ; am BK or (°° 2u+1) sae Oo) 
and that for an electron with momentum k one has 
ou oe ‘ Dhe ft ' ‘ 
FY aaOy ‘ee ut ag | dr. <* ) Rees 


For an s wave function, K? is of the order of r~! at the origin, hence 
u~Kr. The eqn. (5) in w may thus be solved by successive 
approximations, starting with w=Kr. This gives wu at infinity by a process 
in general rapidly convergent. The amplitude a is then obtained by a 
simple integration. 

For large r and if the potential V=}(K?—k?) tends to zero rapidly 
enough at infinity, 

u~kr+y 
where 7 is the phase shift. A development of a for large r shows that the 
conditions are: r"V->+0 for r+0o, with n=2 if k=0 and n—1 if k=£0*% 
and = Ar if k=0 (A integer). 

It is worth noting that eqn. (6) leads, for large r, to eqn. (2) when 
we replace wu by kr+n, a by 1 and K by k. For k0 on the other hand, 
system (5) tends for slowly varying potentials to the WKB approximation. 

We shall now review the main points of interest concerning exchange 
phase shifts and the energy of solution. 

Huchange. The main complication arises from the exchange inter- 
action with the bound electron which introduces a distinction between 
the spins of the Fermi electrons. Fig. 3 gives the potentials acting on 
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them. The dotted line (curve III) represents the Coulomb potential 
provided by the proton and the bound electron. The polarization of the 
Fermi electrons gives an additional barrier ; the potential IT, acting on the 
Fermi electrons with spins anti-parallel to that of the bound electron 
becomes therefore slightly negative at large distances. This barrier 
is masked by the attraction due to exchange in potential I acting on 
electrons with parallel spins. 

The average exchange term AV, due to the bound electron was 
computed using Slater’s approximate formula (1951) 


gees 3 AZ \iu8 
Av——;| %,+3(ssa S) i: ns op eT) 


Fig. 3 


rove) 
? at. units 


=0-5 


Potentials acting on Fermi electrons with spin parallel (I) and antiparallel (IT) 
to that of the bound electron. The punctuated curve III gives the 
Coulomb potential of a hydrogen atom. Curve IV gives the comple- 
mentary potential acting on the bound electron in Slater’s method. 


where Z, is the ‘effective charge ’ in Hartree’s sense. Correlation 
between antiparallel spins and changes in exchange interaction between 
Fermi electrons were neglected in this computation, which does not claim 
to be very precise. 

AV, is somewhat different from the exchange potentials one can 
calculate for various k and J in the last stage of the computations (fig. 4). 
Owing to the approximations made on other terms, formula (7) gives 
nevertheless a satisfactory average value. It has the advantage of giving 
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for the Fermi electrons of parallel spins a potential I which only differs 
by a nearly constant term from the potential acting on the bound state 
(fig. 3, curve IV). The bound state has therefore almost the same form 
as if it were moving in potential I, as suggested by Slater (1951). Its 
wave function is practically orthogonal to those of the Fermi electrons 
with parallel spins, and this simplifies the computation of exchange 


energies. 


Fig. 4 


» at. units 


Exchange potentials obtained in the last stage of computations for various 
land k, compared with Slater’s average value. 


a,: l=0, k=0-725, a,? 1=0, L=0-5; ag: 1=0, K=0:3; b,: 1=1, k—0-725, 
b,: l=1,h=—0-bS Ci lee to 


Phase Shifts. Fig. 5 and table 3 give the phase shifts and the displaced 
charges (2/-+1)m~'y(ky) for electrons of both parallel (I) and anti- 
parallel (II) spins in the Fermi gas. Potential I is purely attractive ; 
and all the spherical components with parallel spin are attracted except 
the s component; this last is repelled owing to the presence of the 
bound Is electron (cf. § 2). Owing to the mixed character of potential IT, 
attractive at short distances only, the s component with antiparallel spin 
is slightly attracted at the top of the conduction band, but the p, d 


9.) eer one 


Hlectrons Round Impurities in Monovalent Metals 167 


components repelled. The total displaced char i i 
ge gives rise to the Ut 
Coulomb potential barrier described above. : a 
Fig. 5 shows that in the part with parallel spin 


1(0)=7. 
One state, corresponding to the bound electron, is subtracted from the 
bottom of the Fermi gas. This still holds, as a computation showed 
when the average potential AV, is replaced by the more exact alee 
obtained for the s component at the last stage of the computations. This 
point will be used in the study of Hume—Rothery’s rules (§ 8). In table 3, 


Fig. 5 


Phase shifts 7,(&) at infinity for electrons with spin. parallel (I) and anti- 
parallel (II) to that of the bound state. 


Table 3. Charge Displaced in Spherical Component / for Electrons of 
Spins Parallel (I) and Antiparallel (II) to that of the Bound Electron. 
The Bound Electron gives a Charge 1 counted separately 


i ee ie aa Pe rs cen oe 


IT |1—0-241) 0-450 0:340 | 0-238 0-170 0-086 0-027 1+1-07 
Il | 0-022 |—0-384 |—0-355 |—0-208 |—0-095 |—0-034 —0-016 |—1-07 
I-+I1I |1—0-219} 0-066 |—0-015 0-030 0-075 0-052 0-011 1+0 


Peeper ee ee 
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the displaced charges have been slightly corrected so as to give no extra 
charge in the Fermi gas. We see that there is about: one additional 
electronic charge of parallel spin and a corresponding missing charge of 
antiparallel spin. 

Energy of Solution. For the energy of solution of a hydrogen atom, we 
computed changes in Coulomb and exchange energies directly, by the 
usual methods (cf. Fock 1933). We thus have : 

(a) The change £ in Coulomb interaction between proton and bound 
electron. 

(b) The change HZ’ in Coulomb energy for the Fermi electrons. They 
are in the field of the proton, the bound electron and the charge displaced 
in the Fermi gas. This gives a Coulomb energy £,’ from which one must 
deduce the self Coulomb energy FH,’ of the displaced charge, counted 
twice in H,’. 

(c) Finally the exchange energy between bound electron and Fermi 
electrons of parallel spins, computed for s, p, d components (Hy", L,", £,") 
and evaluated roughly for the rest (£,”). 

We find, in atomic units* 


Ee -2 0-135 
E,'= 0-068 
—E,'=—0-093; 
E’ =~—0-025," 
Ey" =—0-025 
E," =—0-080 
E,." =—0-030 
E;" =—0-015 
EH” =—0-140 


Hence for the energy of solution per hydrogen atom we have 
E,=E+E'+E"=—0-030=—0-8, ev. 
, is therefore small and mainly due to the difference between the exchange 
energy &” and the change in Coulomb energy H of the bound state. The 
change in exchange energy between Fermi electrons has been neglected 
here. It is certainly small. For, by analogy with the results obtained 
for a vacancy in copper (Huntington loc. cit.), it is the difference of two 
terms of the order of 0-03 atomic units corresponding to the charges 
attracted and repelled of respectively parallel and antiparallel spins. 
7.2 Hnergy of Solution. Corrective Terms due to the Copper Ions 


When hydrogen is dissolved in copper, and not in a simple Fermi 
gas, several other terms must be considered. 


*1 Atomic unit=2 Rydberg units=27-06 ey. 
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We have first to consider the corrective potential (44, . 44,)(¢, ¢.)~4 
in which ¢, is the wave function at the bottom of the conduction band. 
It is constant over a large region round an interstitial hydrogen atom 
(cf. Fuchs 1935), and the corrective potential may therefore be neglected. 
Also the electronic density is nearly constant throughout the metal. 
The approximation of a simple Fermi gas seems therefore satisfactory for 
the Fermi electrons.* 

The bound electron however interacts in two ways with the copper ions. 
Owing to the exclusion principle, it cannot enter the regions occupied by 
the 1s’ shells of the copper ions. These regions are practically restricted to 
small spheres and scatter strongly the wave function of the bound state. 
When the proton is at the common vertex of six atomic polyhedra, it 
seems probable that the 1s? shells are sufficiently far from the proton and 
close enough together practically to stop the bound electron and thus to 
enclose it in a spherical volume of radius ~ 3-2 atomic units. A calculation 
by Sommerfeld and Walker (1938) enables us to estimate that the energy 
of the bound state is thus raised by 

; E£,=0-05 Atomic units 
without appreciable change in the electronic densities. 

On the other hand, the attractive potential wells of the copper ions 
Cu* are neutralized by the charge of the conduction band at the surface of 
the atomic polyhedra only. Assuming that the bound electron is enclosed 
in the sphere described above but is otherwise unperturbed we compute 
that its potential energy is lowered by 

H,=— 0-095 Atomic units. 
7.3 Energy of Solution. Comparison with Expervment and with the 

Energy for Screening by Fermi Electrons 

The energy of solution H, for screening by a bound electron is obtained 
by summing up the various changes in energy ; 

E,=H,+#H,+H,=—0-07, Atomic units=— 2-0; ev. 
The experimental value is deduced from measurements of the solubility of 
H, and from the dissociation energy of the hydrogen molecule (Fowler and 
Smithells 1937). It is 
E' = —0-060 Atomic units=— 1-6; ev. 
- Owing to the approximations made, the agreement may be considered as 
satisfactory. ee: 

Without computing in detail the self-consistent solution for screening by 
Fermi electrons, we may compare with H, the minimum energy of a 
configuration of this type: i.e. when the screening charge is just at the 
bottom of the conduction band (cf. § 5). ‘The Coulomb terms are 

Ey—EH,—E,. 
Here H,=15-4 ev is the energy of the bottom of the conduction band 
(cf. table 1); Hy=13-5 ev is the electron affinity of the proton, and #, is 


eee ee ee eee 
* An effective mass m*41 should perhaps be used (cf. Mott and Jones 1936 c). 
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the self energy of the displaced charge. To evaluate this small term, we 
assume that the displaced charge has a density p(r) such as to give Mott’s 
approximate screened coulomb potential (formula 3) : 


Var3[ pdr +] pr-! dr=Zr-"[1— exp(— qr) ]. 
0 r 


This gives a density p very similar to that of a hydrogen atom. And we 
obtain 
B,=| pV dr=0-13 ev. 
0, 
Hence the total of the Coulomb terms is 
— 2-03 ev.* 


The exchange term E’’ must be of the same order of magnitude as for a 
vacancy, but with opposite sign (Huntington loc. cit.). Hence using 
Huntington’s results ; 
E" ~—0-03 Atomic units=—0-81 ey. 
The corresponding energy of solution (the sum of #’’ and the Coulomb 
terms) is 
EH" =— 2:8; ev. 

Owing to the approximations in the calculation, the exact value of this 
term is devoid of significance. It merely shows that a screening by Fermi 
electrons is more probable than by bound electrons, the additional charge 
being near the bottom of the conduction.band. Conversely, in a screening 
by bound electron, if it exists, the electron must be very close below the 
band. Hydrogen in copper is clearly a limiting case. 


§8. Hume—Rotuery’s RULE 


We have now to reconcile our description of valency electrons of 
impurities in bound states with Hume—Rothery and Jones’s point of view 
(Jones 1934), according to which all the valency electrons are contained in 
the conduction band. We show how this can be done and describe a 
mechanism which may explain some experimental deviations from the 
rules. 

In screening by Fermi electrons the impurity merely displaces states in 
the zone of the conduction band, without creating new ones (cf. § 2).+ 
One more state is occupied for every attracted electronic charge, in accord 
with the rule. 

If a proton keeps a bound electron, it subtracts one state from the 
bottom of the conduction band (cf. §7). Similarly a He atom subtracts 
a state of two electrons (§ 3, fig. 1, point B). Thus every bound electron 


* The width of the Fermi band is constant (§2); its variation gives no 
contribution to the energy of solution. 
nt This was proved for a Fermi gas. Wannier’s approximate theorem (§ 3) 
allows an extension to a general zone structure. No direct general proof has 
been given as yet. 
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corresponds to an electronic state subtracted from the band and hence 
may be said to belong to the band. For this description to be rigorous, 
the bound and Fermi states should be orthogonal. We saw however in §7 
that owing to exchange this condition is almost satisfied for hydrogen in 
the Hartree-Fock description. 

These results may be directly extended to. substitutional impurities 
of the same row as the matrix. Thus Zn, if when dissolved in copper it 
has one 4s bound electron, probably subtracts one electron from the 
4s conduction band of the matrix; Ga, two, etc. As the valency 
of the impurity increases (Ge, As...), the exchange energy of the 
Fermi electrons with the smaller s shell decreases, and it is not a priore 
certain that the states subtracted are not fewer than the additional 
valency electrons. Examples of diminution of valency (Elektronen- 
ricktritt, Umladung) are in fact cited by Schubert (1948). Thus for 
example Sb in Cu forms at high temperature a body-centred phase 
around the Cu,Sb composition, which is readily explained as a Hume— 
Rothery 6 phase if Sb contributes three valency electrons only to the 
conduction band (Hofmann 1941). 

_ Let us now consider two elements of different rows in the periodic table, 
Ag in Cu for instance. Ag has a full bound L shell (4s? 4p6 4d!°) which 
certainly corresponds to states subtracted from the conduction band of Cu. 
In other terms the perturbing potential introduced by Ag is attractive 
enough to give phases shifts 7, 7,, 72 tending to 7 at the bottom of the 
conduction band (cf. §2). But the same also applies with M states at the top 
of the same zone, where the same perturbing potential acts on ‘ positive 
holes’ instead. of electrons, the corrective potential being also slightly 
different (cf. §4). The zone corresponding to the conduction band, therefore, 
gets rid of the L states at its bottom and accepts M states instead at its 
top. Ag may be said to modulate the 4s band into a 5s one round itself 
without altering the number of its states. Similarly, Cd, In, Sn... . act 
as Zn, Ga, Ge... . in Cu, but for the bound L shell introduced below the 
conduction band without changing the number of states. 

Summarizing, the valency electrons of the impurities are in general 
in the conduction band or in bound states subtracted from the band ; 
they are the only ones to have this property. They may indeed be 
considered as a mere extension of the band. 

We may note in conclusion that the configuration corresponding to a 
_ negative ion H~ does not fit with this description. In the Hartree-Fock 
approximation, the conduction electrons move in a repulsive potential 
which does not subtract any state from the bottom of the conduction 
band (cf. § 6.1). This approximation is probably satisfactory. For the 
s part of the conduction band is repelled at the bottom of the band and is 
replaced by the p component at the top, owing to overlap of the s and p 
bands (Krutter 1935). The 1s bound state is therefore nearly orthogonal 
to it and mainly subtracted not from the conduction band but from the 
rest of the system of s bands. Such a configuration requires the formation 


N2 
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of a repulsive ion, i.e. at least a negative interstitial ion or a neutral 
substitutional atom in monovalent metals, or a monovalent substitutional 
ion in divalent metal, etc. . 

This behaviour seems to occur mainly when the components of the 
alloy have very different valencies. Thus Schubert gives numerous 
examples for alloys of polyvalent metals with transitional and rare earth 
elements (ZnNi, pure Ca, MgLa, etc.). In a similar way, Raynor (1949) 
assumes that transitional elements subtract conduction electrons from 
aluminium to build up negative ions with complete d shells. Truly conic 
configurations may be considered as an extreme case, where the 
concentration of negative metalloid ions is large enough to empty 
completely the conduction band of the metal. 

The difference in magnetic properties of Ni, Pd and Pt dissolved in 
noble metals may be explained in the same way. In Hume—Rothery 
alloys, it is well known that the transitional elements normally make no 
contribution of electrons to the conduction band, hence dissolve as neutral 
atoms. One may ask what sort of neutral atom. Ni for instance may 
dissolve as 3d!°, 3d94s or 3d84s?, etc., the 4s states remaining bound. 
As already explained for H~ in Cu, the possible 4s bound states are not 
to be subtracted from the conduction band. It is thus probable that, 
when dissolved in small concentration, these transitional elements have the 
same configurations as in a gas, paramagnetic 3d84s* and 5d%6s? for 
Ni and Pt respectively, diamagnetic 4d!° for Pd. Pd is actually 
diamagnetic in Cu or Au, Ni and Pt rather strongly paramagnetic (cf. 
for instance Vogt 1951). 


§9. APPROXIMATE EVALUATION OF ENERGIES OF SOLUTION IN 
Monovatent METALS 

We devote the last sections of this paper to applications to energies 
of solution (§ 9), paramagnetism of alloys (§ 10) and x-ray spectra (§ 11). 

The assumption of bound screening, necessary for As or Al in Cu, 
probably gives reasonable results for the energy of solution even when 
screening by Fermi electrons actually takes place, as is probable for 
Zn in Cu. In calculating these energies we shall take as an example 
aluminium in copper. 

The energy of solution of Al in Cu is the change of energy of an 
aluminium atom taken from the metal and put in substitution into a 
copper lattice. We obtain it by the following cycle : 


(1) First, take an aluminium atom from the bulk metal. If H, is the 
sublimation energy, the work required is H,. 

(2) Then take one electron from it, leaving an Al+ ion. The work 
done is the ionization energy H,. 
(3) Take a copper ion out of the matrix, putting an aluminium ion 
instead. The work done will be small; we denote it by 6. 


(4) Put the electron on the copper ion. If EZ,’ is the ionization energy 
of Cu, the work done is —E,’. 
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(5) Put the copper atom back in the lattice; the work done is —H’ 
where H,’ is the sublimation energy of Copper. 
The total work done (the energy of solution E,) is 


E,=H,+E,—H,’—E,' +8. 


Table 4 gives the computed values of H, for Cu and Ag alloys, 
compared with experiment. The measurements are accurate for Cu Zn. 
Ag Zn and AgCd* only (Birchenall 1949). The other experimental 
values were obtained by fitting some reasonable hyperbolic curve with 
results obtained at rather high concentrations (cf. Wagner 1940 
Smithells 1949). Values for H, are taken from Seitz (1940 a) antl 
Kubaschewski (1951). 

Discrepancies between computed and experimental values must be 
ascribed to the term 6 which we neglected altogether. 5 contains several 
contributions. First, a very small positive elastic energy, arising from the 
difference in atomic radii between solvant and solute elements. Then a 
term corresponding to Coulomb and exchange interaction between the 


Table 4. Energies of Solution in Cu and Ag (in ev) 


ike EG Tiley x, +H,comp. |+#,exp.| 
Cu Zn | 3-53-+0-09 | 7-68 | 1:365+0-02 | 9-36 —0-485-+0-11 | —0-35, | 
Cu Al 3-30 -+0-22 5-96 —1-:95+0-31 | (—1:5) | 
Cu Sn (3°38) 7:3 (—0-53) - (—0-5) } 
Cu As | © (1-31) 10 (+0-10) (—2-7) 
Ag Zn | 3-:00+0-06 | 7:54 | 1-365-0-02 | 9-36 | +0-165-40-08 | —0-18 | 
Ag Cd 1-16+0-02 8:96 —0-42+0-08 | —0-19 | 


conduction band and the bound states of the impurity. This is small for 
H in Cu (§7), but difficult to evaluate for substitutional impurities in 
noble metals (an evaluation will be made for light metals in § 11). 
Finally, the change in Van der Waals energy between ion cores is probably 
not negligible. 

We must thus not expect any great accuracy. Except for Cu As and- 
_ Ag Zn the agreement is actually reasonably good. It is interesting that 
bound screening gives satisfactory results even for atoms with an even 
number of electrons (Zn, Sn) for which two electrons are more probably in 
the conduction band. We may also point out that to put the three 
valency electrons of Al in the conduction band of Cu would require at 
least +15 ev, as is shown by a simple computation, in contrast with the 
small negative energy of solution of this element. It is therefore impossible 
for elements like Al, Sn, As... to have all their valency electrons in the 
conduction band of the matrix. 
Pe ee ee 

~* Cu Zn means Zn dissolved in Cu. 
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§10. MaGNETISM OF ALLOYS WITH NopieE MetTALs 


Measurements of the susceptibility of alloys would probably provide 
the most direct check of some of our results, and especially on the existence 
of incompletely filled bound states. We analyse in this section the few 
experimental data available. 

We shall not deal here with magnetism of alloys of noble metals with 
transitional elements, the interpretation of which seems to be of a certain 
complexity (Vogt 1951). We refer only to a qualitative explanation 
proposed in § 8 for para- and diamagnetism in alloys of Ni, Pd, Pt, with 
Cu and Au. ; 

From the few measurements available on alloys with polyvalent elements, 
the only reliable ones seem to be those on Cu Zn (Endo 1925) and Cu Al 
(Vogt and Harms 1943). 

For Cu Al, we assume that the contribution of the Fermi electrons of Al 
is the same as in pure Cu, and that the variation of atomic susceptibility is 
due to the replacement of Cut by Al* ions. The same assumption would 
lead, for Cu Au for instance, to a linear variation of y which is actually 
followed within 1 to 2 units « 108 (Shimizu 1932). 

Let then 

X=@Xai t+ (1—%)xXeu 

be the susceptibility for atomic percentage 100z% of aluminium. yey 
is the susceptibility of pure copper, and y,, is the contribution of the 
aluminium atoms, expressed in CGS units. For every aluminium atom, 
contributions to y,,; are from the Al*++* ion, the 3s? shell and the Fermi 
electron. The unknown susceptibility of the two 3s electrons of Al* was 
computed by Slater’s method (1930) using for the screening constant the 
value deduced from the ionization potential of Al. This gives 
X1=—21x10-%. Taking for Alt+++: y,—=—2-5x10-® (Pauling 1927); 
for Cut : y,=— 14-7 x 10-6 (Hartree and Hartree 1936) ; and for Cu metal : 
Xa= —5:3 x 10-6 (Vogt and Harms loc. cit.), we have for Al in Cu 


Xal= Xd Xs Xe ka ete te 
The experimental value is 
Y= aes 12-2 < 10-8, 


Owing to the approximations made, the agreement is satisfactory. 

For Cu Zn, it is interesting to note that the diamagnetism is increased 
when Zn is added. If the additional 4s electron of Zn was in a half occupied 
bound state, one would have expected a strong paramagnetic term to 
appear. We must therefore assume a screening by Fermi electrons 
(cf. § 6). Qualitatively, the effective mass in the Fermi gas increases 
slightly when Zn is added, so as to accommodate the extra electrons in the 
band whose width is unchanged. Hence the paramagnetic term decreases 
and the Landau diamagnetism increases (Mott and Jones 1936 d). 
Quantitatively however, this effective mass varies with the energy in the 
conduction band ; and it is simpler to assume as for the energies of solution 
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(§ 9) that the supplementary electronic charge has nearly the same distri- 
bution in space as a 4s bound electron in the field of Zn++. To compute 
the diamagnetism of this charge, we use the same method as above for the 
3s* electrons of Al*+, and take for the difference between the susceptibilities 
of Zn and Cu : +2-10-6 (Pauling loc. cit.). This gives 

Xen= — 24-1 X 10-8, 
The experimental value is 

X m=— 24-5 x 10-*. 
The agreement is better than the approximations warrant. 

Other measurements (Davis and Keeping 1929, Spencer and John 1927) 
show rather strong paramagnetic increases of the magnetic susceptibility 
for Cu Mg, Ag Pb and Cu Sb. They are probably due to paramagnetic 
properties of Mg*, Pb+ and Sb++ ions. The variations with concentration 
and temperature should however be more thoroughly investigated before 
any definite conclusion could be made. 

One could expect a priori a paramagnetic susceptibility inversely 
proportional to the temperature, as it is the case for transitional metals 
like Co, Fe in noble metals (Vogt 1951). Here however, the paramagnetic 
bound states have much larger radii and a direct coupling is probable 
already for small concentrations. We would also expect a rather strong 
indirect coupling by the electrons of the conduction band, which are 
strongly perturbed (§ 7). These two reasons seem to indicate that the 
paramagnetism should not vary much with temperature. 


§11. X-Rays anD OpticaL TRANSITIONS IN MonovaLent METALS AND 
THEIR ALLOYS 


We shall now use the results obtained in the preceding sections to 
study the absorption and emission of x-rays by monovalent metals. 

The absorption of an x-ray quantum by a metal is usually described as a 
transition of an inner electron to the continuum above the top of the 
conduction band. This process leaves however a metallic ion with an 
excess positive charge. Just as for the case of a Zn ion in copper, the 
conduction electrons of a monovalent metal cannot remain unaltered 
round this charge, but must provide some sort of screening, the most 
stable form being given by a bound electron below the conduction band. 
This electron, together with the ion which has lost an electron from an 
inner shell, forms an excited monovalent ion which replaces a normal 
monovalent ion without much effect on the conduction electrons. In 
calculating the energy of the absorption edge, we may therefore neglect as 
before (§ 9) the changes of energy in the conduction band, and consider 
only the excitation of the monovalent ion. 

The absorption process thus involves several simultaneous electronic 
transitions, and the selection rules only apply to the total resulting 
changes of symmetry. The positive hole in an X-ray level and the 
additional bound electron form together an exciton in Wannier’s sense 
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(Wannier Joc. cit.). We may thus state that in the main absorption or 
emission processes the exciton of lowest energy is respectively created or 
annihilated. We shall see that the presence of the Fermi gas broadens the 
emission into a band and makes possible transitions to an exciton of the 
same symmetry as the initial state. In metals, the successive excitons are 
widely separated in energy *those of energies higher than the fundamental 
one thus give rise to a fine structure on the short wave length side of the 
main emission or absorption edge. 

We develop the simplest case, Li K, then review various other cases 
where simple quantitative checks are possible: Nal, ;, and optical 
transitions in noble metals ; shifting of K absorption edge in Na and NaCl, 
K and KCl; and finally some K transitions in Cu alloys. 


11.1 K Absorption and Emission in Li 


We shall describe the transition process in a qualitative way, then check 
energy results with experiment for both absorption edge and fine structure. 

(a) We show first that the energy of the absorption edge corresponds to the 
excitation of a Lit ion from the state 1s? to 1s 2s. During the absorption 
(fig. 6) an electron is taken from the 1s? shell and sent into the continuum 
in a p state B. The hole A and the remaining Is electron D normally 
go out of the 1s? band into what we may call the 1s’ bound state* of a 
definite atom (A’, D’). This transition, which takes place when the Li 
atom is in a gas, reduces the absorption energy by DD’. According to our 
preceding results the now localized positive charge A’ requires screening, 
the most stable configuration being obtained with a bound electron C’ in a 
2s’ state. The field acting on this electron is mainly due to the Li** ion 
and the energy of the conduction band is not very much affected by the 
excitation.1s*>1s' 2s’ of the Lit ion. Small corrections, corresponding to 
the term 6 in § 9, will be considered later. ) 

In conclusion, we picture the K absorption by a triple electronic 
transition : a 1s electron A, going into B in the continuum, leaves a positive 
hole which goes with the remaining 1s electron D into a localized 1s state 
(A’D’). The charge of the positive hole is screened by an electron C’ 
produced from the Fermi band and put into the 2s state of the excited Li 
atom. ‘The positive hole A’ and the bound electron C’ constitute a Wannier 
exciton (Wannier loc. cit.). The total energy for the absorption is | 
approximately 

AB—CC’—DD’. 


The auxiliary transitions CC’, DD’ are not allowed in a one electron 
theory. They must be considered as induced by the perturbation of the 
photon at the same time as the main transition AB: AB? COSDE: 
represent a unique transition of the whole electronic system. Similar 
processes occur in every polyelectronic system, the simplest being the K 
absorption of helium where the remaining 1s? electron goes into a Is’ state. 


* Dashes for states of the excited atom. 
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(b) It may be shown that the absorption has a non-vanishing probability 
when the hole C is at the top of the Fermi distribution. There is therefore a 
discontinuous increase of absorption at the absorption edge, when the 
energy of the incident photon becomes sufficient to raise pleciron B above 
the Fermi distribution. This point is in accord with experiment. It may 
however, only be proved by a study of all the Fermi electrons, and m 
intend to develop it in a later paper. 


Fig. 6 
EeV 
O 


Or 


K absorption process in Li.* 


We shall only compute here the energy of this absorption edge, when B 
and C are both at the top of the Fermi band (fig. 6). The main term is the 
excitation energy H, of Lit 1s’>Is 2s. This is given by optical data 
(Bacher and Goudsmit 1932). 
We make two corrections. First, an exchange energy — H, between 
the 2s’ electron C’ and the 2s conduction band. We use here Fock’s 


* A’ should be a white circle (as A). 


178 _J. Friedel on the Distribution of 


formula (Fock loc. cit.) with Slater’s approximate wave functions for both 
2s and 2s’ (Slater 1930). A straight-forward computation gives for the 
exchange between (7, /,) and (7, /,), in atomic units 


= Hj Ad,* Gig” (4;+02)7 7. 
with a=2n,+1, y=2n,+1, a=2Z, nr, 


where Z, is the ‘ effective charge ’ acting on the electron, computed fron. 


Slater’s empirical rules. A isa numerical factor given in table 5 for various 
cases. As we used for 2s a bound wave function non-orthogonal to 2s’, 
values computed by (5) are only approximate. 


Table 5. Numerical Values of Factor A 


3s 4s 3s “48 4s 


1,056 3,328 13,200 


As a second correction, we compute the small negative energy —H#, due 
to the attraction of 2s’ by the potential of the lattice in the neighbouring 
atomic polyhedra. We assumed as in §7 that the form of 2s’ was not 
perturbed ; the potential of Lit was taken from Fock and Petrashen (1935), 
and the density of Fermi electrons assumed to be uniform. 

The presence of a polarized charge in the conduction band alters slightly 
the form of the 2s’ wave function, and hence changes its Coulomb inter- 
action with the Li*+ ion (cf. § 7). This is however smaller than for H in 
Cu where the bound electron is nearer to the conduction band. Assuming 
the same displaced charge as in Cu H, we found a maximum value of | ey. 
We neglected here this term altogether, and owing also to the approxima- 
tion in the exchange term, we do not claim an accuracy better than 1-2 ey. 


Table 6. K Absorption Edge and Fine Structure for Li (Energies in ev) 


Ky E. Ky, Total hy abs. hv em. 
28; 58-7 3:1 0-15 55:4, 5 54:5 
2 p’ 61 3.7 joann th 57-1, 585 
3p’ 69 3:00 2°15 63°85 63:5 


ee A ea ee ae el ee ee | hee I 


The total #;—H,—H,, thus computed (table 6) is in Saetartors 
agreement with the experimental value, taken from Skinner and Johnston 
(1951). Two points may be noted. First, the energy of the 2s’ electron 
is about —18-5 ev, definitely below the conduction band which extends 
from —5:5 to —l0ev. Also a screening by two bound electrons (2s?) 
and a positive hole in the conduction band is a priori possible (analogous 
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to configuration (c) §6). Using Slater’s approximation for the energy of 
bound states, one finds however that such a configuration is less stable 
by about 3 ev than the screening by one electron described above 
_ Aconfiguration 1s 2s? could thereforeplay a role only in the fine structure. 
(c) The fine structure to the short-wave side of the absorption edge Fe 
interpreted by Jones and Mott (1937) as a Kronig structure. Assuming 
that the average density of states is the same as in a free electron gas 
reflexions on (110) and (200) planes would correspond to two peaks a 
absorption at about 3-8 and 12-3ev from the absorption edge. Two 
strong lines are actually visible (fig. 7, taken from Skinner and Johnston 
1937) at 3-5 and 8-5 ev from the absorption edge. 


Fig. 7 


_ Absorption . 


260 240 220 200 BO ~——«IBO, 


Li K absorption (after Skinner and Johnston). 


The preceding analysis suggests however an alternative interpretation. 
The screening at the end of the absorption process could be obtained by 
a 2p’ (or 3p’... .) electron instead of 2s’. These configurations require 
a higher energy than the fundamental one with 2s’; but owing to their 
p symmetry, a direct Is+>2p’ transition is possible, without affecting 
the Fermi electrons. Hence we expect the presence of strong absorption 
lines on the short-wave lengths of the absorption edge. Table 6 gives their 
position computed by the same method as for the absorption edge. 
A stronger dependence of the fine structure near the absorption edge 
upon atomic than upon lattice structure has been reported by several 
authors (see for instance Coster and Kiestra 1948, Sanner 1943) and may 


probably be explained in a similar way. 
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(d) The emission process is just the reverse of absorption, starting 
with a bound electron 2s’ screening a hole in the Is shell (fig. 8). The 
edge of the emission band corresponds to the case when B is at the top of 
the Fermi band. Thus absorption and emission edges must correspond 
to the same energy, a well known fact in metals (cf. O'Bryan and Skinner 
1934). 


Fig. 8 


fev 


-0--4- bi 
BY N 


K emission in Li. 


Summarizing, both absorption and emission edges correspond to a 
transition which may be represented by an internal 1s*>1s 2s excitation of 
a Lit ion. The fine structure in the absorption may be attributed to 
further excitations 1s 2p, 1s 3p... . 

A similar quantitative study is possible for other monovalent metals 
as long as the excitation energy H, of the monovalent ion may be deduced 


from some other experiment. We shall now review a few examples 
where this is the case. 
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11.2. Na Ly, py, Emission and Absorption 


The description is similar to that for Li K, a positive hole in the 2p level 
being screened by a bound 3s electron. The small peak at the top of the 
emission band (Skinner 1940) may be due to the allowed direct transition 
2p§>2p°3s. Table 7 (first row) gives the computed and observed values 
for the emission edge. The two values for H, correspond to the two 
states 2p°3s of momentum J=1, which give rise to the doublets Lit um 
Ka;2, Ly, etc.... H, and EH, are computed as before (values “for 
Na* taken from Fock and Petrashen 1934). 

Skinner’s emission curve shows a small satellite at 33 ev followed by a 
more extended one between 42-5 and about 39 ev. The first satellite is 
readily explained by assuming an initial screening by a 3p instead of 
a 3s electron. Its intensity is small, for the 3p—2p transition is possible 
by the intermediary of the conduction band only. The 3d, 4s... . excited 
states seem responsible for the second (extended) satellite. The 3s, 3p 


Table 7. Ly, Absorption and Emission for Na (energies in ev) 


Ky E. Ey, Total hv abs. hy em. 


Bet 32°77 30-1, 30-5 30-25 
Gepaey| 0:37 
> GBPS 30°5 30-7 30°55 
3p’ 36-17 2°75 0:37 33-0; (32-4) 33:0 
Series | 
Limit 47-0; 0 5:16 41:9 42-5 


states are below the conduction band, as is easily deduced from data in 
table 7. The 3d, 4s... . states on the contrary are within the conduction 
band and should be replaced by equivalent screenings by Fermi electrons 
(cf. §9). Their series limits correspond to an electron displaced to 
infinity in the metal, and hence to the top of the Fermi band whose 
energy is computed with the conventions of §5 by 


— HB, = —E,—E,+0-4 Ey=—5-1, ev. 


The second and third rows in table 7 give the corresponding computed 
and. observed energies. 

We may note that the series limit (hv=42-5 ev) corresponds to the 
L emission of a singly ionized (1s2s?2p*3s) atom, and not to a doubly 
ionized one as is usually stated for satellites.. It corresponds in fact to 
the mechanism of emission as usually described. But the screening by 
a 3s electron reduces the energy from 42-5 to 30-2 ev. 

The absorption edge, as the values in table 7 due to Rule (1944) show, 
is in good agreement with the emission edge. To the short wave side of 
it, one sees in O’Bryan’s curves (1940) an absorption peak at ¢ 32:4 ev 
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which may correspond to the satellite at 33-0 ev. Here too, the Kronig 
structure would not give very different results ; one would expect maxima. 
at 2:55, 8-25 and 13-95 ev from the absorption edge. 

‘In conclusion, we suggest that the mechanism may be analogous to 
that for Lik. The satellites and secondary absorption peak may be 
‘explained by the presence of more excited Na* ions. 

11.3. Displacement of K Absorption Edge for Na and NaCl, K and KCl 

The position of the absorption edge of an element depends upon the 
nature of the solid in which the process takes place (ef. Pauling 1929 for 
instance). We give two examples where the major part of the displace- 
ment may be computed, and analyse the fine structure of the K edge of Na 
metal. ; 

(a) In NaCl, the Na K absorption edge (series limit) has an energy hv 
given by (Mott and Gurney 1948 b) 

hv—E—W,—P—y. 
Here £ is the ionization energy of a 1s electron in a Na* ion; —P the 
energy due to the polarization in the crystal round a point charge ; 
—y is the energy of the bottom of the band ; finally — Wy, is the lattice 
energy per ion pair, mainly due to the electrostatic interaction between. 
ions. 


! 


Table 8. Displacement of Absorption Edges of Na and K (ev) 


HK, W h(v—v')comp h(v—v')obs 


NaCl 14-9, 
KCl 11-8, 


In the metal on the other hand, we assume that the point charge is 
screened by a 3s electron; the sodium ion is thus in the state 
(1s 2s? 2p$ 3s). The frequency hv’ of the absorption edge is thus given by 

hy’ ~H—EH,—E,—E,, 
where —#, is the ionization potential of the screening 3s electron and 
—(H£,+,) the usual corrections for exchange and lattice interactions. ~ 


The difference between the two energies, giving the displacement of the 
absorption edge, is therefore 


hv—v') = (B, +B, +E,)—(Wy+P+y). 
We take as an approximation for the unknown £, the ionization potential 


of the 38 electron in Mg*, which is in a similar field, that is to say the 
second ionization potential of Mg. Thus —H,—14-9,ev. The quantities 


x, P and W;, are taken from Mott and Gurney (Joc. cit.). As we consider - 


Na and not Cl, we must however change the sign of the electrostatic term 
in W,. H,and HE, are computed as before. Table 8 gives the computed 
value for h(v—v’), in satisfactory agreement with experiment (Rule 1944) 
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In a similar computation for KCl, the potential of K+ used for H, is 
that of Hartree and Hartree (1938) and —#, is taken as the second 
ionization potential of Ca. The observed displacement is taken from 
Lindh (1930). 

(6) The fine structure on the short wave side of the absorption edge in 
Na metal may be studied as for Lik. The absorption edge itself is 
decomposed into three parts (K,, K,, Kg, in Rule loc. cit.), which may 
correspond to screening by 3s, 3p and 4s electrons. The absorption then 
increases slowly up to a maximum A which is approximately at the same 
wavelength as the series limit. Table 9 gives the computed and observed 
energies. In the first column, the unknown values of E, are replaced by 
the corresponding values for the normal Mg+ ion. As for Li K, the fine 
structure could also be explained by a Kronig structure. One would 
expect a maximum of absorption for the energies listed on the right in 
table 11; the agreement is of the same order. 


Table 9. Fine Structure of Na K Absorption Edge (ev) 


Deno- || Expected | Deno- 
AE, | —AE,|—AKy,,| Total Ahv exp mina- || Kronig | mina- 


tion ||Structure| tion 
BS: 0+0-1 KG 
3p’ 4-4, | —0-4, 0 3-9, | 3-0;+0°8 K, 2°54 (110) 
4s’ 8-6, O27 t= 1 Go| 2172, 5-6+0-1 Ke 8:24 (200) 
Series 
limit | 14:9, 2:1, | —4:2, |. 12-2, | 10-2,+1°4 A 13:94 Sian 


11.4. Optical Absorption in Noble Metals 

Copper, silver and gold show absorption in the blue or ultra-violet ; 
to this is due the colour of copper and gold. As was pointed out by Mott 
(Mott and Jones 1936 e), the main absorption seems due to internal 
photoelectric effect, with excitation of a d electron. We shall show that 
the same model as for x-rays can be used, and shall relate the absorption 
energy to the energy separation between the s and d parts of K and Ly, pry 
emission bands. : 

Optical absorption in noble metals is usually described as a transition of 
an electron from the d band to the top of the Fermi band, without screen- 
ing. A simple computation shows that in copper the top of the Fermi band 
is at about —8-2ev. The first ionization potential of the d shell, given by 
the second ionization potential of Cu, is —20-34 ev ; and the absorption 
edge should therefore be at about 12-2 ev, very different from the actual 
value of about 2:1 ev. ’ 

As for x-rays, we must therefore assume a screening electron, 4s for Cu. 
Thus the energy of the absorption edge should be near to the excitation 
energy H, of the inner d shell (3d" > 3d94s for Cu). Table 10 gives the 
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value of H, for Cu, Ag and Au (Cu and Ag: Bacher and Goudsmit loc. cit., 
Au: Platt and Sawyer 1941). The correction —(H,+4;) is probably 
small, as for Cu H (§ 7); but Slater’s approximation is very poor for such 
heavy atoms and a more exact computation would be necessary. The 
parallelism between E, and the experimental values taken from Mott and 
Jones (1936 e) is however satisfactory. 

As shown by Rule (loc. cit.), the Ly, 7, absorption in elements Na to Cl 
is equivalent to the difference between a Ka, and the normal K&£., 
emissions, and its energy is thus equal to the Kx, ,—Kf, separation. In 
a similar way, the energy of the optical transition in copper must be equal 
to the separation between the s band edge and the secondary d peak in 
the CuK and Ly yr, emissions. For Cu K, the distance between 
K8,(4sp > ls) and K£,(3d > Is) is about 3 ev (Bearden and Friedman 
1940). For Cu Ly yy, recent measurements by Cauchois (1951) give 
about 2-5 to 3ev. Both values agree with the peak 2-5 ev of the optical 
absorption. 

Table 10. Optical Absorption in Noble Metals (ev) 


Ey hv exp limit hv exp peak 
Cu 4s 2-90 2:1 2°5 
Ag 58 5-40 4-0 4-4 
Au 77.68 3°42 2:5 


11.5. ALK Emission Band in Cu Al and similar Emissions in other Alloys 

We shall now extend our results to alloys of noble metals. We start 
with Cu Al as an example, and verify some of our conclusions on other 
alloys. 

Recent measurements by Cauchois (1950, 1951), completing previous ones 
by Yoshida (1935) and Farineau (1939), present some characteristic features: 
the Kf emission band is divided between a rather sharp peak at the edge 
(width about 2 ev) followed by a second peak (width about 4 ev) at 4:5 to 
5 ev from the edge ; the position of the edge is practically independent of 
Al concentration. 

We picture as in § 9 an Al atom dissolved in Cu as an Al* ion plus one 
Fermi electron having the same energy limits and density of states as in 
pure copper, but with an s component repelled and p, d, . . . components 
attracted round the Al atoms. 

Before the K emission of an aluminium atom the positive hole in the 
Is shell of Al is therefore screened by an additional 3p electron, so that the 
configuration is (1s 2s? 2p* 3s? 3p), similar to an excited free ion. The 
electron which fills the positive hole during emission may come from the 
conduction band, the 3p state or the 3s? shell. The three cases are 
described in fig. 9 (a), (b), (c). In (a) the additional 3p electron must 
disappear during emission ; and in (c) the transition is allowed only by the 
intermediary of the conduction band. 


Electrons Round Impurities in Monovalent Metals 185 


; Transition (a) gives what we may call the weak ‘ normal ’ emission band, 
similar in form to KB, of Cu. Superimposed at the emission edge, there is 
a strong peak corresponding to the direct transition (b). Finally, displaced 
towards the long wavelengths by an amount equal to the excitation energy 
of 3s" —> 3s3p in Al‘, there is a second and stronger emission band, also 
similar in form to CuK. Table 11 (first row) gives the computed value for 
this displacement, in good agreement with Cauchois and Farineau’s 
results. Ifnow the Al concentration increases, the form and position of the 
conduction band may vary. But this has very little effect on the energy 
of excitation of the Al* ion, at least as long as the 3s? band is not too much 
widened. Thus the absolute position of the emission edge remains practi- 
cally constant. 


Fig. 9 


(2) (d) ate, 
Three types of ALK emission in Cu Al. 


Table 11. Separation between Secondary Peak and K Emission Edge (ev) 


AE, —AE, — AE, Total Ahv exp 


| Cu Al3s 3p —4:6, +0:0, 0 —4-6, —4, to —5 
Cu Zn 4p —5:9, +0-6, 0 —5°3, ce. —5 


In conclusion, we have satisfactorily explained the various characteristics 
of Al K emission in Cu Al. Some of the results are valid for other alloys. 
huss aS i 

(1) The position of the emission edge is practically independent of the 
(not too large) concentration of solute atom. 
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(2) K emission band of polyvalent elements dissolved in Cu presents one 
or more secondary peaks on the long wavelengths side (e.g. Bearden and 
Friedman loc. cit. for Cu Zn; Skinner and Johnston 1938 for Cu Be). 
These peaks may be interpreted as inner transitions producing excited 
solute atoms. Table 11 (second row) gives the computed displacement for 
Cu Zn corresponding to an excitation 4s > 4p. 4s and 4p are actually 
replaced by equivalent screenings by Fermi electrons (cf. § 9) 

The author is greatly indebted to Professor N. F. Mott, who suggested 
the present problem, for his constant help and encouragement during the 
investigation. He wishes also to thank Dr. D. Polder and Mr. J. 8. 
Plaskett for stimulating discussions. 


APPENDIX 


(1) DisPLACEMENT OF CHARGE IN A FERMI GAS 


Let V(r) be a spherically symmetrical potential such that r?V(r) tends 
to zero at the origin and at infinity. We denote by %, and ¢%,’ two solutions 
of the corresponding Schrédinger equation for two different values HZ and EH’ 


of the energy 
ay. [. _Ui+i 
Sa bake V) 7a v2 ‘| ,=9, 


0 _ Ul+)) 
at te | 2" s ‘| hy =0. 
We multiply the first equation by ¢', the second by ys, and subtract ; then 
Oy OP,’ 


—Efy, = Or = by ore 


An integration gives 


R 
2B B) | dé dr =|’ Set y, 
Or Jo 
Let now the energies tend towards the same value 
EK’ +> H=3k?. 


We may replace y by ,+(0p,/0k) dk, and (E’—E) by k dk. This 
substitution gives, after dividing both sides by dk 
R Op Os 07s R 
2k | 2 dp =| tL OEE ny OM 
Fated ysis Fst | 
We replace now %, by its asymptotic values for 0 and R. When 
r—>0, we have ae (A a constant). When r-> oo, we have 


%~sin [kr+-n(k)—}lr] (mq the phase shift). Hence, when R tends to 
infinity 


“R r d 
2h | pedr=k (e+ it) —4 sin 2(kR+7,—4Uz). 
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When potential V vanishes, we obtain a Bessel function 
-R ; 
2k | An J? (kr) dr ~ kR—4 sin 2(kR—Mn). 
0 


A subtraction leads therefore to the equality 


2 /R 
2 eld H 

x 2A akrJ?, (kr)| dr ae — cos n, sin (9 us 

mae [b> —$akrJ7, (kr)] dr iT 4. Fy 08m Sin (2kR+-7,—I7). 

The left member of this equation gives the amount of charge attracted 
by potential V in the /th spherical component of an electronic state with 
momentum # in a Fermi gas. The second term on the right is a sinusoidal 
function of zero average value. This gives no contribution if R tends to 
infinity ; and we have 

pe eee i dn 
ie peel krJ?., (kr drase 
7) [b?—gakrJ7.,(kr)] dr ede | 

This reasoning uses an approximate sinusoidal form for jj, at infinity. 
Using the WKB approximation, one can show that the terms neglected do 
not alter the final result. 


(2) WANNTER-SLATER’S THEOREM 


We make the following definitions: Vj, is the potential of a perfect 
lattice and V,, a perturbing potential, equal to zero at infinity. k being 
an arbitrary constant vector, 4,=¢, exp (tk . r) isa wave function solution 
of the perfect lattice for the energy H# and ¢, exp (tk.r) is the wave 
function of an electron with momentum k in potential V,. A simple 
mathematical derivation gives 


V*[ bibs exp (7k . r)] 
+ 2[E—V_—V,]¢i¢2 exp (tk . r)=2(Vd, . Vee) exp (Uk. r). 


bib exp (tk.r) is therefore an approximate solution in the perturbed 
lattice if the potential (Vd, . Vd2)(¢,2) 1 given by the corrective term is 
small. When nothing is known about V, the average value of the corrective 


energy is given by i a ¢,*V¢, dz which is easily proved to be equal to 


0 
‘| a [os 5 (hte | do—ik | or%, exp (—tk . r) dr. 


Both terms disappear in this expression if , is periodical inr. The best 
choice for k is therefore such that ¢, is a Floquet function. 

We may therefore state the following : for an electron of given energy in 
a perturbed lattice, we denote by ¢, the Floquet function in the perfect 
lattice corresponding to the same energy and by k the corresponding 
momentum: the wave function of the electron is approximately the 
product of ¢, and the wave function po=do(ik.r) of the electron of 
momentum k travelling in the perturbing potential. 


02 
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The corrective potential (V¢, . V¢2)(¢1¢2)~* is small in an s s band and 
for perturbing potentials V,, possibly Coulombian at the origin but varying 
slowly at large distances. Near the edge of a band, ¢, may be replaced 
by the wave function at the bottom of the band, with the use of an effective 
mass, as in the Wigner—Seitz method. 
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SUMMARY 


Cosmic ray stars in light and heavy elements have been investigated 
by the photographic plate method. The angular distribution of the 
particles from stars in heavy elements is shown to be isotropic, while 
that of those from light elements is not. Possible explanations of this 
are discussed. The numbers of stars with various numbers of 
alpha particles and protons are shown to be in accordance with a random 
distribution. This indicates that the emission of one charged particle 
has no influence on the nature of a subsequent charged particle. The 
emission of excited beryllium eight nuclei is also investigated and their 
frequency found to agree with the evaporation theory. 


INTRODUCTION 


Iv has been shown by a number of investigators (Harding, Lattimore and 
Perkins 1948, Le Couteur 1950, Page 1950, Hodgson 1950a, b) that 
the main features of the disintegration of nuclei after they have been 
highly excited by collisions with energetic particles of the cosmic radiation 
may be adequately accounted for by the evaporation theory. 

Phenomena have, however, been recorded (Harding, Lattimore and 
Perkins 1948, §5; Perkins 1950a) which do not agree with the 
evaporation theory. In this paper an account is given of a search for 
such anomalies in the angular distribution of the particles (Part I) and 
in the number of particles of various types (Part II). Part III is devoted 
to an investigation of the possibility of emission of unstable beryllium 
eight excited to various energy states. 


EXPERIMENTAL PROCEDURE 

The observations were made on disintegrations in Ilford C2 emulsion 
exposed horizontally to the cosmic radiation on the Jungfraujoch and 
processed in the usual manner. C2 emulsion was chosen rather than the 
more sensitive G5 so that the tracks of alpha particles could be readily 
distinguished from those of protons. The energetic particles emitted in 
the initial impact of the cosmic ray particle with the target nucleus are 
not therefore considered here, but only the lower energy particles emitted 
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during the subsequent de-excitation. Stars so near to the surface of the 
emulsion that one or more tracks could not be identified were omitted 
from the analysis. 

Throughout this work, all singly charged particles were included in the 
class of protons and all particles of charge greater than two in the class of 
alpha particles. 

Particles crossing the emulsion sometimes collide with hydrogen 
nuclei, producing a forked track which simulates a three track star. The 
relation between the spatial configuration of the tracks and the masses 
of the particles enabled them to be recognized and excluded from the 
analysis. 

The stars were classified according to their numbers of tracks, and, for 
those of less than seven tracks, according as they had a recoil or not. 
A recoil is a short, thick track due to the residual nucleus. Tracks less 
than four microns in length were counted as recoils. 

Stars of more than six tracks are nearly all due to the disintegration of 
the heavy elements silver and bromine. Of the stars with six tracks or 
less, Harding (1949) has shown that those with recoil fragments are 
mainly due to the heavy elements also. The majority of the small stars 
without recoils are formed in the light elements carbon, nitrogen and 
oxygen. The remainder are stars in the heavy elements for which the 
recoil is invisible owing to its short range or unfavourable orientation. 
These considerations enable a partial separation of the disintegrations in 
light and heavy elements to be made. 


Part I—THE ANGULAR DISTRIBUTION 


§1. INTRODUCTION 


The evaporation theory predicts that particles will be emitted at 
random in the centre-of-mass system of the excited nucleus. It is of 
interest to examine the extent to which this is fulfilled for nuclear 
disintegrations produced by cosmic rays. 

In an investigation of the angular distribution of emitted particles, it 
is important to distinguish between the departures from isotropy caused 
by a movement of the emitting nucleus with respect to the apparatus, 
and those caused by an anisotropic mode of disintegration in the centre 
of mass system. These two sources of anisotropy will be referred to as 
extrinsic and intrinsic respectively. 

The extrinsic angular distribution of the particles emitted from nuclear 
disintegrations caused by cosmic rays has been investigated using 
photographic emulsions by Brown, Camerini, Fowler, Heitler, King and 
Powell (1949), by Bernardini, Cortini, and Manfredini (1950) and by 
Manfredini (1951). They measured the angular distribution of the star 
particles projected on a vertical plane and found that the fast particles 
emitted in the initial impact were strongly collimated downwards. Those 
lower energy particles emitted during the subsequent de-excitation process 


192 P. E. Hodgson on Some Features of 


were slightly collimated downwards in the stars of more than six tracks, 
and more strongly collimated downwards for the smaller stars. 

This downward collimation of the evaporated particles can be 
attributed to the motion of the excited nucleus. The increase in 
collimation for the small stars is in accordance with expectation since, 
for the same momentum transfer, a light nucleus will be given a greater 
downward velocity than a heavy one. 

Perkins (1950 a) has used the numbers of alpha particles and protons 
emitted upwards and downwards from large stars to calculate the down- 
ward velocity of the excited nucleus. For stars of up to eleven tracks he 
found that the downward velocities calculated separately from the 
alpha particles and the protons are the same. This is what would be 
expected if all the particles were emitted isotropically in the centre of 
mass system. For stars of more than eleven tracks, however, there was a 
marked discrepancy, which he attributed to an anomaly in the angular 
distribution of the alpha particles due to an abnormal mode of 
disintegration which has not yet been explained. 

In order to examine the mechanism of disintegration, the intrinsic 
angular distribution of the star particles has been investigated by several 
authors using different methods. 

Perkins (1948) and Harding, Lattimore and Perkins (1949) have 
associated a unit vector with each track and evaluated the resultant for 
each star. The observed distribution of resultant vectors was then 
compared with the distribution expected on Chandrasekhar’s random 
walk theory (1943). A significant difference between the two distributions 


indicates anisotropy. This was found only for alpha particles from very _ 


high energy disintegrations in which fifteen or more charged particles were 
emitted. However, this method is laborious since it involves measuring 
the orientation of every track in the star and calculating the resultant 
of unit vectors laid along them. 

Another method, developed by Lovera (1949), consists in measuring 
the distribution of all the angles, projected on the plane of the emulsion, 
between adjacent pairs of star particles. If the particles are randomly 
emitted from a massive nucleus, the theoretical distribution of these 
angles, for an n-track star, is given by 

W—T n—2 
P (¢)dd= oe (1- £) dd 0<¢ <2n, 


“7 


A variant of this method, developed independently at Imperial College, 
consists in measuring the least projected angle ¢, on the plane of the 
emulsion, between any two tracks of the star. The corresponding 
theoretical distribution for random emission is given by 
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In the following sections the results of measurements of this quantity 
are described and discussed. 

This method, like the others described here, does not provide any 
a priort way of distinguishing between intrinsic and extrinsic anisotropy. 
This must be done from a consideration of the conditions of exposure 
or of the type of disintegration studied. (See §3.) Thus, for example, 
extrinsic anisotropy would have a much greater effect on the distribution 
of ¢ for stars on plates exposed vertically than for those on plates exposed 
horizontally, whereas little difference between the two distributions would 
be expected from intrinsic anisotropy. Comparison of groups of stars on 
plates exposed in both ways can therefore be used to distinguish between 

the two kinds of anisotropy. 


Fig. 1 
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Distribution of least projected angle for three track stars with recoil fragments. 


§ 2. EXPERIMENTAL RESULTS 


The least projected angle was first measured, considering all the tracks, 
for each star in a randomly selected group of a thousand stars. The stars 
were then grouped according to their numbers of protons, and the least 
projected angle measured for each, only the proton tracks being 
considered. This was repeated considering only the alpha particle tracks. 
Thus a particular star can be included in up to three groups. The stars 
were classified according to the number of tracks N considered in the 
measurement. The distributions of least projected angles for the various 
groups of stars were compared with those expected on random emission. 
Typical results, for three track stars (with and without recoils), are shown 


in figs. 1 and 2. 
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The agreement between the two distributions was examined in each 
case by the x? test, and the probabilities that the observed distributions 
are chance fluctuations on the theoretical are given in table 1. The results 
in the columns headed ‘NV proton tracks’ refer to measurements in 
which only the proton tracks were considered, and similarly for the 
columns headed ‘NV alpha tracks’. Thus the three measurements on 
a five-track star of three protons and two alpha particles appear in three 
columns of the table. The results for two alpha tracks only refer to 
measurements on stars with just two alpha particles, irrespective of their 
numbers of protons, and similarly for the other groups. Results for the 
large stars are omitted because of the small number available. 


Fig. 2 
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Distribution of least projected angle for three track stars without recoil 
fragments. 


Table 1. Probabilities, according to the ,?-Test, that the Observed 
Distributions of € are Random Fluctuations from the Theoretical. 
Distributions 


pub Stars without recoils 7 Stars with recoils m~-meson 
) 
stars 

tracks | All N |Nproton|N alpha} ANN |N proton} NV alpha (all 

N tracks | tracks | tracks | tracks | tracks | tracks tracks) 
“3 0:8 0-1 
7 0-9 0-008 
5 0-2 0:07 


0-06 0-02 
0-003 0-2 0-01 


0-1 0-2 
0-2 Sees *¢ peste. 


0 

: (): 

0-2 0.2 r 0: 
Q- 
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: These results show that the agreement between theory and experiment 
= aaa Lage stars with recoils. Therefore the particles from 
e disintegrations of hea i are 1 i i 
gee as nee sea are randomly emitted, in agreement 
The probabilities for the small stars without recoil fragments, which 
are mainly the disintegrations of light elements, are in many oe less 
than 0-05. This is commonly considered to indicate a significant difference 
between the two distributions. Since emission from stars in heavy 
elements has been shown to be random, this result shows that emission 
from the stars in light elements alone is even more anisotropic than 
from the recoil-less stars considered as a whole. 


Fig. 3 
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Distribution of least projected angle for three track stars initiated by 
am—-mesons. 


This conclusion is supported by the results of similar measurements on 
two-, three-, and four-track stars initiated by 7~-mesons, which are nearly 
all formed in light elements. (Perkins 1949, Menon, Muirhead and Rochat 
1950.) These results are included in table 1, and the distribution of least 
projected angle for three-track 7~-meson stars given in fig. 3. These were 
shown by the x? test to be consistent with the corresponding distributions 
for small stars without recoils. 

The results show no significant difference in angular distribution 
between alpha particles and the protons. 

Possible explanations of this anisotropy of particles emitted from stars 
formed in light elements are discussed below. 
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§ 3.. Discussion 


The anisotropy cannot be ascribed to motion of the excited nucleus, 
since this would increase the number of stars with small least projected 
angles, contrary to what is observed. Also, since the plates used in the 
present experiment were exposed horizontally, motion of the excited 
nucleus could only have a small effect on the least projected angle 
distribution. The effect would be zero if all the star producing radiation 
was incident vertically. Therefore the observed anisotropy must be 
intrinsic, and due to the mode of disintegration of the light elements. 

In order to see whether stars produced in one specific nuclear reaction 
would show an anomalous distribution the method was applied to the 
C12(, 3a) stars observed by Goward and Wilkins (1951), who kindly 
sent their measurements. The results are plotted in fig. 4, and indicate 


. 


Fig. 4 
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Distribution of least projected angle for three track “C (y,3a) stars observed 
by F. K. Goward and J. J: Wilkins. 


a marked anisotropy. This supports the view that the emission of 
particles from a light nucleus is anisotropic. The results for small cosmic 
ray stars of course represent the average over a large number of different 
reactions. It is likely that many of the anisotropic features cancel out 
when such a group of stars is considered, but, nevertheless, some evidence 
of anisotropy remains. 

A possible explanation of this is that a particle is often emitted when 
the nucleus is recoiling from the emission of a previous one. This causes 
the angle between the directions of motion of the two particles to be larger 
in the laboratory system than in the centre-of-mass system, and this 
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increases the least projected angle. The effect is quite small owing to the 
disparity between the velocities of recoil and emission, and also because 
a neutron must often intervene between successive charged particles, 
destroying any angular correlation between them. 

The observations may also be accounted for by assuming that a light 
nucleus, when struck by an energetic particle, is shattered in one act, 
instead of evaporating via a number of intermediate nuclei as it doés 
when struck by a low energy particle. The mutual obstruction of the 
particles would cause them to be more symmetrically emitted than in the 
case of random evaporation. This also would result in an excess of stars 
with large least projected angles. 

It is not possible to choose between these alternatives without a 
detailed calculation of the effect of the nuclear recoil on the distribution 
of least projected angles. This calculation is now in progress and will be 
reported on completion. 


Part [I—TuHe FREQUENCIES OF STARS WITH VARIOUS 
NUMBERS OF ALPHA PARTICLES AND PROTONS 


§1. IyrRopUCTION 


Several authors (Perkins 1950 b, Crussard 1950 a,b, Hodgson 1950) 
have shown that pairs or triplets of particles due to the emission and 
subsequent disintegration of unstable fragments are sometimes observed 
in nuclear disintegrations caused by cosmic rays. In general these 
fragments are only seen if their energy of break-up is small, since then 
the disintegration products continue nearly in the direction of the original 
fragment, giving a narrow cone of particles which is readily recognizable. 
If, on the other hand, the energy of break-up is large, the disintegration 
products are widely separated and are lost in the background of 
evaporated particles. It is therefore of interest to examine the distribu- 
tion of alpha particles and protons among the stars for any anomalies 
which might be ascribed to unstable fragments. . 

If the emission of one particle does not affect the probability of emission 
of another similar particle, then the frequencies of emission will follow 
Bernouilli’s Law. Thus, if p, is the probability of emission of an alpha 
particle from a disintegration in which charged particles are emitted, 
then the proportion of a large group of such stars with just m alpha 
particles is 

n! m —m 3 ee enj-1 
Taare oe Cp)" P= l1—(a[p)*T*, 
where («/p) is the ratio between the numbers of doubly and singly charged 
particles emitted from the disintegrations considered. 


§ 2. COMPARISON OF RESULTS WITH THEORY 


One thousand stars were examined and the numbers of alpha particles 
and protons in each recorded. 
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The distribution given by the above formula were calculated for each 
group of stars and compared with the experimental by the x” test. 
Typical results, for three-track stars with and without recoils, are shown 
in fig. 5. The probability in each case that the observed distribution is a 
random fluctuation on the theoretical is given in table 2. 

Since stars without recoils are a mixture of stars in light and stars in 
heavy elements, the relevant calculations should take this into account, 
and the sum of two terms of the above form with the «/p ratios of stars 


in the light and heavy elements evaluated. Using approximate values 
Fig. 5 
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Table 2. Probabilities, According to the y? Test, that the Observed 
Distributions of Particles among the Stars are Random Fluctuations 
from the Theoretical Distributions 


rors Stars Stars 
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3 0-9 0-7 
4 0:2 0-1 
5 0-1 0-8 
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of the z/p ratios and of the relative proportion of the stars in light and 
heavy elements, this calculation was made, and the results are shown in 
fig. 5 by the dotted lines. It can be seen that this refinement makes 
little difference to the final result. 

Table 2 shows that in no case is there a significant difference between 
the observed and calculated distributions. These results therefore give 
no support to the theory that.some alpha particles come from the 
disintegration of excited fragments and consequently tend to come out 
in pairs or triplets, but, of course they do not exclude the possibility 
of the occurrence of such fragments in a small number of cases. 


Part IIJ—Tue Emission or Beryutivum EIGHT FROM 
NvucLEAR DISINTEGRATIONS 


Perkins (1950 b) and Crussard (1950 a,b) have observed close pairs 
of alpha particles emitted from nuclear disintegrations caused by cosmic 
rays. Their frequency of occurrence is greater than would be expected 
from random emission, and they are ascribed to the emission of $Be 
nuclei which almost immediately decay into two alpha particles. 

The energy of disintegration of the beryllium nucleus can be calculated 
from the energies of the two alpha particles and the angle between their 
directions of emission. The values obtained are about 90 kev in 
agreement with those obtained by other investigators (Hemmedinger 
1948, 1949; Tollestrup, Fowler and Lauritsen 1949) using artificially 
produced beryllium nuclei. 

In order to extend this work to other possible excited states of 
beryllium eight, measurements were made on all the pairs of alpha 
particles, both of which ended in the emulsion, from five hundred 
disintegrations. The energy of excitation, assuming the alpha particles 
to be from the disintegration of beryllium eight, was calculated for 
each pair. The distribution of energies of excitation was compared with 
that expected if all the particles were randomly evaporated. If any 
of the alpha-particle pairs are due to the emission and subsequent 
disintegration of beryllium eight, then the corresponding energy of 
excitation should appear as a peak on this distribution. It was found 
that the two distributions were consistent, except for an excess of eight 
pairs with small excitation energies due to the ground state of beryllium 
eight. Correcting for those pairs passing out of the emulsion, this gives 
a mean frequency of emission of 0-024+0-008, in good agreement with 
the approximate value 0-02 calculated by Le Couteur (1950 b) on the 
evaporation theory. The mean value of their disintegration energies Is 
82 kev, in good agreement with the results of other authors. The results 
give no evidence for the emission of beryllium nuclei excited to higher states, 
but they do not exclude their occurrence in a small number of cases. ; 

Recently, Millar and Cameron (1951) have reported similar events in 
which the track of the beryllium nucleus is visible. They find a value 
of (5-3-+1-1)« 10-4 sec for its lifetime. Wilkins and Goward (1951), 
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however, have examined their stars and find that it is only possible to 
‘put an upper limit to the lifetime of about 5x 10714 sec. The more 
energetic particles emitted from the disintegrations considered in this 
paper were therefore examined, and the distance from the star centre to 
where the two alpha particle tracks were resolved was measured and 
divided by the velocity of the particles as determined from their ranges. 
These results, in agreement with those of Wilkins and Goward, showed 
that it is only possible to put an upper limit of about 2 10~™ see on 
the lifetime of beryllium eight. ; 
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ABSTRACT 


The equilibrium position of an edge dislocation near certain types of 
lattice irregularities is investigated quantitatively, using the Peierls 
model of a dislocation. The external shear stress required to maintain 
equilibrium is calculated and the result extended to an array of like 
dislocations piled up against a barrier. The integral equation for the 
atomic displacements is linearized by a perturbation technique, the 
resulting equation being solved exactly. 


§1. INTRODUCTION 


It is well known that dislocations in a crystal lattice interact with each 
other, and with other structural defects and discontinuities. Of these, 
commonly occurring examples are solute atoms, vacant lattice sites, 
grain and crystal boundaries, micro-cracks and changes in elastic constants 
or lattice parameter. Interaction effects can be of great importance 
for the understanding of the mechanical properties of crystals, as is 
shown, for instance, by the work of Taylor (1934), Cottrell (1949), Mott 
and Nabarro (1947) and Seitz (1950). The main purpose of the present 
paper is to attempt to give a quantitative analysis of the effect on a 
dislocation of certain types of barriers and discontinuities which are 
thought to exist in actual crystals. More precisely, we set up a general 
method for calculating the external shear stress 7 required to hold 
a dislocation in equilibrium at a given distance s from a specified lattice 
inhomogeneity. 

At the very outset of the work it is necessary to exclude from 
consideration all barriers which involve a discontinuity in elastic 
constants, for these have so far proved to be mathematically intractable. 
Thus, for instance, it appears to be quite impossible to calculate 
+(s), even approximately, for the approach of a dislocation towards 
an ideal surface of separation in a bicrystal as illustrated in fig. 1. 
Fortunately, this type of barrier is not of very great interest from the 
point of view of mechanical properties. Most of the inhomogeneities of 


* Communicated by Professor H. Jones. te : 
+ Koehler (1941) has analysed the equilibrium of a dislocation in a finite 
single crystal, and since this is an extreme case of a bicrystal (u,>0, 4,=0) his 


method may possibly be generalized. 
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physical interest can be effectively represented by a model in which there 
is a change in the interatomic law of force across the slip plane, but no 
change of elastic constants elsewhere in the material (which is assumed. to 
be elastically isotropic). This means that whilst there occurs an alteration 
in the nature of the ‘glue’ which binds together the two dislocated. 
half-crystals, there is no effective alteration elsewhere. A direct example 


Fig. 1 Fig. 2 


< —Tt 
Gs 
The rigidity modulus jumps from ,, to A micro-crack in the slip plane; the 
ty on crossing the barrier. external shear stress 7 is negative 


in this case, for the crack attracts 
the dislocation. 


A change in lattice parameter blocking A dislocation near a grain boundary. — 
the path of a dislocation. The broken line represents the — 
hypothetical extension of the slip 
plane. 


of such a situation is provided by a Zener micro-crack (1946), this being 
defined as a region of the slip plane across which there exists no law of 
force between the atoms (fig. 2). Another example is provided by a 
material in which there is a change of lattice parameter over a certain 
region, as illustrated in fig. 3. A third and important type of 
inhomogeneity is provided by a grain or twin boundary. <A dislocation 


Interaction of a Dislocation with a Lattice I nhomogeneity 203 


is obviously repelled from such a barrier, and this can be regarded as 
being effectively due to a ‘hardening of the ‘ glue’ across a hypothetical 
extension of the slip plane (fig. 4). The barrier is defined mathematically 
as the point on the slip direction at which the law of force across the slip 
plane alters. 


§2. THE Prterts MopeL or a DisLocatTion 


In order to calculate the interaction between dislocations and barriers, 
we adopt the Peierls model of a dislocation (Peierls 1940, Nabarro 1947), 
in which the atomic structure of the lattice is taken into account along the 


Fig. 5 


Cross-section of a positive edge dislocation. 


slip plane. Applied to a dislocation in a homogeneous single crystal, 
this is successful in leading to stresses and strains which do not diverge 
in the neighbourhood of the origin, and which at the same time 
approximate to those of the Burgers classical dislocation elsewhere. 
The extension of this analysis to the case where the atomic structure 
along the slip plane is no longer constant but varies in a specified manner, 
constitutes the essential mathematical problem of the present paper. 
A diagram of the cross-section of an edge dislocation is shown in 
fig. 5, the elastic problem being a two dimensional one of plane strain. 


P2 
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The dislocation is constructed by cutting the material along x’Oz, inserting 
an extra row of atoms along Oy and applying a system of equal and 
opposite shear stresses to the surfaces A’A and B’B until the corresponding 
atoms in each move as far as possible into alignment. These surface 
stresses are the horizontal component of the interaction forces between 
corresponding atoms in the two surfaces, the normal component being 
neglected. 

The displacement w(x) of an atom in A’A will be approximately equal 
and opposite to that of the corresponding atom in. B’B, so that their 
relative displacement 7 is given by 


(1) 


the x-axis being the direction of positive u(x). If for convenience we adopt 
the non-dimensional measure of displacement %, where 


p2u/b, NO ee 
then (1) becomes 
7 (ac) O==ah(a@) he re coe Se, icy 


For a dislocation in equilibrium in a homogeneous single crystal there 
must be perfect alignment of corresponding atoms at infinity, so that 
r(-+- cc)=0, leading to 

ys( 00) = —y(— 00) =—}. oe eS res 
Ifthe crystal is not homogeneous then an external shear stress 7 is required 
to maintain equilibrium. The shear stress in the slip plane accordingly 


reduces to 7 at +00, so that if this is small, we have by Hooke’s law 
r(-+ 00)/b=7/u leading from (3) to 


yb( 00)= —4-+ (7]u), | 


— 
oO 
— 


where p is the shear modulus. 
Assuming the two half-planes are isotropic elastic continua,* the 
relation between the shear stress p(w) and du/dx along A’A is given by 


eee © -du(t) dt 
pel= | eit ein 2) le hee (6) 


where v is Poisson’s ratio (Nabarro 1947) ; the Cauchy principal value is 
taken at t=a. If distances 2, y, s are measured in units of f, where 
f=6/2(1—v) and the notation (2) is adopted, then (6) becomes 


peep SO Tsp cok Ms iphone aT 


ee eee 
*The present theory may be easily modified for an anisotropic crystal, 
following the method of Eshelby (1949). : 
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it being noted that p(+co)=r. An alternative integral relation, which is 
quite equivalent to (7) is ; 

des (x [ime . 

oe | GOR... 
this being determined either as the Hilbert transform (Titchmarsh 1948) 
of (7), or again directly from elasticity theory (Leibfried and Lucke 1949). 
Since the shear stress distribution is provided by the interatomic forces 
across the slip plane, it follows that p is a function of 7 and therefore 
of , ie. p(x)=f[p(x)], where f is a function either known or assumed. 
On now substituting for p(#) in either (7) or (8), we obtain the integral 
equation required. 

The function p(r) has not yet been calculated, but over any homogeneous 
region of the crystal it must be periodic with period b and alternately 
positive and negative in each half-period. Furthermore, it must 
everywhere have the limiting form p=yr/b for small relative displace- 
ments, in accordance’ with Hooke’s law. A simple function of this 
type, viz., 

p=(u/27) sin 27r(r/b)=—(p/27) sin 2b 6 ww. (9) 
was put forward by Peierls as the most convenient assumption as to 
the law of force in a homogeneous single crystal. On setting this into (7), 
we obtain the integral equation 

eal 5sin mp (a) —7= = | oa ate es (10) 
which, although non-linear and falling into no standard form, has the 
exact solution 

(x)= — (1/7) tan~1 2, ee ees «(11) 
for s=0, as may be verified by substitution. Since this also satisfies 
the boundary condition (4), it represents a single dislocation in 
equilibrium in a crystal free from external stress. _ Another exact solution 
for non-zero + is also known (Nabarro 1947), this representing a pair 
of dislocations of opposite sign held apart against their mutual attraction 
by the external shear stress. 

If the crystal is not homogeneous, this means, as pointed out in the 
introduction, that the p(r) function changes over the appropriate region 
of the slip plane. Thus, for instance, the presence of a Zener micro-crack 
of length / at the point s along the slip direction is represented by 


p=(~/27) sin 2x (7/b) ;  — OSS, 
=0; s<a<s-+l, 
=(p/27) sin 2n(r/b);  s+l<a<o. 
As another example, a change in lattice parameter from 0 to 6’ at s means 


that 
p=(u/27)[sin 2n(7/b);  — oes, 


= (1/27)(b'/b) sin 2n(7/b’);  s<w< OO. 
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As a third example, the effect of a change in lattice orientation at s may 
be described by — 
p=(p/27) sin 2n(r/b); — ©Nx<S, 
=prlb; 8s<“%<00. 
This is plausible since the absence of periodicity on the extension of the 
glide plane is represented by a restoring force which increases steadily 
with relative displacement (Hooke’s law instead of a periodic one). 


Fig. 6 


P (ee) 


The general form of the function p(r) when a dislocation lies near to a Hooke’s 
law barrier. 


The two laws of force in this case are illustrated by the discontinuous 
P(r) function in fig. 6, r= +(b/2) being the relative displacement at the 
origin (i.e. the centre of the dislocation), and r=r(s) being the relative 
displacement at the barrier, where the law of force changes. A further 
example is provided by a solute atom lying in the slip plane. Owing to 
| the difference in size and electronic structure between solute and solvent, 
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the solute will interact with the dislocation. The former effect has 
already been investigated (Cottrell 1950), whilst the latter can be dealt 
with to a rough approximation by the methods of this paper. Thus 
a solute atom* at x=s can be regarded as a region along the slip plane 
such that 
P(r)=(u/2m) sin 2n(r/b), —co<e<s—(6/2), 
=(1'/2m) sin 2n(r/b),  s—(b|2)<e<s-+(b/2), 
=(u/2n) sin 2n(r/b), s+ (0/2)<ax<o, 


the parameter yu’ being suitably chosen. 


§3. THE EXTERNAL SHEAR STRESS 
The integral equation for ;; may be written from (8) as 
ap ple) 

—a = [ oe Es mate | alles ses eee) 
where p is a given function of % for each separate region of the slip plane : 
more concisely, 

cfs Ly? pools 4, 

= Sah el ae ets 13) 
where the positions of the discontinuities are now parameters in the 
functional relation between p and #%. It is convenient to take the Hilbert 
transform of (13), i.e 


piiley—r= Ef SE. 


which equation, together with the boundary conditions (5), will be adopted 
for the purposes of calculation. 

A simple expression for the external shear stress 7 may be derived as 
follows. Multiply through Oo) ) by db/dx and integrate over all values of 


2x to give 
= ap 
| @-ngeail 


pes! dis dt 

where Ji | Fede [ eee ar 
pe! = ae dis dx 
=-| dt wo dx x—t 


assuming that it is valid to interchange the order of integration. Since 
I,= —I,, it follows that I, is zero, giving 


[° pga] un] ye ee eater LD) 


* More correctly we should talk of a row of solute atoms in the slip plane and 
parallel to the dislocation line. a, 
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On now making use of the boundary conditions (5) on ws, we have 


—r=r[yo |” =|" pga ees fe 76755 
=| pty 
=|) ptm | id 


~1/2 . rie 
=) pdtr/d) +| pd(r/b) 
t/u 1/2 
from (3), i.e. 
1/2 
r= | p(r/b)d(r/b). iS a ae 
—1/2 


From the result* (17) it will be seen that the magnitude of the external 
shear stress is that of the shaded area shown in fig. 6. Although in this 
diagram the position 7(s) of the discontinuity is not known without 
calculation, it can usually be estimated in cases where the interaction is 
not too large. For example, if a Peierls dislocation lies outside a Hooke’s 
law region, then it is a sufficiently good approximation to write 


r(s)= (6/er) tan=* (1s) 2 ee 


this being the relative displacement at x=s in the absence of the barrier. 
The resulting curve of 7 against s, derived from the equation 


1/2 Y , *1/atan-! 1/s , ee taf , 
vi JP (G)4Gs) =e], G~ e827) 465) 


SES (Prete: \ eee ee >0 19 
mei bGery) wae fcc 
is illustrated in fig. 7. For s>1, 7(s)=p/67s4, i.e. 
4 b4 
nite = (20) 


67234 967?(1—v)4s*’ 


on reverting to ordinary units of length. This gives s=5b for 
7~10-*u(v=1/3), which means that a typical yield stress pushes the 
dislocation to within five interatomic spacings from the barrier. It is 


*F. R. N. Nabarro has pointed out that this result may also be derived by 
the following virtual work argument. If the whole pattern of the displacement 
curve is moved bodily forward by da, i.e. the dislocation is moved forward, 
then the external shear stress does work tbda, which must be equal to the 
virtual change in the internal energy of the system. The strain energy of the 
two elastic continua (fig. 5) is unchanged, but the relative displacement of 
atoms in A’A and B’B changes by 


dr=1(% —Sa) —1(v) = —8a(dr/da) = —bda(dib/dx), 
so that the stresses p(x) do work 


—b3a i *  Pidiidr) de. 


On equating this to rbda we are led to (16), from which (17) is readily obtained. 
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interesting to note that when the centre of the dislocation coincides with 
the barrier, the relative displacement there is always 6/2, in agreement 
with the result obtained by setting s=0 in (18), which estimate is in this 
instance therefore strictly correct. An immediate consequence is that 
7(0) can always be calculated exactly, e.g. for a Hooke’s law barrier 
7(0)=(/27)[(7?/4)—1]=7-4 x 10-2u on setting s=0 in (19). 

It is obvious on physical grounds that, as illustrated in fig. 6, 7 will 
increase* as s decreases, so that 7(s)<7(0). Whether the dislocation 
can proceed further and penetrate the barrier depends on the nature of 


Fig. 7 


cS) 
SLL 


0-005 


ia BARRIER 


Alen 


= 


SMALL CRACK 


LARGE CRACK 
—0-005 


Graph of the external shear stress 7 required to hold a single dislocation at a 
distance s from a Hooke’s law barrier, a large crack, and a small crack 
(l=5B=~4b). The separate plots © are based on the more accurate values. 
of r(s), derived from the perturbation method. 


the law of force in the region concerned. If it is not periodic with 
period b, then + increases so rapidly that the dislocation is effectively 
blocked, as with the Hookean barrier above. If, on the other hand, the 
periodicity of the law of force is maintained, then 7 attains its maximum 


* In some cases, e.g. a micro-crack, the dislocation will be attracted towards 
the barrier, rather than repelled ;_ and if it can proceed through then 7(0) will 


be the minimum. 
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value at s=0, ie. Tmax=7(0), as can be seen by reference to figs. 8 (a) 
and 8 (b). In fig. 8 (a) the dislocation is at the barrier, and the shaded 
area represents 7(0). In fig. 8(b) the dislocation is at a distances beyond the 
barrier and the discontinuity occurs as shown, i.e. 7(s) is the sum of the 
shaded areas (of opposite sign) and therefore less than 7((). 


$4. ARRAYS OF DISLOCATIONS 


We now consider the more general problem of an array of n like 
dislocations piled up against a barrier by the stress 7 (fig. 9). The integral 
eqn. (14) still applies here, but has not yet been solved. However, the 


Fig. 8 


PO) 


; (a) (d) 
The general form of the (7) function (full line) when the centre of the dislocation 


lies (a) at the barrier, and (b) inside the barrier. 


applied shear stress may be estimated in much the same way as before, 


the only difference arising from the altered boundary conditions which 
are now 


Mp comes tniayeert: : 
cola eee ou 


It follows from this that [(x)]_2=—n so that the result (15) is 


modified to 
Lael (ae "\ alt ib 
i— n)}_ A (5) (5) = a are Vie) xe, Mas ch eins (22) 
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the integral J having the same geometrical interpretation in the (p, 1) 
diagram as before. It is still necessary to know the relative displacement 
at x=s, but in most cases an estimate of the type (18) may again be used 
provided that the dislocations behind the first are relatively widely ppaecds 
This is usually the case if 7 is small, because pairs of like dislocations 
at a distance R apart repel one another as (1/R) (Cottrell 1949) whilst 
a barrier usually repels dislocations as (1/s*) or (1/s*) ; accordingly s will 
be small compared with the spacings of the other dislocations. 

It will be noted that the integral J in (22) is identical with that occurring 
in (17),* which gives the external shear stress required to hold a single 
dislocation at a given distance s from the barrier. It follows therefore 
that 67 (using the Mott—Nabarro relation F—7rb) is the repulsive force 
exerted by the barrier on a dislocation at a distance s from it, i.e. bI is 
_ the repulsive force exerted on the first dislocation of the array above. 


ce y 
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Like edge dislocations piled up against.a barrier by an external shear stress 7. 


Since this dislocation is in equilibrium, the force exerted on it by the 
external shear stress + plus the other n—1 dislocations of the array is 
given by bI, which is equal to nbz from (22). This is a well-known result, 
first derived by Cottrell (1949) from a virtual work calculation. 

The equilibrium positions of an array of n Burgers dislocations under an 
external shear stress 7, relative to the first ‘locked’ one, have been 
calculated by Eshelby, Frank and Nabarro (1951) as the n roots of the 


equation 
2 
| al,' {=a} 0, PAGE ar ago) 


* Provided that the same estimate of 7(s) is used in both integrals. 
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where L,, is the Laguerre polynomial of order and is in B-units. These 
results, together with the preceding calculations, suffice to determine the 
equilibrium positions taken up by an array of dislocations when their 
path is blocked by some type of barrier in the lattice, i.e. the position of 
the first one is determined from (22) while the positions relative to that 
of the others are given by (23), the interaction between those and the 
barrier being insignificant. 


§5. APPROXIMATE SOLUTION OF THE INTEGRAL EQUATION 


In order to find the relative displacement 7 at any point on the slip 
direction, and in particular at the discontinuity, it is necessary to 
solve eqn. (14) for the displacement function %. When r(s) is known 
then r(s) can be found from the relation (17) without reference to 
approximations of the type (18). 

The eqn. (14) is very difficult to solve exactly, since it is non-linear 
and does not fall into any standard form. Many of the inhomogeneities, 
however, which are met with in practice (apart from Zener micro-cracks) 
constitute ‘ weak’ barriers, involving only a small change from one law 
of force to another. This suggests accordingly that we adopt a 
perturbation technique, regarding the Peierls case as a first approximation, 
since this has been solved exactly. This procedure leads to a linear 
. integral equation of which a series solution has been found in terms of — 
a certain complete set of functions. 

To effect the linearization, we write 


Pl¥]=pol¥]+pil¥, . 
h(a) =yo(w)+4(@), (24) 
where 
Pol |= —(/2c7) sin 2a, 
Po(*%)= —(1/77) tan-} a, | (25) 


these constituting the zero order solution and satisfying, from (10), 
the relation 
ree dah. (tye ade 
plda(ey=_ [et 
vo dt t— 


7 


(26) 


The quantities Pil%] and ¥,(x), where the perturbation p, is of course a 
given function of ys, are assumed to be first order, so that we may write 


PLYI=Polhot+Pilt+rilbot.] 
=Pol%ol+¥1Po [bolt+Pilbolt4p1' (vol, 


on applying Taylor’s theorem and retaining only zero and first order 
terms. The product £4p,'[%o] is also neglected as being effectively second 
order, for, although in general the derivative of a first order quantity is 
not necessarily first order, p,’ is of the same order as Py, aS can be seen 
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on physical grounds by the aid of fig. 6. Accordingly, on substituting 


P=Potp'o+p, into (14) and cancelling the zero order terms on each 
side, we obtain 


dis, at 
1 (X) Po [o(x)] +p [ypo(x) Ef ee see ec wate) 
which is a linear integral equation for the unknown function %,(~). The 
boundary conditions are . 
wy (+ c0)=(c/p), ee tOaY shes 128) 
b1(0)=0, a Meee tech eee 29) 


the former following from (5) and (24), and-the latter from the condition 
that the centre of the dislocation is to be fixed at the origin. A series 
solution for this equation with these boundary conditions is derived in § 6. 

In cases where the perturbations are not everywhere first order, it is 
necessary to be able to estimate the magnitude of the errors involved 
in the above method. On adding eqns. (26) and (27), and comparing 
with (7), it will be seen that the stress p(x) corresponding to the 
displacements 


h(x) =fo(x) +(x) 


where ,(x) is ee by (27), is simply 
=Polo(x N+ Pilpolx + 41(x)p0 [o(x) 


On Faas, x between p(x) and (x) we can find the function p(%), 
and hence p(r), for which e above (x) is the exact solution. The 
deviations from the required laws of force in each region of the slip 
plane provide a measure of the errors involved in the method ; a typical 
example is shown in fig. 11. 


§6. SERIES SOLUTION 


A solution for 4, may be obtained by expanding it as a Fourier series 
of argument 2 tan-1z, so that 


W@l=aot ZF {a,0,(2)+0,8,)  - + (80) 


n=1 
where* C,,(a)=cos (2n tan-tx) ; O,(-00)=(—1)", 
S,,()=sin (2n tan-'-x) ; S,,(- c0)=0, 
* These functions do not form an orthonormal set in the range [—%, « | 
unless multiplied by the factor 
% 1/2 
(saa) 


as follows from the properties of trigonometrical functions of argument 2 tan-z. 
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and where a, and b, are constants to be determined. The condition (28) 
is automatically satisfied by re-writing (30) as 


hb, (x)= (7/e)+ 2 {a,[C;,,(x)—(—1)"]+6,8,,(2)}, . 7 (31) 
since the terms in the summation vanish at infinity. 
Three important properties of C(x) and S,,(x), which we shall make 
use of, are as follows :— 
(1) C,,(#)—(—1)” and S,,(#) are conjugate functions in the sense 
defined by Titchmarsh (1948), i.e. they satisfy the integral relations 


ee dt 
-{  {C,()—(—1)"} — =—S,(@), 

aye ys 
1; dt Sel 


where m is any positive integer (see Appendix I). 
(2) The derivative of any function of the set may be expressed as a 
linear combination of other functions of the set, i.e 


S,'=(n/2)(Cp_r +20, +Cns1); 
C, = a (n/2)(Sp_4 25, sass J 
as may be easily verified. 


(3) The product of any two functions may also be expressed in a 
similar manner, i.e. A 


(33) 


Behm t Cnemaneels 
CG, Cer tiG a Ges). «is pee 
S,Cn= +n Spee); 

from elementary trigonometrical formulae. 


In order to expand p,[%,] in terms of the above functions, we write 
P,[~] as a Fourier series, so that 


Pilbl=p 2 {A, cos 2nmb+B, sin 2nmp}, ... . «+. (35) 
n=0 

=p 2' (—1)"{A, cos 2nmr/b+B,, sin 2nmr/b},  . . (36) 
n=0 


from (3). We note that the stress is zero when 70, i.e. when the atoms 
are in alignment, giving 
& (—1)" A,=0, Se cathe. Le: 


n=0 


and hence from (25), (35) and (37), 
Pilol=e 2 {A,[C,,(x)—(—1)"]—B, 8, (x)}. ha ane (38) 


On substituting into the integral eqn. (27) for (x) from (31), 
Pilvo(x)] from (38), and po'[Yo(x)] from (25), it is readily seen by virtue 
of the properties (32-34) that all terms in (27) can be expressed as linear 
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combinations of the functions C,,(~)—(—1)" and S,(x). On comparing 


their respective coefficients on both sides of the equation, we are led to 
the following set of relations ':— 


ne Neo , 
—2(r/)+a,+2 2' a,(—1)"=—24,, | 
[pe ear] Cee ( ) 1 ' : (39) 
NOn 11+ 2NA, + (N—2)a,_;=—2A,, (n>2), 
26,+6,=2B,, 7) 
ne 2 if (40) 
1 444+2nb,+(n—2)b,_;=2B,, (n>2). 


Nn? 


J 
These include all the coefficients a, and b,, except a,, which is found from 
(29) and (31). 

It will be seen that a, and 6, are specified by two quite independent 
sets of relations, but except for n=1 these are both second order difference 
equations of the type 

NY ni 2M Yn + (N—2)Yy 1 =Xy(NS2).  . . « (41) 


Since y, does not appear in the relation for n=2, the general solution 
contains only one arbitrary constant; this can be found from the condition 
that y,, must tend to zero as noo, as otherwise the series solution (21) is 
divergent. It may be verified (Appendix II) that the required solution 
of (41) is 
[ce] iI & ° 

Yn=—YVy © oor a. © (t= 1) 4, ¥;, (n 2), ee ee 
. Yn io Rea ae ) Gat ( a ) ( 
where Y,, is the complementary function, i.e. the solution of (41) when 
X,=0. A standard method (Boole 1946) yields the result 


L, 1 4 es ae = 
Y.-H ay — Te (n>2), 2 Dp @ (43) 


where L,,(x) is the Laguerre polynomial of order », as may be directly 
verified by substitution. 

The external shear stress 7 can be calculated directly from the relations 
(39) in the following manner. Multiplying through each relation by 
(—1)”, including the first, summing from n=1 to n= 00, and making use 
of the identity 
& n(—1)"a,4i:+ 2 (n—2)(—1)"a,1=—2 

n=2 


n=1 


leads to the result 


i MSys 


(ea 1 )( ne Nib 


N=2 


Pe A 


n=1 


=[LAG from (37) 


ne r r 
2: -|d (5) from (36). 
(ee () b 
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As the zero order law of force is an odd function, we conclude that the 
external stress may be calculated by integrating the function p(r) over 


a complete period, i.e. 
1/2 NTE 
= : = Me ee de 
[sph )aG)aaas ws 


which is in exact agreement with the result (17). 


§7. HaRMONIC PERTURBATIONS 
As a check on the theory consider a homogeneous single crystal where 
the law of force is everywhere given by 
p=(l-+e)(u/2m) sin 2n(r/b), 
there being no barrier and no external shear stress. The exact solution 
in this case is readily shown to be 


p= —(1/m) tan-! {a(1+e)}, 
= — (Im) tan} 2—(¢/m)(«]1-+2%) + 0(e2) 
from Taylor’s theorem. The above law of force, however, can be regarded 
as a perturbation of the original Peierls law with 
Pp =e(p/27) sin 277(r/b) 
so that from (36), 


By = (€/2r),.3 Ba— De 
A =A, = As 
The perturbation theory gives for this case 
Oe Dg. e 0, — Use eee 
= Ot ge eee 


so that the approximate solution is 


p= y—(e/27) sin 2 tan-t x 
= —(1/m) tan~! x—(e/m)(x/1+2), 
which agrees with the exact solution to the first order in e. 

In the simple example given above the series solution happens to 
terminate but, in general, this is not the case. Computations can be 
considerably simplified, however, by noting that owing to the linearity 
of the integral eqn. (27), the effects of different perturbations are additive. 


For this reason we have considered the displacements due to harmonic 
perturbations of the form 


e(—1)" sin 2na(r/b), 4(—1)"[cos 2na(r/b)—1], 
the corresponding #,’s being shown graphically in fig. 10. The effect 
of the general perturbation 


Pip = (—1)"{A,[cos 2n7(r/b)—1]+-B,, sin 2nz(r/b) } 
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is now simply obtained by taking the appropriate linear combination 


of the separate harmonic y,’s. A typical displacement curve, which was 
obtained in this way, is shown in fig. 11. 


Fig. 10 


2, -“ Cl)" [cos2n a - |] 


-1-0 


Graphs of :,(x) for harmonic perturbations to the Peierls law of force. 


APPENDIX I 


The pair of integral relations (32) may be most conveniently derived 
by considering the contour integral 


1--12’\" 2 dz’ 
$4 (=e) ==) jee 
where the path of integration is the real axis from —R to k(k>2) 


together with a semi-circle of radius R in the upper half-plane. The 
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(a) Displacement curve for a dislocation at one lattice spacing from a Hooke’s 
law barrier (the broken curve is the zero order solution). (b) The law of 
force. Curve (i) is the function p(r) for which the perturbation method 
yields the approximate solution (a), whilst (ii) is the function which 
corresponds exactly to (a). The latter differs from Hooke’s law by not 
more than 5°, and with negligible errors in the sinusoidal region. 
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integrand has a simple pole on this contour at z’=a (at which the Cauchy 
principal value is taken), so that from the calculus of residues 


1% n dz’ 2 1+ ix ‘ n 
§{ (SEY Cn} A 88) ar} 


The semi-circular contour C gives no contribution to the integral in the 
limit when R--0o, since 


ils 


Furthermore, since tan~!w =(1/27) log (1+-ix)/(1—ia), we have 


1 . n . 
( =) =e" 2G (nr) 468, (a). 


Lon; 


7R 
Max |R—|z| ||| 
(0 as Roo. 


Thus on separating out real and imaginary parts, it follows that 


aay dx’ 

vs eee (eet) Tt sas Be 

a “Re {C(x ) ( 1) Yon! =p BAe) 
1;2 : di’ : 
=| Same =C,,(x)—(—1) : 


which are the integral relations required. 


APPENDIX II 


We write (n—l)y,=z, and (n—1)Y,=—Z,, where y, is the complete 
solution and Y,, is the Beam function (43), so that from (41) 


Pnvit part =X; 
; Ves 
| 


Fh ote B+ )1=0. 


Multiplying the first of these equations by Z,, and the second by z,,, and 


subtracting gives 
, CAR Lie ee Oe Zn }=X An; 


VA 
ae 12, 4(7= :)} Ava n> 


where 4 is the finite difference operator (4f,=f,,1—f,,). This integrates 
to 


and hence 


Zi Ln a(7 t) "5 Kec, 


Z 
where A is a constant which is found to be zero on setting »=2(z2,=Z,=0), 
and further to 


nat 


$= 1 
Ste eet SP Geese 
Liga : Ph OME IS so 


Q 2 
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From this it follows, on substituting for z, and Z, and stepping up s and 
n by unity, that 


3 (t—1)¥,X,+-BY,. 


) il 
aoe ae 
Yn Y (s—1)s YY 41422 


s=n 
This contains one arbitrary constant and is therefore the general solution 
of (41). From (43) we find that the function Y,, oscillates and diverges 
for large values of n, so that we choose B=0, giving the solution (42). 
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XVII. The Application of Proportional y-ray Counters to the 
Determination of the Decay Scheme of 111 


By P. E. Cavanacu 
Atomic Energy Research Establishment, Harwell + 
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ABSTRACT | 


A naphthalene—anthracene crystal has been used to separate the high 
energy y-rays in I? from the main low energy component. Coincidence 
absorption measurements with f-rays and the high energy y-rays show 
that the latter are in cascade with B-rays of maximum energy 315 kev. 
Consideration of the bias curve of these coincidences however shows that 
the coincidence f-spectrum is complex and may be analysed into com- 
ponents of maximum energy 328 and 246 kev in the ratio (2 : 1) expected 
if they are in cascade with the 638 and 720 kev y-rays respectively, 
which are known to be emitted. This suggests that the latter are emitted 
in transitions to the ground state of Xel®1, 


\ 
§1. LyrRopucTION 


In recent years I'*! has been investigated by many workers, but, although 
there is a large measure of agreement as to the energy of the main f-ray, 
and the number and energy of the various y-rays, it has so far proved 
impossible to fit them together into a completely self-consistent decay 
scheme which is generally acceptable. The present paper describes 
coincidence measurements made using proportional scintillation y-ray 
counters, which it is believed resolve many of the difficulties associated 
with the establishment of such a scheme. 

Briefly recapitulating results which are generally accepted, there is 
a main f-ray of energy very close to 600 kev, and there are y-rays of 
energy 363, 283, 80 and 640 kev, (Metzger and Deutsch 1948 and others). 
There is an isomeric transition of energy 163 kev with a half-period of 
12 days (Brosi, Dewitt and Zeldes 1949). More recently there have been 
reported low intensity y-rays of energy 177 kev, (Cork et al. 1951) and 
720 kev (Zeldes, Brosi and Ketelle 1951). The latter has also been 
observed in both the primary and secondary electron spectra in this 
laboratory,* and is of about half the intensity of the 640 kev y-ray. In 
addition the last-mentioned workers have reported the existence of a 
1% branch of a B-ray of maximum energy 810 kev, also observed in this 


+ Communicated by the Author. 
* This, with other data which will be quoted in this paper, and also marked 
with an asterisk, is to be made the subject of a further paper. 
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laboratory, and possibly associated with the isomeric level in Xe. As 
for intensities, data on which have recently been tabulated by Feister 
and Curtiss (1950), the 363 kev y-ray is by far the most intense, the 
283 and 80 kev are estimated to be emitted in about 5°, of the transitions, 
while the 640 and 720 kev y-rays together are emitted in about 15% of 
the disintegrations. The other y-rays are each of about -1% intensity. 

As a result of the coincidence measurements of Metzger and Deutsch 
(1949) and the high precision measurement of the energies of the 363, 
283 and 80 kev y-rays by Lind et al. (1949), it can now be taken as 
established that the 283 and 80 key y-rays in cascade form a parallel 
branch to the 363 kev y-ray. As was to be expected the most intense 
y-ray is in cascade with the most intense B-ray. In addition Metzger 
and Deutsch showed that no appreciable fraction of the 80 kev y-rays 
could be in coincidence with B-rays of energy much less than the 600 kev 
f-ray, and this rules out the possibility that the 80 kev y-ray could be 
in one to one coincidence with one of the higher energy y-rays. That part 
of the decay scheme about which there is no dispute is indicated in fig. 5. 
The chief point of controversy is the measurement of the energy of a low 
energy B-group of about 15°, intensity, which according to Kern, Mitchell 
and Zaffarano (1949) is 250 kev and according to Metzger and Deutsch 
(1949) is 315 kev, all measurements being made on the basis of deviations 
from the straight line Fermi plot of the primary spectrum. We have 
already briefly reported the result of a measurement using a proportional 
scintillation y-ray counter in coincidence with a fast B-counter, which is 
in agreement with that of Deutsch, Cavanagh (1950). We will first 
describe this work more fully. 


§2. THE MEASUREMENT OF THE ENERGY OF THE B-RAYS IN COINCIDENCE 
WITH THE HicH ENERGY y-RAYS OF [191 


The use of organic crystal scintillators as a crude spectrometer for y-rays 
was originally reported by Jordan and Bell (1949) and this work was 
done shortly afterwards. The crystal used, of naphthalene with a small 
percentage of anthracene, took the form of a disc 2 in. in diameter and . 
3 in. thick, which was sealed onto the end window of an E.M.I. multiplier 
type VX. 5032. This was followed by a pre-amplifier and a 2 Me/s 
amplifier. The pulses were analysed by means of a simple discriminator. 
Measurements made in this way of the pulse distribution due to the 
y-rays of various isotopes such as Na*4, Au!®8, and [31 did in fact show 
structure due to the complexity of the y-ray spectra in very much the 
same way as that reported by Jordan and Bell. 

The object of the present experiment with 13! was to choose a value 
of bias such that none of the main group of 363 key y-rays was recorded, 
so that coincidence measurements with B-rays and the high energy 
y-rays alone might be made. In order to determine the point on the bias 
curve at which all the lower energy y-rays were eliminated, a series of 
lead absorption measurements in a well-collimated geometry, was made 
at various values of the bias. This set is shown in fig. 1. It is evident 
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that at low values of the bias the recorded absorption curve is complex, 
and that as the bias is increased the soft component diminishes with 
respect to the hard component, until at 20 volts bias it has completely 
disappeared. The two components appear to have exponential absorption 
curves with coefficients close to those expected for the 363, 640 and 


Fig. 1 


U2 ABSORBER (Goer 86) 


Lead absorption curves for various values of 
discriminator bias for y-rays of ['**. 


Fig. 2 


COUNTING RATE 


10 {e) 
BIAS (VOLTS) 


Bias curves for the main y-rays of 1. 


720 kev y-rays. So it is possible to extrapolate the curve for the hard 
component to zero absorber thickness, and thus determine the relative 
proportion of the two components present at any value of bias. This 
makes it possible to analyse the complex bias curve into the bias curves 
due to the two components and this is done in fig. 2. To obtain some 
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idea of the efficiency of the crystal for counting these y-rays, it was 
replaced by a copper-walled geiger counter, which gave a counting rate 
of 20 per sec as against 5 000 per sec for the crystal. The ratio of 250/1 
fits in well with the calculated efficiency of the crystal for 360 kev 
y-rays of about 50%. That for the high energy component is about 
35%, so even after the effect of the 360 kev y-ray has been completely 
suppressed, the higher energy y-ray is still counted with an intrinsic 
efficiency of 5%. 


§3. B-y CoIncIDENCE MEASUREMENTS 


In the coincidence measurements the response from the crystal was 
biased off just beyond the end point of the 360 kev component. The 
B-rays were recorded by means of an end window geiger counter operated 
at about 50 volts below the threshold, in the semi-proportional region, so 
as to eliminate dead-time, and working into a 2 Mc/s amplifier. We have 
used such counters at rates up to 50000 per sec. Both amplifiers were 
followed by discriminators, feeding in turn into a conventional coincidence 
unit. All rates were recorded simultaneously with type 200 scalers. 

An effectively weightless source was deposited on mica and sealed in 
with 4} mg/cm? nylon sheet, and its intensity was such that the genuine 
coincidence rate with no absorber was: about twice the accidental rate. 
Characteristics of the genuine coincidence rate versus delay, and also 
coincidence resolving time, were taken, and such a value of the latter, 
0-2 « sec, was chosen, that not only were all genuine coincidences recorded 
to within 1%, but some margin was left for small changes in delay caused 
by variable settings of discriminators, etc. The coincidence resolving 
time was measured at intervals throughout the experiment at the single 
rates actually used, and was found to remain constant to within the 
measured accuracy of 1%. 

The source was placed almost in contact with the crystal and a short 
distance away from the B-counter, so as to permit of the introduction of 
aluminium absorbers. Coincidences were measured as function of thickness 
of absorber placed before the B-counter, at least 4 000 coincidences being 
recorded at each point. The genuine coincidence rate decreased by a 
factor of over a hundred over the range of absorber thickness used, being 
twice the accidental rate at zero thickness and falling to about one-third 
of the accidental rate at the maximum thickness of absorber used. 
Corrections were made for the cosmic coincidence rate and also for 
y-y and scattered coincidences. Even at the largest thickness of absorber 
used these corrections were small. The absorption curve of the coincident 
B-rays was derived in the following way: at each absorber value the 
ratio genuine/accidental coincidence rate was calculated, with the genuine 
B-y rate corrected as described above. The source was replaced by a 
much smaller one, and the absorption curve of the primary f-rays 
measured in the same geometry as the ‘coincidence experiment. Each 
value of the genuine/accidental rate was then multiplied by the counting 


es 


aE 
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rate obtained for the same thickness of absorber in the subsidiary 
experiment, corrected for backgrounds. This gives a quantity proportional 
to the detection efficiency of the B-rays in coincidence with the high 
energy y-rays at each value of absorber. The reason for this procedure 
is that it 1s necessary to eliminate a change in intrinsic efficiency with 
counting rate in the semi-proportional B-counter, which amounted to 
some 20% over the range of counting rates used. The absorption curve 
derived in this way is shown in fig. 3. It is essentially exponential over 
the range of measurement. A comparison was made with the absorption 
curve of Co® £-rays, measured in the same geometry, under the same 
back-scattering conditions, and at small counting rates. It is interesting 
to note that the estimated errors in this measurement are comparable 
with those of the coincidence measurement, arising in this case from the 
large and somewhat uncertain value of the y-ray background. 


Fig. 3 
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(a) Absorption curve for Co® B-rays (310 kev). 

(b) and (c) Absorption curves for 328 and 246 kev f-rays.derived from (qa). 

(d) Experimental points from f—y coincidences for high energy y-rays of I* 
together with composite absorption curve of 328 and 246 kev f-rays. 


Inspection shows that the energies of the two f-rays are very nearly 
equal, and a Feather analysis gives the value for the maximum energy of 
the f-rays in coincidence with the high energy y-rays as 315 kev, if the 
maximum energy of Co® B-rays is taken as 310 kev. This is in agreement 
with the result of Metzger and Deutsch (1949) and the estimated error 
of ~5% is such that it is believed to be incompatible with the value 
of 250 kev derived by other workers. On the basis that the 603 kev 
B-ray and the 363 kev y-ray are in cascade, the total energy of the 
transition is 966 kev; the energies of the two high energy y-rays as 
measured in this laboratory * are 638 and 720 kev. The coincidence 
measurements of Deutsch require that the 80 kev y-ray is not directly 
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in cascade with a higher energy y-ray, hence it is very probable that 
these two y-rays are emitted in transitions from levels at 638 and 720 kev 
above the ground state of Xe!*! and are associated with B-rays of maximum 
energy 328 and 246 kev respectively. An energy of 328 kev is undoubtedly 
compatible with the measured energy of 315 kev, but the 720 kev 
y-ray is certainly in coincidence with a f-ray (the small conversion 
coefficient of 0-27°/ ensures that it does not belong to an isomeric 
transition) and this B-ray must be of energy 246 kev, or less (in the case 
where a further y-ray is in cascade with the 720 kev y-ray). As will be 
shown below, by consideration of the bias curve of B-high energy y-ray 
coincidences, it must in fact be in cascade with a f-ray of energy not 
very much less than 250 kev. The low energy f-group should then 
consist of two f-spectra of energy 328 and 246 kev in proportion to the — 
relative intensities of the 638 and 720 kev y-rays, with a correction made 
for the difference in detection efficiency for these two y-rays. Absorption 
curves for the two B-spectra were constructed from the experimental 
one for Co® B-rays, using the Feather rule, and taking the maximum 
ranges of 246, 310 and 328 kev f-rays from interpolations in the table 
given by Gupta and Ghosh (1946). An estimate of the relative detection 
efficiencies of the 638 and 720 kev y-rays was made by first comparing 
the bias curves for the soft and hard components, and from this, estimating 
the relative bias curves for the two hard components. The relative 
intrinsic efficiency of the crystal for the latter (at zero bias) was calculated 
on the basis of Compton absorption. This together with the spectrometer 
measurement of a relative intensity of 2: 1 gives the ratio of the mixture 
of the 328 and 246 kev f-spectra as 1-6: 1 at the value of bias used for 
the coincidence measurement. These two absorption curves, together 
with the composite, constructed on the foregoing basis, are also 
reproduced in fig. 3. The energy derived from the latter is 310 key, 
in excellent agreement with the experimental value. 


$4. Bras CuRVE oF B~y COINCIDENCES 


Further information may be derived from an examination of the bias 
curve of B-y coincidences, in this instance with no absorber inserted 
before the 6-counter but only the zero value of 9-0 mg/em? aluminium 
equivalent due to the counter window, the source covering, and the air 
between source and B-counter. If we consider the ratio of the genuine 
coincidence rate to the product of the single counting rates, then at 
low values of bias the great majority of the coincidences will be due to 
363 kev y-rays and 603 kev B-rays, the latter not being greatly attenuated 
in the zero thickness of absorber. At high values of bias however, 
coincidences between high energy y-rays and low energy b-rays will 
become increasingly predominant. The latter are heavily attenuated in 
the zero absorber thickness, and hence the recorded y-rays will produce 
relatively fewer coincidences. Hence one would expect the ratio of 
genuine coincidences to the product of the single rates, as the bias is 
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increased from small values, to fall slowly over the region where 
coincidences due to the main group of low energy y-rays are predominant ; 
ee Seales the intermediate region, and to level off beyond the 
pomt where the coincidences due to the high - 
eran gh energy y-rays become 
Such a characteristic was in fact obtained and is shown in fig. 4, all 
quantities having been corrected for backgrounds. It becomes important 
to know if this measured characteristic is in fact compatible with the 


mixture of low energy f-spectra required to be in coincidence with the 
638 and 720 kev y-rays. 


Fig. 4 
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(a) Experimental points and composite bias curve for 640 and 720 kev y-rays 
in coincidence with 328 and 246 kev f-rays respectively. 

(b) Calculated bias curve for y-ray of 670 kev in cascade with B-ray of 315 kev. 

(c) Calculated bias curve as (a) except low energy f-ray 210 kev. 


At this point it becomes necessary to consider an aspect of the bias 
curves due to y-rays of different energies which we have only touched 
on very briefly before. If we consider those values of bias for which 
the initial counting rate falls by successive factors of 2, then the ratio 
of corresponding bias values for y-rays of different energies is constant 
over a range of counting rates of at least one hundred, and is approximately 
equal to the ratio of the y-ray energies. This property 1s unexpected 
because the resolution of the crystal would certainly be expected to 
depend on the energy at least as far as electrons are concerned. However | 
the absorption characteristics of y-rays, including secondary processes, 
certainly enter into the question, and make it extremely difficult to handle 
from the theoretical view point. We have therefore been content to take 
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this as an empirical property, which has proved extremely useful in 
separating the components of a complex spectrum, and estimating y-ray 
energies. : ; 

We may illustrate this property by reference to the bias curve | of 
y-rays in the experiment under consideration. The relative counting 
rates of the two components is taken to be that determined in the original 
absorption measurement ; the relevant bias values are given in the 
following table. 


Table 1. Values of Bias Required to Reduce Initial Counting Rate 
by Successive Factors of Two 


Bias value for 
low energy | 4:1 6-8 8:8° 110-07 (AL-1 )12-0 
component 


Bias value ~ for 
higher energy | 8-2 | 13-7 |16-8 | 19-2 | 21:6 | 23-5 
component 


Ratio high/low 2-00 | 2-00} 1-91] 1-92} 1-95 | 1-96 |average—1-96 


The bias curve of the high energy component should correspond to a 
mean energy of 670 kev, and therefore its ratio to the energy of the 
363 kev y-ray is 1-85, which compares with the 1-95 derived from the ratio 
of bias values. For our present purposes however, it is sufficient to take 
it as strictly proportional and to construct, on this basis, bias curves 
for the two high energy components separately. We may also from fig. 3 
read off the relative efficiencies of detection of the 328 and 246 kev 
B-rays for 9-0 mg/cm? aluminium absorber, and from an experimental 
curve determine the attenuation of the main 603 key component in the 
same absorber. The ratio of the intensities of the two high energy 
components is taken, as previously, to be 2:1. All this information 


enables us to construct a curve (a) corresponding to the experimental | 


points of fig. 4; and the fit is extremely good. On the other hand, a 
similar bias curve (b) calculated on the assumption of a single y-ray of 
energy 670 kev in coincidence with a B-ray of 315 kev is not compatible 
with the experimental points, the B-rays being attenuated faster in the 
first 9-0 mg/cm? of absorber than is appropriate to an energy of 315 kev. 
Taken in conjunction with the previous absorption measurements this 
implies that the low energy B-spectrum is in fact complex, and it may be 
resolved into components of energy 328 and 246 kev in the proportion 
expected if they are in coincidence respectively with the 638 and 720 kev 
y-rays. If the lower energy component is assumed to have an energy 
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as little different as 210 kev, curve (c) in fig. 4 results, which again is 
inconsistent with the experimental results. This strongly suggests that 
the 720 kev as well as the 640 kev y-rays correspond to ground state 
transitions. 

The decay scheme shown in fig. 5 also includes this new evidence. If it 
is assumed that the known high energy B-spectrum leads to the 163 kev 
isomeric level, the only known y-ray not included in the scheme is the 
177 kev y-ray. 

These results are not in agreement with the published decay scheme of 
Bell et al. (1951) but later work by the latter group is essentially in 
agreement with us (private communication). 


Fig. 5 
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Decay scheme for I'*!. The part enclosed within the 
broken lines is already fairly well agreed. 


§5. CONCLUSIONS 


B-y coincidence absorption measurements have shown that the high 
energy y-rays in [1*! are in coincidence with f-rays of maximum energy 
315 kev. A consideration of the bias curve of these coincidences, however, 
shows that the B-spectrum is complex and may be analysed into components 
of maximum energy 328 and 246 key in the ratio (2: 1), expected if they 
are in cascade with the 638 and 720 kev y-rays respectively. This suggests 
that the latter are emitted in transitions to the ground state of Xe?*!. 
These results are compatible with the latest measurements of P. R. Bell 
using a sodium—iodide crystal spectrometer (private communication). 
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ABSTRACT 


A study has been made of the nuclear disintegrations produced by z- 
mesons of energy between 50 and 1100 mev, and they have been compared 
with those produced by protons of energy up to 800 Mev. The results are 
consistent with the hypothesis that a 7-meson commonly interacts with a 
complex nucleus by making two or more elastic collisions with its nucleons. 
By studying the characteristics of the meson ‘ stars ’, estimates have been 
obtained of the relative probabilities of the meson being absorbed in a 
nuclear collision, of undergoing charge exchange, or of being elastically 
scattered. 


$1. [LvTRoDUCTION 


THE application of the photographic emulsion technique to cosmic-ray 
studies and the construction of the great synchrocyclotrons, have made 
possible detailed investigations of nuclear reactions at energies above 
the threshold for meson production. Considerable progress has thus been 
made in our knowledge of high-energy nuclear processes, and new problems 
have emerged, of which one of the most important is the following : 

When a nucleon with an energy above | Bev collides with a nucleus, a 
large fraction of the available energy is absorbed in the creation of 7-mesons. 
These mesons interact strongly with the nucleons, so that they may be 
involved in secondary interactions within the parent nucleus and 
contribute to its disintegration. How important are such interactions and 
what are their detailed characteristics ’ 

This problem may be approached. in a number of ways, one of which is 
to study the disintegrations produced by z-mesons which have escaped 
from the parent nucleus. Such studies have been made by Bernardini 
et al. (1950, 1951) using artificially generated 7-mesons of energy less than 
110 mev. Although similar experiments using 7-mesons in the cosmic 
radiation are more arduous, they allow the study to be extended to particles 
of much greater energy and this paper describes preliminary results in the 
interval up to 1-1 Bev. In several instances the statistical weight of the 
observations is small, but the analysis has been made in order to indicate 
the method of approach which may be adopted. 


* On leave from the Hebrew University, Jerusalem. 
+ Communicated by Professor C. F. Powell, F.R.S. 
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§2. IDENTIFICATION AND MEASUREMENTS ON THE STARS PRODUCED BY 
Mrsons AND PROTONS 


In the matter of photographic plates exposed to the cosmic radiation, 
nuclear explosions are produced which are accompanied by the creation of 
a-mesons. Many of these particles escape from the assembly, but some 
make secondary nuclear disintegrations within it. The problem of 
recognizing these particular events can be solved in one of two ways : 

It sometimes happens that a nuclear collision involving the production of 
many ‘ shower particles’ occurs in the emulsion, and that one of these 
‘shower particles ’ makes a second collision in the same plate. These 
events are referred to as ‘ double stars’. It is possible to identify the 
particle linking the two disintegrations as a 7-meson if the track is 
sufficiently long and if the mean angle of multiple scattering, ~, is greater 
than 0-0226° per 1004. In such events the direction of motion of the 
a-meson is generally unambiguous—the probability is very small that 
the particle assumed to produce the second disintegration really emerged 
from it—but they are rare. Thus, in an examination of 125 cm? of 
emulsion, the observers in this laboratory have recorded 30 000 stars, from 
which a total of 14 800 shower particles were ejected. Of these, only 40 
were observed to make a nuclear interaction before leaving the emulsion, 
and only 15 of them were identified as 7-mesons by scattering and grain 
density measurements. These numbers are so small that we cannot 
expect to obtain observations of sufficient statistical weight by this method 
alone. A second method was therefore employed. 

In the second method, each isolated star recorded by the observers was 
examined for associated tracks of minimum ionization. Scattering and 
grain-density measurements were made on those tracks which were (a) in 
the ‘ upper hemisphere ’, and (b) of sufficient length in the emulsion to 
allow « to be determined. From the tracks thus selected, only those with a 
mean angle of scattering, z, >0-0226° per 100, were accepted for further 
study, i.e. only those which gave sufficient evidence to allow the nature and 
energy of the primary particle to be determined. (Fowler 1950.) A 
proton or 7-meson so identified was regarded as having produced the star 
with which it was associated if it possessed sufficient energy to cause the 
disintegration.* 

Results of the scattering and grain-density measurements, for these 
‘primary’ particles and for those identified from the double stars, are 
represented in fig. 1. The resolution between protons and mesons is 
satisfactory. The identification was ambiguous in only two cases, and 
they were excluded from the observations. It will be convenient to 
refer to the ‘stars’ produced by the two types of particle as ‘ meson 
disintegrations ’ and * proton disintegrations ’ respectively. 

The value of the scattering parameter <=0-0226° per 100 u corresponds 
to a kinetic energy of 800 Mey for a proton, and 1100 Mev for a z-meson. 
The investigation was therefore confined to particles with energy less than 


* Since none of the emulsions searched were exposed under dense materials, 
the number of stars produced by u-mesons is very small. 
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these values. The numbers of disintegrations available for analysis, 
produced by protons and mesons, were 70 and 86, respectively. 

The second of the methods described above, will lead, in a small pro- 
portion of the events, to errors of identification; and some of the ‘ primary ’ 
particles will, in fact, be secondary particles which originate in stars 
produced by fast neutrons. In such events there is a certain probability 
of a shower particle being ejected ‘ backwards ’; and, although rarely, it 
may satisfy the criteria which lead us to regard it as a ‘ primary ’ particle. 
The proportion of such cases of wrong identification has been estimated by 
Camerini ef al. (1951, VI), and shown to be less than 4) Oe 


Fig. 1 
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Distribution in the values of x° (100 .) and grain density of mesons and protons 
identified as the primary particles of nuclear disintegrations. 


This second method also suffers from the disadvantage that the observers 
only record stars with N,>3, and they miss a substantial fraction of the 
stars with N,,—3 and 4 (Rochester and Rosser 1951). Some of the smaller 
stars produced by 7-mesons therefore escape observation. ‘The disintegra- 
tions observed by the first method are, however, recorded with an efficiency 
which is independent of N,. A comparison between the results of the two 
methods indicates that about 70° of the total number of meson stars 
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occurring in our emulsions have been observed, the missing 30°, corre- 
sponding to about half the stars with V,,<4. 
Measurements on ‘ Stars’ with Identified Primary Particles 

The stars produced by protons of energy less than 800 Mev or mesons of 
less than 1100 Mev, were analysed in the following manner. The length of 
each track associated with the star, its grain-density and angle in space 
with respect to the primary particle were measured. In favourable cases, 
scattering measurements were also made on the secondary particles. In 
addition, any particles which reached the end of their range in the emulsion. 
were identified by inspection. 

From these observations, the variation of a number of parameters 
characteristic of the disintegrations could be determined as a function of 
the primary energy, for stars produced by both mesons and protons. Fig. 2 


Fig. 2 
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0G Results of Germain (1951) for artificially accelerated protons. 
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shows the variation of the average number of heavily ionizing branches 
N ,,; fig. 3 the average number of ‘ grey ’ tracks N ervey} and fig. 4 the average 
energy of the ‘ grey’ particles, Eyoy. 

Values of the mean-free-path for interaction of both z-mesons and 
protons were obtained from an analysis of the ‘ double stars ’, the calcula- 
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tions being made by the method described by Camerini et al. (1950, IV). 
The results are shown in tables 1 and 2, together with data obtained by 
other workers. The results for mesons are consistent with a geometrical 
cross-section for the nuclei of the emulsion. Some results for shower 
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Variation with primary energy of the average number of ‘ grey’ tracks, Nyroy, 
from meson ( x ) and proton (@) disintegrations. 
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particles obtained subsequently to those given by Camerini et al. (1950, IV) 


are summarized in table 3. 


Table 1. 


Energy 


range Mesons 


pes 
AS 
— 
100-110" =a 
30-40 


7 


m+ 


FS) 
| Ht 


{fre PS ald ot ete | 


YAMIIIIAAAA 


80-1100 


3 
H- 


Table 2. 
Stars 


N,21 

N,>3 
Identified 

Protons 


SNnzl 
VY, 23 


They are consistent with an interaction 


The Mean Free Path for Interaction of 7-Mesons 
M.F.P. in 
gem? a Technique Authors 
of material Geom. value 
95-+12, emulsion 1:5 Emulsion ae 
78+10, emulsion 1-0 | ( Bernardini, 
86-117, emulsion 1-1 > do. Booth and 
72-- 7, emulsion 0-9 J Lederman (1951) 
: ; Jf Bradner and 
86+ 9, emulsion 0-93 do. Rankin (1950) 
Cloud Camac et al. 
aU A { chamber (1951 a) 
~ 100, C 7 do. fae. (1951 b 
125-410, H 0-22 
48-42: Ti 1-08 
59+ 4, Be 0-96 
Boa 2 C 1-08 i Te een 
570 1-2 ‘ Counters d Stel : e 
724+ 3, Al 1-13 an bee, erger 
106-£ 5, Cu 1-0 Naeet 
117.5, Cd 1-12 
143-L 6, Pb 1-2 
‘ . : . Lock and 
95-+28, emulsion 1-05 Emulsion Yekutieli (1952) 


The Mean Free Path of Profane for Star Production 


Energy « _ 
range E 
(Mev) 

40-100 60 

40-100 60 

100-500 180 

100-500-175 

100-500-195 

100-500-195 

240 

— 175 
— 245 
— 340 
350-400 375 


Path 
length 
(cm) 


Estimated M.F.P. 


number of cm of Authors 
interactions emulsion 
8°75 35-+13 
5:75 54420 
10 30-L10 Lock and 
6 51+20 ( Yekutieli (195 
21 36-+ 8 
16 47+10 | 
281 82+12 Perry (195] 
91 160+17 
90 100+11 } Germain (1 
206 85+ 6 
Bernardini 
Booth, anc 
ge bo. 9 Lindenbau1 
(1951) 


The second group of Bristol results, lines c and d, refer to double stars caused by ident 
protons of track length greater than 3 000 4, while the first and third groups refer to all do 
stars with an interconnecting track greater than 400 in length and of grain density wi 


specified limits. 


The number of interactions due to protons has been estimated from the obse. 


fraction of protons amongst the identified interconnecting particles of that grain density inte1 
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Sol as for shower particles between 0-6 and 1-0 times the geometrical 
7 O . 

va me e 7-mesons, which comprise some 80°, of the shower particles, 
must have an interaction cross-section of the same order of magnitude. 


§3. CHARACTERISTICS OF THE PROTON DISINTEGRATIONS 


Several studies have been made of the characteristics of ‘ stars ’ produced 
by protons of energy up to 385 Mev (Bernardini, Booth and Lindenbaum 
1951, Germain 1951, Perry 1951). The general features of the disinte- 
grations appear to be in satisfactory agreement with the theories of 
Goldberger (1948) and Wouthuysen (1949). These theories assume that 
with protons of energy below 1 Bev, the majority of the collisions are 
elastic in the sense that no energy is released in the form of mesons. 
Because of this feature it has been convenient to compare and contrast the 
meson induced stars with those due to protons. 


Table 3. Some Recent Determinations of the Interaction Mean Free 
Path of Shower Particles ; 


M.F.P. in ‘ Ener 
Method Authors 2 cm or reese ee 
material Geom. value (mesons) 
Cloud chamber at 
mountain alti- Hartzler 230+60, Au 1°5 >1000 Mev 
tudes 
Gregory and |172-+30, Pb te 3 
Tinlot  |164-50, Al 2-0 Ree 
Lovati, Mura, 150 | 
Succi and |200+50, Pb 1-25 mA aye 
Tagliaferri ee 
. Photographic ; ; 
emulsion ae ok ns 0, 1-4. >100 mev | 
technique ekutieli emulsion 


(a) Meson Production by Fast Protons 

In 60 of the 70 disintegrations due to protons the energy of the particle 
was greater than the threshold for meson production, viz. 170 Mev. The 
ejection of a charged 7-meson was observed in only two cases. Applying 
a correction for those mesons which escape identification,* it is found that 
only 0-067 +-0-045 charged mesons are emitted per disintegration. 

In an attempt to estimate the total number of mesons produced, both 
charged and uncharged, a study was made of the ‘energy balance’ in 
these disintegrations, for different values of the primary energy #.0 be 
average values of H,, V,, N,, and Ngrey were computed for each group, 


pee et ies oe SE eee 
* For details see § 5 (i). 
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together with the average energy for each of these types of particle. V,, 
represents the number of ejected protons, and NV, the number of ejected 
z-particles of energy less than 25 Mev, i.e. those ejected during the 
evaporation process. In each group, except N,, allowance was made for 
the emitted neutrons. The results for different values of H, are shown in 
table 4. Although subject to large statistical erorrs, they strongly suggest 
that in the energy-range under consideration the fraction of the available 
energy which appears in the form of mesons is very small. It is, of course, 
possible that mesons are created and subsequently absorbed within the 
parent nucleus. 


(b) Recoil Protons 


The tracks of 52 secondary particles with a grain-density less than 80 
grains per 50. were found associated with the 70 stars due to protons. This 
value of the grain-density corresponds to a proton of 25 Mev, a deuteron of 
50 Mev and a triton of 75 Mev, respectively. Scattering and grain-density 
measurements showed that 2 of these particles were mesons, 17 were 
protons, 4 were deuterons and 1 was a triton. The others could not be 
identified. 


Table 4, The ‘ Energy Balance’ for Proton Disintegrations : 
Average energy of Average energy 
primary particle in star 
(Mev) (Mev) 

123415 133438 
286+11 204441 
500 +12 410 +50* 
681 +16 672+188 


* Allowance has been made for the two z-mesons observed to be ejected. 


Let us assume that the primary proton, p, makes a succession of elastic 
collisions with individual nucleons in the nucleus. A cascade of recoiling 
nucleons will then be produced in the way illustrated schematically in fig. 5. 
For each ejected particle there is a trajectory through the nucleus which is 
inked with the primary collision. For example, particle 4 is ejected as a 
result of collisions a, 6 and ¢, i.e. after 3 collisions. 

Suppose that the energy of the primary proton, H,, the energy of each 
ejected charged nucleon, #,,, and the angle, «, between the direction of 
motion of the primary particle and that of the ejected nucleon are all 
known. If a nucleon is ejected after only one collision (e.g. particle 1), 
the relationship between £#,, H,, and « is eh 


cos a= 4/{(1+ 2u/H,)/(1+2u/E,,)}, cet eS 
where ee is the proton rest-mass. If, in any particular case, the observed 
value of « is greater than that derived from the above expression, using the 
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measured values of E,, and E,,, the particle must have been ejected as a 
result of more than one nucleon-nucleon collision. If the number of 
collisions leading to the ejection of a particle is m, the maximum value of « 
is obtained when the magnitude, 8, of each of the m individual collisions is 
the same in each case, and equal to 


cos” B= (1+ 2u/E,,)/(1+-2u/E,,) ee es eee) 
and 
eee 00 aes Cy ers 3) 
Fig. 5 


Schematic diagram of a nucleon initiated cascade in a nucleus. 


Inserting the measured values of H,, H,, and 0(=%max), in Geb eae o 
and (3) we can find the minimum value of m, i.e. the minimum num eo . 
collisions required to cause the ejection of the observed particle with the 
2 d values of H,, and «. : : 
Sree of the disintegrations from which 17 identified sae abe 
emerged showed that in 3 cases the proton could be regarded as ee 
made only one collision in the nucleus (i.e. relation (1) was satisfied), ae 
in the remaining 14 cases at least two collisions took place. Similar 
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calculations have been made for the 28 unidentified secondary particles, 
assuming them all to be protons, and the results indicate a minimum of 
one, two and three collisions within the nucleus in 4, 17 and 7 cases, 
respectively. The average number of internal collisions—as deduced 
from all the 45 secondary particles—was 2-0, a value which corresponds 
to a cross-section for the interaction between two nucleons of 3 x 10~?° cm?. 


(c) Angular Distribution of ‘ Grey’ Particles 
The angular distribution of the ‘grey’ tracks associated with the stars 


produced by protons is shown in fig. 7. The median angle of the distri- 
bution is ~40°. For 240 mev protons, Perry (1951) finds that the median 


Fig. 6 
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angle of the “ grey’ particles is 30°. These figures are in qualitative 
agreement with the theoretical predictions of Goldberger (1948) and 
Wouthuysen (1949), (See also Bernardini, Booth and Lindenbaum 1951.) 


§4. CHARACTERISTICS OF MEson DISINTEGRATIONS 


The simplest model for the production of a nuclear disintegration by a 
m-meson has been proposed by Heitler and Janossy (1950). The incident 
meson ‘knocks-on’ a nucleon in the struck nucleus, and this recoil nucleon 
causes the observed disintegration by secondary collisions within the 
nucleus. We have attempted to test the validity of this model, 


Interaction of Fast 7-Mesons with Nuclei 241 


oeic , haoiyy ‘ 
cool ta A a lel ea 
(ii) by studying the ejected protons ; and 
(iii) by studying the ejected mesons. 
(i) Variation of N,, with E,. 
if the meson stars are indeed produced by a single recoiling nucleon 
there should be a strong similarity between the stars produced by neh are 


a ; , 
nd by mesons. The maximum energy which a meson of total energy, 
w,, can transfer to a nucleon at rest is 


Enax= 2M )(w,?—moc*)/{mec?+ Mo(2w,,+ Myc?)} be Pee poe ql) 
(Janossy, 1950), 


Fig. 7 
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where m, and M, represent the rest-masses of the 7-meson and proton, 
respectively. The nuclear disintegration produced by a recoil-nucleon of 
this energy will, on the average, be a little less energetic than one produced 
by a proton of the same energy incident from outside the nucleus, for the 
recoiling nucleon has to traverse a smaller thickness of nuclear matter. 

- For each disintegration produced by an identified meson, the maximum 
energy E,,,, which could be transferred to a recoil nucleon was computed 
from equation (4), and in fig. 6 the size distribution (N,) of these stars is 
plotted as a function of H,,,,. The distribution for disintegrations produced 
by protons is included in the same figure. There appears to bea significant 
difference between the two distributions, more energy being transferred to 
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the nucleus in the meson stars than is compatible with the assumption of a 
single recoiling nucleon. This result is particularly striking since we have 
assumed the recoil nucleon to receive the maximum and not the average 
energy. For low energies (w,<450 Mev), the energy transfer is lar ge 
compared with that of a proton of the same energy (see fig. 2) and it is 
reasonable to assume that a considerable fraction of a mesons are 
completely absorbed in the nucleus. (See also §5, iii.) 

(ii) Recoil Protons 

69 particles of which the grain density in the track was less than SO 
grains per 50. were emitted from 66 stars produced by 7-mesons of energy 
greater than 160 Mev. 23 of these particles were protons, 2 were deuterons, 
and the others could not be identified. 

If we assume that the incident meson makes only one collision with a 
nucleon when it interacts with a nucleus, then all the ejected nucleons will 
arise from the nucleon ‘ cascade’ produced by the recoil particle. This 
cascade will be similar to that initiated by each of the recoiling particles 
arising from the first nucleon-nucleon collision in the ‘disintegrations 
produced by protons. 

Now the number of collisions which link the ejection of a charged 
nucleon with the first meson—nucleon collision may be calculated in a way 
similar to that described in $3(b). It is thus found that the average 
number is about two. On the other hand, in disintegrations due to 
protons, the first generation of recoil particles were found to make, on the 
average, only one collision before escaping from the nucleus. It is 
reasonable to suppose that the physical processes responsible for the 
nucleon cascades, in the disintegrations produced by protons and mesons, 
are similar. The above result therefore appears to be inconsistent with 
the assumption of a single meson—nucleon collision in the disintegrations 
produced by mesons. The results suggest that only about 30% of the 
emerging nucleons are ejected as a consequence of the first meson—nucleon > 
interaction. 

Let A be the average number of collisions which the incident meson 
makes in traversing the nucleus. The probability that the meson makes 
one or more collisions is (l1—e~’), and the fraction of the total number of 
meson collisions which correspond to the first meson—nucleon interaction 
is (l—e~*)/A. Since the results given above show that about 30% of the 
meson—nucleon collisions correspond to these encounters, it follows that the 
average number of collisions which the incident meson makes in traversing 
the nucleus is about three. 


(iii) Angular Distribution of the * Grey’ Particles 

Fig. 7 shows the angular distribution of the ‘ grey’ tracks associated 
with the meson disintegrations. Any tracks identified as due to mesons 
have been excluded. The distribution is much broader than that for 


proton disintegrations, and a considerable number of the particles are 
ejected at an angle of more than 90° to the primary meson. The fraction of 
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such particles is called the ‘ backward fraction ’, B,, and its observed 
value is 0-29+0-05. The median angle of the distribution is 60°. It is 
possible that some of the particles ejected backwards are unidentified 
7-mesons, but their number cannot be very great, for the value of B. for 
the identified protons amongst the ‘ grey’ particles is 0-21 -40-10. "The 
ejection of these particles at a large angle to the direction of the primary 
meson cannot be explained in terms of a single meson-nucleon collision. 
(iv) Summary 

The results show that the meson disintegrations cannot be interpreted 
in terms of a single nucleon, recoiling from the incident meson, and 
producing a nucleon cascade. The average number of meson—nucleon 
collisions within the nucleus, for nuclei of average atomic number, appears 
to be between two and three. In collisions with nuclei, mesons of low 
kinetic energy, <160 Mev, have a high probability of being absorbed. 


§5. THe Fats or tHE [NcIDENT Meson 


As a result of its interaction with a nucleus an incident meson may 
(i) emerge as a charged meson, having lost some energy ; 
(ii) emerge as a neutral meson, having suffered charge exchange ; 
or (ili) be completely absorbed, the energy corresponding to its rest- 
mass contributing to the disintegration of the nucleus. 
(i) Ejection of a Charged Meson 

In 16 cases a single charged 7-meson has been observed to emerge from a 
meson-disintegration, and the results are summarized in table 5. In 
most of these events the energy of the outgoing meson is much less than 
that of the primary particle. There appears to be no simple correlation 
between the loss of energy and the angle of scattering of the meson.* A 
similar observation has been reported by Bradner and Rankin (1951) for 
38 Mev 7 -mesons, and by Bernardini et al. (1951a, b) for 85 and 110 Mey 
7 -mesons. 

Energy and momentum considerations show that the first three cases 
given in table 5 may be attributed to a single scattering of the primary 
meson by a nucleon. Similarly, the fourth may be attributed to the 
scattering of a meson by an «-particle. In each instance, the kinetic 
energy of the primary meson was less than 350 Mev. 

Assuming that a 7-meson can only lose energy in elastic meson—nucleon 
collisions, the observed energy loss in the remaining 12 cases can only be 
accounted for if it is assumed that the average number of collisions, made 
by the meson’ with nucleons, is about three. The minimum number of 
collisions in each individual case is given in table 5. : 

The 16 cases given in table 5 represent favourable examples, in which, 
by chance, the track length in the emulsion was sufficiently great to permit. 
identification. An estimate of the fraction of ejected mesons which are 


not identified has been made in the following way: of tracks with 
OES? CER AP tie aR A Ai a TON Aidaaiale ete BIRRIE NE AES Te 
* For plates see end of issue. 
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grain-density less than 80 grains per 50, a certain number have been 
identified as mesons, protons, deuterons or tritons. It is reasonable to 
assume that the proportions of different particles amongst those not 
identified are the same as for the identified particles. An approximate 


Table 5. The Scattering of Energetic 7—Mesons 


Kinetic Kinetic 
energy of | energy of : ‘ Minimum 
Number | primary | secondary pes : of Eat) » | NV, | number of 
mesont mesont Mee REPTGN aaa ee collisions 
(Mev) (Mev) (mev) (angle) 

1-07 79 67 138° — — 2 1 
1-19 50 39 63° —- — 1 1 
2-07 252 180 62° —- — 4 1 
2°16 340 28 28° —_- — 1 1 
3-03 460 310 26° —- — 1 1 (2) 
3-14 640 39 18° 380 Yaad Cad 3 
3:15 375 85 19° —- — 6 2 
4-07 755 70 167° — 6 2 
4-11 755 450 = LD 43°) AS aha 1 (2) 
4:17 825 17 125° 100482 tro 4 
4-19 670 110 22° —- — 9 2 
4-21 785 56 96° —_- — 6 3 
5-01 900 29 43° 45 104° | 6 3 
5-08 1040 120 5° 52 20°; 3 ae 
5-11 960 225 eis —- — 3 2 
5-12 1050 120 90° fs 4 2 


* In this column the energy of the ‘ grey ’ particle is given first, assuming it to 
be a proton, followed by the angle which it makes with the primary particle. 
+ The error in the values of the energy is <20%. 


Table 6. The Frequency of jection of z~Mesons from Meson 


Disintegrations 
Primary No. No. of | — . 
Group energy of mesons | "” i aie 
(Mev) stars observed Cres 
I 0-160 20 2 0-265 +0-115 
i 160-360 18 2 0-21 +0-15 
inal 360-660 15 3 0-43 -+-0-17 
IV 660-860 21 5 0-55 +-0-16 
V 860-1060 11 4 0-61 40-25 
II-V 160-1060 65 14 0-44 +0-12 
All 0—1060 85 16 0-39 -+40-10 


value can then be obtained for the total number of charged mesons, 
n,, ejected from meson disintegrations, for different values of the primary 


energy. The results are shown in table 6. The value of 7”, for all the 
stars is 0-39 -L0-10, ; 
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{ii) Conversion 


Tt follows from the above result that ~60°, of the incident mesons are 
_ either absorbed in the target nucleus, or emerge from it having lost their 
charge. Evidence for the relative importance of these two processes can 
be obtained by studying the energy-balance in the meson disintegrations. 
If the average energy associated with a sample of such disintegrations is 
greater than the average kinetic energy of the primary particles, # this is 
evidence for the absorption of mesons; if it is less, it is reasonable to attri- 
bute the © missing energy ’ to charge-exchange, for the neutral 7-mesons 
commonly escape detection ; and if there is a balance of energy it may be 
that both absorption and charge-exchange occur, the energy liberated by 
absorption compensating that ‘lost ’ through conversion. 

The study of energy-balance has been made by the same methods 
employed for the proton stars. Allowance has been made, however, for 
the ejection of charged ‘mesons, using the results given in table 6. The 
results are set out in table 7, from which it may be seen that for mesons of 


Table 7. The © Energy Balance‘ for Meson Disintegrations 


Average energy of | Average energy 


primary particle in star 
(Mev) (Mev) 
95+ 7:5 129+ 22 
270-+12 206+ 37 
493 +20 568 +190 
761+13 605+ 91 
1000 +15 7714183 


energy less than 160 Mev there is some evidence for absorption. This result 
will be considered in more detail in the next section. For mesons of energy 
greater than 160 Mev there is some indication of ‘ missing energy * and 
therefore, of charge exchange. 

The conservation of energy may also be expressed, for these disinte- 
grations, by the relation: 


140+ BE gar t1400,+7(140+H,),  . - » - (5) 


where H,, represents the kinetic energy of the primary meson in Mev; 
Egor, the total kinetic energy in the star including that of the neutrons 
and charged mesons; #,, the average kinetic energy of the ej ected card 
a-mesons, and also, it is reasonable to assume, of the ejected neutral 
mesons; and 7,=average number of ejected neutral mesons. 


Solving for 7, we obtain 


— Estar = Ea a 
i= {(Fae) +07} / (1+ Es): ie. eer GC) 
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We have calculated 7, for three ranges of primary energy, viz- 
0-160 Mev, where it .is already known that some absorption occurs, 
160-800 Mev and 800-1100 Mev. The results are given in table 8. 
Although subject to large statistical errors the figures give an indication of 
the relative probabilities for the different processes which can take place 
in meson—nucleon collisions. It can be seen that approximately 50%, of 
the incident mesons suffer charge exchange, irrespective of their energy. 

It is possible, in principle, to observe the electron pairs arising from the 
gamma rays which are the decay products of the neutral 7-mesons. How- 
ever, the frequency with which such conversion electrons appear near the 
parent star must be small, and a search for them would be very arduous in 
the present experimental conditions. 


Table 8. The Relative Probabilities for the Scattering Conversion or 
Absorption of Energetic ~—Mesons 


Energy range | Scattering | Conversion | — ,— | Absorption | Multiple 
3 = on NaN ms a= 
(Mev) Nx No ° | 1—(m,+7p») | processes 
0-160 0-26 0-38 0-64 O37, — 
160-800 0-40 0-57 0-97 — — 


800-1 100 0-75 0-70 1-45 — ~0-45 2 


The statistical errors associated with the values of n, have been given in 
table 6 ; those associated with the values of n, are of the same order of magnitude, 
ie. ~30%. 


(iii) Absorption 

It is difficult to prove the absorption of mesons of great kinetic energy 
because the contribution from the rest-mass is relatively small. The 
results shown in table 8 indicate, however, that no absorption takes 
place at energies above 160 Mev. For mesons of kinetic energy less than 
160 Mev, the absorption of a meson can be definitely established in three 
disintegrations as a result of the ejection of a proton with a kinetic energy 
greater than that of the primary meson. Assuming that an equal number 
of fast neutrons are ejected in a similar manner, the fraction of mesons 
absorbed in this energy range is ~50°%. This is in agreement with the 
value derived in table 8, and with the data of Bernardini et al. (1951), 
who have shown that when a beam of 110 Mey 7 -mesons passes through 
matter about 65°, of the particles are completely absorbed by the nuclei 
which they encounter. 
(iv) The Multiplication of Mesons 


We have observed one event in which two mesons, of kinetic energy 
365 and 375 Mev respectively, were ejected from a nucleus disintegrated 
by aa-meson, of 1050 Mev kinetic energy. Energy considerations suggest 
that the production of the extra meson (one presumably being the primary 
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particle with reduced energy) cannot be explained in terms of meson 
production by a nucleon ‘knocked on’ by the primary meson*}. However, 
the event may be explained if it is assumed that 7-mesons can create 
further mesons directly in collisions with nucleons. The cross section for 
this process of ‘ meson splitting ’, at energies of the order of 1 Bev, has 
been estimated by Heitler and J anossy to be about 1/10 of that for the 
large angle scattering. 

Further evidence for the reproduction of mesons in meson—nucleon 
collisions has been put forward by Rosser and Swift (1951). They have 
exposed photographic emulsions above and below 30cm of lead, at an 
altitude of 2860 m, and find that for showers of n,=2 the ratio of charged 
to neutral primary particles is much greater for the showers observed under 
the lead, than for those above. This result is confirmed by similar 
observations made from plates exposed at an altitude of 3540 m 
(Camerini, Lock and Yekutieli 1951). Rosser and Swift suggest that the 
increase in the number of showers with charged primaries is caused by 
7-mesons created in the lead assembly, which produce showers of 7-mesons. 


CONCLUSIONS 

If the disintegrations produced by fast 7-mesons are to be explained in 
terms of elastic meson—nucleon collisions, the average number of such 
collisions in the nucleus must be between two and three. This result is 
obtained from, (i) a comparison of the size distribution of meson and 
proton stars; (ii) an analysis of the recoil-protons in meson and proton 
stars; and, (ili) a study of the characteristics of the secondary mesons 
which emerge from the collision. These observations refer mainly to 
stars with V,>3. However, the geometrical cross-section for the meson- 
nucleus interaction, derived from an analysis of the double stars, is 
consistent with this average number of collisions. On the other hand, the 
meson—nucleon cross section of 1-3 x 10-26 em? obtained by Chedester et al. 
(1951) for 85 Mev 7~-mesons corresponds to only one collision per nucleus. 

In the low energy range (£,<160 Mev) there is direct evidence for the 
absorption of about half the mesons which collide with nuclei. Absorption 
does not appear to occur at higher energies. At all values of the primary 
energy about 50% of the incident mesons suffer charge exchange. A 
definite case of ‘ meson splitting ’ has been observed at an energy of the 
order of 1 Bev, but the available data does not allow an estimate to be made 
of the frequency of occurrence of this type of multiple process. 
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SUMMARY 


This paper describes a scintillation spectrometer, which measures the 
energy of the pairs produced by gamma-rays. Each gamma- energy gives 
a symmetric peak in the energy distribution. The width at half-maximum 
is 7% at 4mev. The greatest advantage of this apparatus is the high 


sensitivity. It has been possible to measure sources with a strength of 
about 0-1 pe. 


§1. INTRODUCTION 


In measuring high-energy gamma-rays with a one-crystal scintillation 
spectrometer a complicated pulse height distribution is obtained when 
the crystal is a substance containing elements of high atomic number, as, 
for instance, thallium-activated sodium iodide. For every gamma-energy 
a photoelectron peak and a continuous Compton distribution are 
obtained. The electron—positron pairs of each gamma-energy cause 
several peaks in the distribution. When the annihilation radiation from 
the positron is captured by the crystal the entire energy of the gamma 
' quantum is given off to the crystal. If one or both annihilation quanta 
escape the crystal, the energy is reduced by 0-51 and 1-02 Mev, respectively. 
Pair production thus gives rise to three peaks for each gamma-energy in 
the pulse height distribution. One of thése peaks coincides with the 
photoelectron peak. Gamma-rays with several energies therefore produce 
a very complicated pulse height distribution and in such cases accurate 
measurement of the energy can be very difficult. Naturally a substance 
with elements of only low atomic number may be used as the crystal, for 
example anthracene. The photoelectric effect and pair production have 
then a very low cross section and the pulse height distribution consists 
solely of Compton distributions. The difficulty then arises of determining 
the end point of a continuous distribution, which is especially hard in 
this case, on account of the comparatively low resolution of the scintillation 
spectrometer. 

It is possible, however, to employ the two annihilation quanta emitted 
from the positron to separate the pair production from the other processes 
in a sodium iodide crystal. The counter which is used for measuring the 
energy is surrounded in two opposite directions by two other counters 
eas erasure Fete es ee 
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designed to capture the annihilation radiation. Coincidences between 
these two counters are used to gate the pulses from the central counter. 
The pulse height analyser thus gives the distribution of those pulses in 
the central counter which are simultaneous with pulses in the two other 
counters, so that the distribution of the pulses due to the electron—positron 
pairs is obtained. The pairs from a certain gamma-ray all give off the 
same energy to the crystal, and each gamma-energy thus corresponds to 
a symmetric peak in the pulse height distribution. This facilitates a 
careful measurement of energy and also a study of gamma-rays with 
several energies. 

A few measurements undertaken with such an apparatus were described 
in a previous report (Johansson 1950b). The present paper describes 
more fully an improved apparatus constructed according to the same 
principle. The idea of using annihilation quanta to select the pair 
production pulses has also been advanced by Hofstadter and McIntyre 
(1950). Bair and Maienschein (1951) have quite recently published 
measurements with a similar apparatus. 


§2. APPARATUS 


Fig. | shows a schematic diagram of the apparatus. A central crystal 
(C,) with the dimensions 131320 mm is surrounded by two side 
crystals (C, and C;) with the dimensions 35x 35x 20mm. Each of the 
three sodium iodide crystals is enclosed in a container. The containers - 
are nickel-plated on the inside and closed at one end by a window of perspex, 
which lies close to the top of the photo-multipliers. The multipliers are 
of type E.M.I. 5311. The pulses from the multipliers are coupled to 
three cathode followers, which are fastened to the box and transmitted 
further by means of cables to other parts of the apparatus. The pulses 
from the central counter are amplified in a linear amplifier with negative 
feedback and are coupled to a gate. The side counters are connected to 
a coincidence circuit and the coincidence pulses trigger a univibrator 
which opens the gate. The gated pulses, finally, are analysed in a pulse 
height analyser. 

- The size of the crystals depends on several factors. The central crystal 
should be comparatively large in order to obtain a high degree of sensitivity. 
Another reason which speaks in favour of a large crystal is that the range 
of the electrons and positrons must be only a small part of the linear 
dimensions of the crystal. Otherwise it may occur that a fairly great 
number of pairs give off only part of their energy to the crystal, thus 
causing a deformation of the peak in the pulse height distribution. The 
central crystal must not, on the other hand, be made go large that a 
greater part of the annihilation radiation is captured in the crystal, as 
this would produce reduced sensitivity. The two side crystals should 
be made as large as possible so as to capture the largest possible part of 


the annihilation radiation. This is essential for the sensitivity of the 
apparatus, 
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Only those pulses from the central counter which are due to pair 
production are accompanied by simultaneous pulses in the two side 
counters. The pulses which reach the pulse height analyser, therefore, 
are a measure of the energy of an electron—positron pair. They give a 
peak in the pulse height distribution. The two side counters naturally do 
not capture only the annihilation radiation but also scattered rays. hy 
the incident rays are not collimated these counters also count the 
radiation direct from the source. They may therefore attain a rather 


Pig. 1 


photomultiplier 
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Schematic diagram of the scintillation spectrometer. 


photomultiplier 


high counting rate. Accidental coincidences may then open the gate 
for all kind of pulses from the central counter. In this way one may 
obtain a continuous distribution in the pulse height distribution as a 
background to the pair peak. In order to avoid this, weak sources should 
be used or the rays be collimated. The coincidence circuit should also 
have a short resolving time. The circuit used in the present work is of 
a type described by Garwin (1950). In order to make the apparatus as 


$2 
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stable as possible the coincidence circuit was not set at the highest possible 
resolution, 10-7 sec proving quite sufficient. The pulses from the cathode 
followers are shaped with delay-lines, short-circuited at one end. The 
coincidence circuit is driven directly from the cathode followers without 
amplification. This is possible since the voltage of the two multipliers 
of the side counters is comparatively high, about 2000 v. The voltage 
of the multiplicator of the central counter, on the other hand, is low, 
only about 1000 Vv. The relation between energy and pulse height has 
been shown not to be linear at higher voltages (Johansson 1951). 

The gate and the pulse height analyser may be constructed in three 
different ways. The simplest one is, perhaps, to pass the pulses from the 
central counter to the vertical deflection plates of a pulse oscilloscope. 
The pulses from the coincidence circuit are then made to start the sweep. 
The oscilloscope then acts as a gate. On the screen only the pair pro- 
duction pulses are seen and they start from one and the same point. 
They have a very short rise time, and what is seen is the exponentially 
decreasing part. If the counting rate is sufficiently high the pulse height 
_ distribution can be studied directly on the screen. It is best, however, to 
have a photographic record of the screen. The time of exposure is chosen 
so that a large number of pulses is photographed. Since all the pulses 
start from the same point of the screen the images on the plate are super- 
imposed together and on a suitable adjustment of the brightness a picture is 
obtained where the density is a measure of the number of pulses. The 
oscilloscope is of the type Tektronix 511A. The sweep speed was 
0-1 cm/usec. As collimator lens in the photographing a Zeiss Sonnar 
f: 1-5 was used. The counting rate was usually about 1 pulse/second 
and the time of exposure 10 minutes. Each gamma-ray energy gives a 
diffuse band on the plate, corresponding to a peak in the pulse. height 
distribution. In addition, the zero line is reproduced, due to the fact 
that the accidental coincidences trigger the sweep without the simultaneous 
occurrence of a pulse on the vertical deflection plates. If an unknown 
gamma-spectrum and a line of reference are photographed on the same 
plate the pulse heights can be compared directly. Since all the pulses 
have the same time constant, this comparison can be made along an 
arbitrary line at right angles to the zero line. It is better, however, to 
determine the density along a line at right angles to the zero line by means 
of a microphotometer. Fig. 2 shows such a densitometer trace for gamma- 
rays from a RdTh-source with a strength of 0-1 me and placed at a 
distance of 20 cm from the crystals. The rays were collimated through 
16 cm of lead. In the figure the zero line (A) together with a well-defined 
peak (B) are seen, This peak is due to the gamma-ray of 2-62 Mey from 
ThC”. From a densitometer trace of this kind the pulse height: can be 
determined with good accuracy. The pulse height is proportional to the 
energy and therefore the energies of gamma-rays can be measured after 
calibration with some known gamma-energy. Measurements of a few 
well-known gamma-energies agreed well with the accepted values. 
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; On many occasions it may be necessary to determine the height of 
single pulses from the spectrometer, for instance, when extremely weak 
sources are used and when a quantitative comparison of the intensity 
of various energies is desired. The spectrometer is therefore equipped 
with an electronic pulse height analyser. The pulses from the coincidence 
circuit trigger the univibrator, which gives rectangular pulses with a 
length of lysec. These pulses open an electronic gate for the pair 
production pulses. The rise time of the latter is very short and the time 
constant comparatively long (100 usec). When the gate is closed after 
I ysec they are cut off, becoming rectangular in shape with an almost 
flat peak. This form enables the use of a fairly simple single-channel 
differential discriminator as the pulse height analyser. The latter consists 


Fig. 2 
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Densitometer trace showing the pulse height distribution of a RdTh-source. 


of two Schmitt triggers together with an anti-coincidence circuit. The 
bias voltages are set by means of a potentiometer. A constant difference 
between the triggering levels is obtained by means of batteries. This 
device thus counts the number of pulses in a certain pulse height interval 
of variable size and position. The pulses from this analyser, finally, are 
counted with a conventional scaler and a register. 

The oscilloscope may be used as a second means of measuring the 
height of single pulses. The pulses from the central counter, which in 
this method have a time constant of 1 usec, are streched by means of 
a pulse stretcher to 5 psec. The stretcher consists of a delay-line, short- 
circuited at one end and charged via a number of crystal rectifiers (Craib 
1951). The pulses obtained are almost rectangular in form with a flat peak. 
They are coupled to the vertical deflection plates of the oscilloscope and the | 
latter is set for recurrent sweep. The sweep speed is kept fairly low 
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(0-1 mmjusec). The peaks of the pulses are then seen as bright spots 
on the oscilloscope screen. The brightness is then lowered by means of 
the potentiometer regulating the grid voltage of the cathode-ray tube so 
that pulses are no longer seen on the screen. The rectangular pulses 
from the univibrator are set at a length of 5 usec and are coupled to the 
cathode of the cathoderray tube. They have a negative sign and the 
brightness of the screen has therefore an increased value for a time of 
5 sec. Thus on the screen only those pulses from the central counter 
are seen which are accompanied by a simultaneous pulse from the coinci- 
dence circuit. The brightness of the cathode-ray tube acts as a gate and 
only the pair production pulses are visible on the screen. As in the 
preceding case, a number of accidental coincidences are also obtained 
with the result that a number of bright spots appear along the zero line 
when the intensity is increased without the simultaneous appearance of 
a pulse on the vertical deflection plates. If the screen is photographed, 
the peak of each pulse appears as a spot on the plate. The height of the 
pulses is obtained by measuring the distance between these spots 
and the zero line. The time of exposure is made so long that the greatest 
possible number of spots are obtained, not so many, however, that they 
coincide. In order to get a sufficient number of pulses for the measure- 
ments it may be necessary to make several recordings. As a rule a 
satisfactory picture of the pulse height distribution can be obtained by 
measuring 3 to 4 plates. 

The last two methods give results that agree well. Fig. 3 shows a 
pulse height distribution that is representative for both methods. It 
shows the number of pulses as a function of the energy. The central 
counter was calibrated prior to the measurements with some well-known 
gamma-rays as previously described (Johansson 1950 a, 1951). The 
energy seen in the figure is the energy of the pairs (the energy of the 
gamma-rays minus 1:02 Mev). The peak A is due to a RdTh-source. 
The strength of the source was 0-1 me, and its was placed at a distance 
of 20 cm from the crystals. The radiation was collimated through 16 em 
of lead so as to avoid too high a counting rate in the side counters. When 
weaker sources are used measurements can be made without collimation 
although they require more time. The main part of the gamma-rays 
from the source has an energy of 2-62 Mev. Other lower energies do not 
cause any peaks in the dicts vHeE partly on account of lower intensity 
and partly on account of the great decrease in the cross section of the pair 
production at lower energies. The peak B comes from a Po—Be-source. 
In the nuclear reaction of *Be(«,)C high energy gamma-rays are 
emitted. According to the latest measurements their energy is 4-45 Mev 
(Pringle e¢ al. 1950, Bradford and Bennet 1950, Bell and Jordan 1950). 
The source emitted a total of approximately 104 quanta per second. It 
was mounted close to the central crystal without collimation. The energy 


was determined to be 4:45-+0-05 Mev, which agrees well with previous 
measurements. 
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The spread of the pulse height around a mean value is due largely to 
statistical fluctuations in the number of electrons emitted from the 
photo-cathode of the multiplier. The statistical fluctuations in the gain 
of the multiplier also contribute to the spread. Theoretically, the peaks 
should have the form of Gaussian distributions with the widths at 
half-maximum proportional to the root of the energy. The peaks in 
fig. 3 are fairly easily approximated as Gaussian distributions. The 
widths are 10 and 7%, respectively. The relation between these figures 
~ Is on the whole what might be expected theoretically. 


Fig. 3 
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Pulse distribution of a RdTh-source (A) and a Po—Be-source (B). 


§ 3. DISCUSSION 


The apparatus described in the foregoing measures the energy of ae 
electron—positron pairs produced by gamma-rays and is thus a type 4 
pair spectrometer. It may be of interest to compare it with i ees 
pair spectrometer, for example the one described by Walker and McDanie 
eae to what has been said above the resolution Se CEC ie the 
width at half-maximum of the apparatus described here is 7% at 4: wee 
improving with increased energy. At 10 Mev it should pig ly 
4%, The value for the magnetic pair spectrometer is given as ' a e 
Des types of spectrometer have thus approximately the same resolution. 
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The lower limit of the scintillation spectrometer is determined by the 
cross section of the pair production. At too low energies the measurements 
take a long time. The decrease in the resolution, on the other hand, is 
less significant. The peaks in the pulse height distribution retain their 
symmetrical shape, allowing a fairly accurate energy determination also 
at low energies. The lower limit for the spectrometer probably lies 
around 1-5 Mey. | 

The greatest advantage of this type of scintillation spectrometer is its 
great sensitivity. No collimation is necessary and the crystal captures a. 
large part of the incident rays. Even very weak sources can therefore 
be measured. The Po—Be-source used in the present measurements 
emitted a total of approximately 104 high-energy gamma-rays per second. 
The counting rate of the spectrometer was then 1 pulse per second, 
i.e. 10-4 of the gamma-quanta emitted were counted. At higher energies 
the sensitivity is increased, due to the increasing cross section of the pair 
production. For the magnetic spectrometer the peak counting rate 
per single counter pair is reported to be 10~* of the total emitted quanta 
at 17 Mev. The scintillation spectrometer is thus considerably superior 
with regard to sensitivity. For practical reasons it may sometimes be 
impossible to place the source as close to the crystal as in this measurement 
with a Po—Be-source. If the source is moved 5 cm from the crystal the 
sensitivity decreases with a factor of 10. 

The gamma-spectrometer described in this paper is fairly simple to 
construct. Its great sensitivity makes it useful in studies of high-energy 
gamma-rays of low intensity, for example from nuclear reactions. 
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XX. CORRESPONDENCE 


Some Observations on Slip Lines in Iron 


By H. W. Paxton, M. A. Apams and T. B. MassaLsKI 
University of Birmingham* 


[Received November 26, 1951] 


GREENLAND (1937) was the first to observe that the appearance of slip 
lines on mercury single crystals depended on previous mechanical 
straining. Brown and Honeycombe (1951) have recently studied the slip 
lines in single crystals of the face centred cubic materials aluminium, 
copper, lead and gold, and shown that the form of these lines is sensitive 
to the state of preparation of the surface. 

The purpose of this note is to point out that slip lines on the polished 
surface of body centred cubic iron single crystals behave similarly. The 
crystals were grown from decarburized mild steel (Si 0:1%, S 0-05%, 
P 0-05%, Mn 0-3%) in the form of wire 2 mm in diameter using the 
strain anneal method. About 0-003° C was introduced afterwards, and 
so the crystals showed the yield phenomenon (Cottrell and Churchman 
1949, Schwartzbart and Low 1949). The gauge lengths were approximately 
5 cm and the crystals were extended with a strain rate of 2x 10-4 per 
minute in a rigid beam loading: apparatus. For testing, each crystal 
was soft soldered axially into mild steel endpieces attached to the loading 
shackles through flexible Bowden cable and ball seatings. The upper and 
lower lengths of Bowden cable were of opposite twist to eliminate 
rotational torque. 

Fig. 1 (Plate VII +) is atypical illustration of the type of slip lines produced 
by an extension of 0-5° ona mechanically polished surface or on an electro- 
polished surface which had been lightly rubbed with fine polishing alumina. 
This extension was sufficient to take the specimen through the yield point. 
The separation of lines is of the order 5-10. If, however, the crystal 
face was given a preliminary electropolish (using a mixture of perchloric 
acid and acetic acid) after grinding the metallographic surface, the slip 
lines that can be produced by 0-5°% extension are extremely difficult to 
‘gee with an optical microscope and much closer together, as shown im 
fig. 2 (Plate VII). The separation is approximately 2 pu. : 

It is necessary to polish electrolytically to beyond a certain depth 
(about 10 in the present specimens), presumably to remove the cold 
worked layer, if the coarse slip lines are to be avoided. When the electro- 
polishing is not as deep as this the specimen appears as in fig. 3 (Plate VII). 
De cen ee 

* Communicated by Professor A. H. Cottrell. 
+ For plates see end of issue. 
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Here both types of slip occur on the same surface, the coarse slip being 
limited to following a straight band across the crystal which is the trace 
of a polishing scratch. It should be noted that the two types of slip are 
continuous, the coarse lines branching at the boundary of the band to form 
the fine ones. 

Fig. 4 (Plate VIII) is a further photograph of this effect, taken slightly 
out of focus to show the traces of two scratches intersecting. The coarse 
slip lines appear on the grooves ; the fine slip lines are not in focus. 

In fig. 5 (Plate VIII), the approach of coarse slip lines to a grain boundary 
isshown. Here the coarse slip branches about 20 from the boundary and 
no visible slip is seen at all within about 5 y from it. Slip took place in the 
other crystal, but was not in the right orientation to be shown in this 
photograph. 

The authors wish to express their appreciation to Professor A. H. Cottrell 
for helpful discussion. They are also grateful to Professor D. Hanson 
and the University authorities for provision of facilities and financial 
support. 
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An Attempt to Produce a Thermo-Nuclear Reaction in Deuterium by Means 
of a High Current Spark Discharge 


By P. Reynowps and J. D. Craces 
University of Liverpool* 
[Received November 21, 1951] 
THE two reactions in deuterium given below have a high probability 


of occurrence, and are therefore suitable for an attempt to produce a 
thermo-nuclear reaction. 


2+ H?-> He?+n! +3-10 Mev 

H?+ H?-—> H® + H!+3-75 mev. 
According to Gamow and Critchfield (1949), in order that these reactions 
should proceed at an observable rate, a temperature in the region of 
106 °c, and a degree of ionization of the gas of nearly 100°, are required. 
A generator had been built in this laboratory for a general investigation 
of the physical conditions in high current transient discharges, and it 
was thought to be worthwhile to investigate the effect of passing sparks 


* Communicated by Professor H. W. B. Skinner, F.R.S. 
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carrying currents in the region of 300 ka through deuterium. It was 
hoped that the self magnetic field of the high current would produce 
a constriction of the spark channel, and that this would so increase the 
voltage across the spark, and reduce the heat losses from it, that much 
higher temperatures, and almost 100% ionization of the gas would result. 
It was realized that the conditions necessary for efficient thermo-nuclear 
reactions could hardly be reached in a spark. In view, however, of the 
uncertainty of the shape of the positive ion energy distribution, especially 
for high energies (although deviations from a Maxwellian distribution 
in the present discharge conditions are extremely unlikely due to inter- 
action effects considered by Cahn (1949)) and the unknown effect of high 
self fields on high pressure discharges, it was considered worthwhile to 
attempt the experiment. 

The high current generator, a description of which will be published in 
full, consisted of a bank of capacitors of a total capacitance of 58-5 pF, 
charged to a voltage of 25 kv, and then discharged in parallel through 
a circuit of low impedance. A damped oscillatory discharge was obtained 
in which the current had a maximum value of 285 ka, and a period of 
29 usec. The tungsten electrodes of the spark gap were 2 cm apart, and 
were enclosed in a spark chamber containing deuterium at atmospheric 
pressure. 

It would have been impracticable to have detected the emission of 
protons from the reaction, since the detector would then have been inside 
the chamber and so exposed to the heat and the shock wave from the spark. 
Attention was therefore concentrated on the detection of neutrons from 
the reaction. Since the yield was expected to be extremely small, a very 
sensitive detector was required. One using a Szilard-Chalmers reaction 
in potassium permanganate, and a technique developed by Rotbiat 
(1941) were chosen, since neutrons could be detected over a large solid 
angle and the effect of many sparks integrated. Also unlike counter 
methods, it was free from problems caused by the large electrical 
interference from the spark. 

The volume of 22 litres of an almost saturated solution of potassium 
permanganate was contained in glass vessels placed around the spark 
~ chamber, so that a solid angle of approximately 2 7 radians was subtended 
by the detector at the spark. A calibration with a Ra—Be neutron source 
showed that a total of 3 x 10° neutrons incident upon the detector produced 
a manganese activity which gave a counting rate equal to the background 
count. ee 

A blank experiment, performed with no neutron sources in the vicinity, 
showed that the activity of the manganese increased the background by 
25%. This activity was attributed to cosmic ray neutrons. 

A total of 100 sparks, each carrying a peak current of 285 ka were 
passed through the deuterium in the space of 1 hour. No activity in 
excess of that in the blank experiment was detectable. The total number 
of neutrons emitted by the sparks must therefore have béen less than 
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6105, Assuming a Maxwellian distribution of velocities amongst the 
ions, the ion temperature in the spark channel must have been less than 
LOS <G. 

Measurements are now being made by Allen (1950) on the current 
densities and energy dissipation in high current sparks in various gases and 
experimental conditions. Recent relevant work, theoretical and 
experimental, on the effect of self magnetic fields of discharges has been 
carried out by Thonemann and Cowhig (1950). 
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Film Transfer in Helium IT below 1° k 


By E. AMBLER and N. Kurri 
Clarendon Laboratory, Oxford* 


| Received December 20, 1951] 


We have recently carried out experiments on the rate of transfer in the 
helium {I film at temperatures down to about 0-15° K and although the 
work is still in its preliminary stages it seemed worth while to report 
briefly not only the results obtained but also the experimental technique 
used which has enabled visual observations to be made on liquid helium 
in the range of temperature attained by the magnetic cooling method. 
From our experience we believe that measurements on other characteristic 
properties of helium II such as film thickness and viscosity can be extended 
into this range of temperature without too much difficulty. 

The transfer rate has hitherto only been measured down to about 
I’ k, that is a reduced temperature of 7/7',=0-4. While most deter- 
minations have indicated a temperature variation roughly proportional 
to [1—(T/T,)"| where o=6 (Mendelssohn and White 1950), and a 
corresponding flattening of the transfer-vs-temperature curve below 
about 1-5° kK we thought it worth while to extend these determinations 
to cover lower temperatures. Quite apart from the desirability, on 
general ground, of testing this law to smaller values of the reduced 
temperature, the range below about 0-7° kK is physically, in one respect, 
different from the higher temperature range inasmuch as the vapour 
pressure of helium becomes so small (e.g. 1-6 10-5 mm Hg at. 0-5° K) 


a Ee ee ee) ee ee 
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that-any possible influence of the vapour phase on the transfer rate by 
evaporation and recondensation or by heat conduction must become 
negligible (Kurti and Simon 1938). It should also be mentioned in this 
connection that the velocity of second sound increases sharply at about 
0-5°k (Atkins and Osborne 1950) and the temperature variation of 
specific heat (Hull, Wilkinson and Wilks 1951, Kramers 1951) and of 


the heat conductivity (de Klerk 1946) changes fairly abruptly at about 
0-6° K. 


Fig. | 
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The method of cooling liquid helium in an ‘ open ’ container by adiabatic 
demagnetization given by Hudson, Hunt and Kurti (1949) was het 
these experiments. The container used for measuring the Cae er ie 
is shown in fig. 1. It is suspended by means oF nylon threads * the 
vacuum jacket of a glass demagnetization cryostat (not aro senting 
liquid helium at about 1° kK. The right-hand side sae ) is a 
longitudinal section and the two left-hand side digs 4 oe @) be 
cross sections. Film transfer takes place from the nese se are 
to the outer glass vessel O and is limited by the constriction ©, 1-2 mm 
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diameter and 5 mm long in the upper part of I. The inner beaker contains 
a compressed pill P, of manganous ammonium sulphate in the form of a 
cylinder split in half longitudinally by a 1 mm wide slit §,. The space 
between I and O is filled with a closely fitting annulus made from a pill 
of the same salt having a 1-5 mm wide slit S, in the lower half in line with 
the slit S,. In order to be sure that the constriction C is at the same 
temperature as the salt, P, completely surrounds I over this portion 
of tube and is stuck to it by means of a cold setting plastic. Helium 
may be condensed into I through the fine capillary L. A further pill 
attached to L above the transfer vessel acts as a thermal shield (Hudson, 
Hunt and Kurti 1949). 


Table 1 
Run | Temperature | Transfer rate Run | Temperature | Transfer rate 
No. ok (cm*/cm/sec) No. ae (cm3/cm/sec) 
] 0-15 15-954 10s" 4 1-46 13-4 10> 
0-17 15-9 1-39 13-0 
0-17 15-9 vA: 12-6 
1-03 12-6 : 12-6 
1-23 11-0 
1-66 10-2 f 16-4 
2-02 4:8 - 16-9 
‘ “TE 12-8 
2 0-20 15:8 | )- 12-1 
0-99 10-9 : 16-7 
1-20 10-8 : 15-7 
2 : 14-4 
3 0-16 16-4 “9% 12-2 
0-5 14-4 “! 13-0 
1-42 12-2 
1-03 12:3 


The demagnetization cryostat was enclosed in a vacuum jacket 
surrounded by a Dewar vessel containing liquid hydrogen. Apart from 
vertical slits through which the illumination and viewing of the helium 
took place, the glass walls of the cryostat, vacuum jacket and Dewar 
vessel were silvered to minimize the radiation falling on the pills. The 
slit 8; was illuminated by means of a narrow parallel beam of light 
from a small mercury lamp, suitable filters ensuring that only the green 
light was used ; stops limited the cross section of the beam to about the 
size of the slit 8,._ In this way we minimized the amount of light falling 
actually on the salt, a precaution made necessary by the apparently 
strong absorption by the salt of visible light, which, as was shown in 
preliminary tests would have given a prohibitively large heating up rate. 
For the same reason the slit was only illuminated during the actual time 
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of observation which, by practice, was reduced to a few seconds. In these 
conditions the transfer vessel remained below 1° K for about half an hour, 
corresponding to an average heat influx of 2 000 ergs/min. 
__ The experimental procedure was as follows. A calculated quantity of 
helium was condensed into the inner beaker so that after demagnetization, 
which, from an initial H/T of 10 kilo-oersted/degree led to about 0-15° x, 
the meniscus M, would be near the top of the slit S,. After demagnetization 
the position of M, was determined every two or three minutes and, during 
the intervening periods, the susceptibility of the salt from which the 
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absolute temperature could be estimated (Cooke 1949) was measured by 
the ballistic method. From the rate of fall of M, and the free cross-section 
of I the creep-rate through C was calculated.* ~ 

The results are given in table 1 ; it will be seen that the rates above 1” k 
are somewhat higher than those usually observed for glass eee and 
Mendelssohn 1939), and there are differences of up to 15 /o between 
different runs. This may have been caused. by slight contamination of 
the constriction in the inner beaker (Bowers and Mendelssohn 1950). 
For the purpose of comparison each set of results from different runs 


* To a transfer-rate of 10 x 10-5 cm3/em/sec corresponds a meniscus movement 
of 0-30 mm/min. 
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‘was ‘ normalized ° (by multiplication by a constant factor) so as to give 
identical transfer rates at 1:-2°K. The results so obtained are plotted 
in fig. 2 where it will be seen that after the flat portion between 1-5° k 
and 0-8° K the transfer-rate shows a further increase which amounts to 
about 30% at 0-2° K., 

We have considered whether this increase may be due to some spurious 
effect such as a preferential heating of the inner beaker caused by the 
illumination and consequent distillation into the outer vessel, but the 
effect seems to be too large to be explained in this way. Before drawing 
any conclusions as to the fundamental implications of this result we 
should prefer to have more measurements made on the film both on 
other surfaces and under different experimental conditions at these very 
low temperatures. It is worth noting, however, that the results of. 
Burge and Jackson (1951) do show an increasing film thickness as the 
temperature falls below 1-4° kK. 


We wish to express our grateful thanks to Professor F. E. Simon for 
his stimulating interest, to Mr. W. E: Gardner for his untiring help with 
the experiments and to the Department of Scientific and Industrial 
Research for E. A.’s maintenance grant. 
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Fine Structure in the U238 Decay 


By Barpara Zasac 
with Appendix by N. MILLER . 
Department of Natural Philosophy, University of Edinburgh* 


| Received December 27, 1951] 


Inrorp G5 nuclear emulsions have been impregnated with natural 
uranium in solution as the complex ammonium citrate, adjusted to pH 8. 
Prior to conversion into citrate the uranium was freed from UX, and other 
decay products and impurities by extraction of the crystalline nitrate into 
ether, followed by washing of the ethereal solution with a little water. 
The plates were exposed for 72-5 hours. 


* Communicated by Professor N. Feather. 
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The developed plates were examined for the occurrence of electron 
tracks in association with the tracks of the a-particles of the three isotopes 
of natural uranium. <A statistical survey showed that such electron tracks 
were associated with 29-1+1:3% of all uranium a-particle tracks. 
Measurements were made on 776 «-particle tracks to give distribution of 
ranges for ‘all’ «-particles, and on the tracks of 635 a-particles having — 
associated electrons. The two range distributions were very similar. 


Longer tracks 
(not included in the figure) 


SIS 
52-54 
wn 55:5 ua 
oe 58 5ue 
: #2 
4a 5 
5 ey ye 
ro) 
i 
o 
a 
= 
= 
Pee 
10 


10 20 30 


Range (microns). 
Distribution of ranges of electrons associated with a-particles of ranges 
smaller than 17-55 wp. 


In both two main peaks occurred at about 16-5.u and 19-5, emulsion 
range and there were signs of a small (incompletely resolved) peak at 
about 17-5 p range due to the 4-4 Mev «-particles of U235 (Baldinger and 
Huber 1949, Ghiorso 1951). The similarity of the two «-particle range 
distributions showed conclusively that secondary electrons follow the 
x-disintegration of U238 as they do that of U234. 

The results of measurement of electron tracks associated with U234 
«-disintegrations were in general agreement with those reported previously 


by Teillac (1950). 
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The diagram represents the histogram (and the corresponding smoothed 
curve) of ranges of electrons associated with «-particles having ranges 
smaller than 17-55, i.e. with the «-particles of U238 and U235. The 
total number of events represented in the figure is 232. A careful esti- 
mation showed that less than a quarter of these could involve UX, and 
UY disintegrations, following the U238 and U235 disintegrations respec- 
tively, or internal conversion of the known y-rays of U235. The high 
energy tail of the distribution might well belong to these events, but the 
main part of the low-energy component, that is fully three-quarters of the 
electrons, must be a feature of the U238 decay. 

Sixteen events are included in which more than one (two, three or four) 
electrons were associated with a single «-particle. All of these could be 
interpreted in terms of subsequent UX, and UX, (Heershap, Ong Ping Hok 
and Sizoo 1950), or possibly UY, disintegrations. 

The main component of the electron range distribution shows two 
maxima at 5-24 and 12, range, corresponding to electron energies of 
27-5 and 45 kev (Zajac and Ross 1949) respectively. Electrons having 
_ these two energies may be interpreted as L and (M- NV) internal conversion 
electrons of a y-ray of about 48 kev energy. 

A calculation taking into account all the facts leads to an estimate that 
23+3%, of U238 «-disintegrations are associated with an L or (M ae 
internal conversion transition of 48 kev total energy. 


APPENDIX 


For emulsions of this type to remain fully electron sensitive it is essential 
that any impregnations be carried out in neutral solutions (Picciotto 1949). 
If impregnations with heavy elements are desired, this usually means in 
practice that a complexing agent has to be added. Citrate ions are found 
to be satisfactory for a wide range of heavy elements, and also prevent 
the precipitation of carbonate within the emulsion during development 
when barium carriers are added. 
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A disintegration initiated by a 7-meson of kinetic energy 640 Mev from which 
a a-meson of energy 40 Mev is ejected, together with a proton of 380 Mey. 
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OBITUARY 


ALLAN FERGUSON, 1880-1951 : 


AN APPRECIATION 
BY H. R. Roprnson 


ALLAN HrrcHEN FERGUSON was born at Entwistle in May 1880, and 
his boyhood was spent in and near mid-Lancashire. He got his first 
systematic instruction in science at the Harris Institute, Preston, where 
he began his studies for a London External degree. In 1902 he entered 
the University College of North Wales, Bangor, as an Exhibitioner, and 
in due course took the M.A. degree of the University of Wales in 
Mathematics, and the London B.Sc. degree in Physics. He obtained the 
London D.Sc. degree in 1915. 
_ From 1905 to 1919 he served on the physics staff at Bangor; during 
the first world war he was also for some time in charge of the applied 
mathematics courses, and acted as. lecturer in electrical engineering. 
During his later years in Bangor, he devoted much attention to hydro- 
electric problems, with particular reference to the possibility of developing 
water-power resources in North Wales. He left Bangor in 1919, and was 
for two years lecturer in physics at the Manchester Municipal College of 
Technology (in the Faculty of Technology of Manchester University). 
_ From Manchester he came in 1921 to join the staff of East London 
College (now Queen Mary College) in the University of London, as second- 
in-command to Professor C. H. Lees, F.R.S. He was made Assistant 
Professor of Physics in 1934, and he remained in the service of Queen 
Mary College until his retirement in 1945. While the College was in 
Cambridge in 1939-45, he assisted the depleted staff of the Cavendish 
Laboratory by lecturing to Cambridge Tripos candidates, in addition to 
his work for London students. 

This bare list of his teaching appointments, possibly of no great interest 
in itself, does lead to one point which needs more than a passing mention 
This is that for the most part his teaching life was. spent in departments 
that were far from lavishly staffed, and over the years he must have 
lectured at one time or another in all branches of classical physics, and 
at all stages of university degree courses. His special interests lay mainly 
in general physics and thermodynamics, but he was always ready to take 
over at the shortest notice an advanced honours course in optics or 
electricity—and he did it all supremely well. He was one of the best 
university teachers in the country, and there are today scores of practising 
scientists who hold him and his teaching in grateful memory. 

At Queen Mary College he took a large share in organizing the teaching 
of physics and in supervising postgraduate research. He was moreover 
an outstanding figure in the general life of the College ; he enjoyed to a 
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high degree the confidence of his colleagues on the teaching staff, and 
from 1934 to 1944 he represented the Academic Board on the Governing 
Body. In 1946 the College fittingly recognized his many services by the 
award of the highest honour in its gift, a Fellowship. : 

He had an unusually sound and comprehensive knowledge of physics 
in general, and of thermodynamics and the properties of matter in 
particular. Indeed, as early as 1906, G. H. Bryan, in the preface of his 
treatise on Thermodynamics (published Teubner, 1907) acknowledged 
the help given to him in its preparation by “Mr. Ferguson, BSe.”,, 
Ferguson’s own research work lay mainly in the fields of the physics of 
liquids and the difficult borderland between pure physics and physical 
chemistry. He was a recognized authority on surface tension, and he 
contributed extensively to the study of electrolytic conduction. He also, 
during the five or six years which preceded the outbreak of the second 


world war, gave much time to perfecting and applying an ingenious and: 


powerful electrical method which he had devised for measuring the specific 
heats of liquids and their temperature variations. All his research 
work was characterized by great elegance and thoroughness; he excelled 
in devising methods which were at once simple and accurate, and his 
published descriptions of them were models of clear and graceful exposi- 
tion. Simplified forms of some of his original experiments—notably one 
for the determination of the surface tensions of liquids available only in 
very small quantities—have been incorporated into the routine of under- 
graduate laboratory courses, and they never fail to impress and interest. 
the better students. 

His contributions to the advancement of physics were, however, far 
from being confined to work done within the walls of the laboratory. 
He was an active member of a number of scientific societies, and he 
rendered valuable services to these societies, and through them to the 
general body of physicists. This was particularly evident in his work 
for the Physical Society and the British Association for the Advancement of 
Science. He joined the Physical Society when he came to London in 
1921, and became a member of its Council four years later. In 1928 he 
accepted the office of Honorary Papers Secretary, which he retained until 
he was elected to the Presidency of the Society in 1938. It was about 
half-way through his period as Papers Secretary that the publication 
was begun of the Physical Society’s annual Reports on Progress in Physics, 
and it is no secret that Ferguson, if not indeed the “ onlie begetter ” of 
the series, was mainly responsible for its inception. He personally edited 
the first five volumes (1934-1938), and raised the series at once to the 
rank of a scientific periodical of international repute. 

The original intention was ‘ to present articles which shall discuss the 
latest developments in physical science, and which shall at once make 
clear to the non-specialist the meaning and extent of these advances, 
and provide the researcher in the particular field under discussion with 


a technical résumé helpful to him in his own work’, The need for a 
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series of authoritative summaries of this kind was obvious in 1934, for 
many branches of physics were then proliferating with bewildering speed ; 
the main difficulties to be expected in fulfillmg the plan were equally 
obvious—that of making the best choice of fields for reporting, and that 
of persuading active workers in those fields to help with the reports. 
Ferguson succeeded in the first place through his wide outlook and clear 
vision as a physicist, and in the second place quite simply because he was 
so widely known and well liked in scientific circles that he found it easy 
to get the right kind of help. The Reports were originally planned as 
no more than a service to the Fellows of the Physical Society, and as 
Ferguson modestly wrote in 1938, “....it is by a happy chance that 
they interest a public outside the membership of the Society’. It was, 
of course, through no chance at all that this happened, and the success 
of the Reports, and the added reputation they brought to the Physical 
Society, may well be regarded as a memorial to Ferguson’s qualities as a 
physicist and his gift for getting the best out of others. 

He was President of the Physical Society for the three years 1938-1941. 
He filled the presidential chair with distinction, and the Society owes 
much to his wise guidance during the critical early years of the war. 
From September 1939 onwards he was working in Cambridge, but the 
records of the Society show that he was regularly in the chair at.its 
London meetings—not only during his Presidency, but also frequently 
as Acting President after his own term of office had expired. This was 
done in spite of increasing physical disabilities, and especially impaired 
vision, which made travelling under war conditions very burdensome, 
and many of his friends felt and said that he was undertaking too much. 
He had his reward in seeing the Society’s activities vigorously maintained, 
under conditions in. which more timid leadership would have cut them 
down. 

He gave equally distinguished services to the British Association for 
the Advancement of Science. After being first a Secretary and then 
Recorder of Section A, he became an ordinary member of Council at the 
end of 1932, and from 1935 to 1946 he served as one of the two General 
Secretaries of the Association. He was President of Section A at the 
Blackpool meeting in 1936, when he tdok as the subject of his Sectional 
Address: ‘‘ Trends in Modern Physics”. It was a wise and witty 
address, and it’ is full of interest as a frank statement of the reactions of 
a philosophic mind, steeped in the physics of Kelvin and Maxwell, to 
the impact of a. newer synthesis. 

Ferguson also took a part in the activities of the Institute of Physics, 
serving as a member of the Board from 1928 to 1932, and again from 
1934 to 1938, and as a Vice-President from 1939 to 1943. He was. a 
member of the editorial board of Endeavour from its beginning in 1942, 
and for a short period in 1939-1940 he was an Advisory Editor to Nature. 

He had frequently expressed to me his affection for the Philosophical 
Magazine, and I remember the pleasure with which he received, in 1937, 
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an invitation to join its Editorial Board. For some years he was 
Executive Editor, until increasing physical disabilities compelled him to 
retire. One of his last direct services to the Philosophical Magazine was 
the editing of the Special Commemoration Number of 1948 (Natural 
Philosophy through the 18th Century, and Allied Topics). 

Another important task which he carried out after retiring from 
university teaching was that of advisory editor of the physics and 
engineering sections of the new Chambers’s Encyclopaedia. 

He was an excellent popular lecturer, and a good broadcaster on_ 
scientific and general topics. Many members of the Royal Institution 
will remember with pleasure a Friday Evening Discourse he delivered 
in 1934; his subject was “ Surface Tension’, and he showed (I think, 
for the first time at a public meeting) a beautiful slow-motion film of the 
formation of drops and splashes, which had been made with a high-speed 
camera taking nearly 2000 exposures a second. This development of 
high-speed cinematography had a special interest for him, and he 
lectured frequently on it, as in the British Association’s Alexander Pedler 
Lecture (Cardiff, 1940). 

To come now to the more personal side of his life and character: he 
was physically a very big man, standing several inches over six feet in 
height, and rather more than proportionally broad. In fact he tended to 
corpulence, and he would at times express a mild disdain for men who 
fell far short of the 20 stone which he affected to regard as the canonical 
weight for a physicist. He was, however, before illness overtook him, 
active and muscular, and in no way handicapped by his bulk. 

He was one of the friendliest souls alive, and his habitual expression was 
one of lively good-humour—which largely accounts for the ease of his 
relations with colleagues and pupils. Even the most diffident students 
were quick to realise that this large, impressive figure was infinitely 
approachable and benign, and he was particularly effective as an oral 
examiner, in getting the best out of nervous candidates. 

He was a voracious and omnivorous reader, and an incorrigible collector 
of books—especially old books. He had a very pleasant and commodious 
house in Bishop’s Stortford, where he lived for the last 22 years. In 
this, as an energetic amateur carpenter of the rough-and-ready school, 
he fitted open shelves into every available alcove and corner, and he 
finally had a large book-room, shelf-lined, built as an annexe, but his 
shelving always lagged far behind his acquisitions. Books were stacked 
all over the house, in a way that would not have been allowed in any but 
a tolerant and thoroughly book-loving household. His tastes were 
catholic, but naturally with a bias towards old science. He made no 
claim to a systematic knowledge of the history of science, but he was 
extremely well-informed on many sections of it; an instance of this 
may be seen in his delightful ‘ Introductory Survey” to the British 
ae ne Mapeah tA a eae ee of Science 

Sant storia science, he had the advantage of a 
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solid background of literary and general history ; he saw scientific 
discoveries not as isolated events, but in their relation to the general 
conditions of the times in which they were made. 

Although he did not confine himself to it, the 18th was his favourite 
century in literature, and Samuel Johnson his chief literary hero. In 
fact I think he liked to dwell in imagination on resemblances, real and 
fancied, between himself and Touneate There was at least one real 
difference between them, which was a blessing to his friends—Allan 
Ferguson was a patient and sympathetic listener, as well as a good talker, 
and he loved the give-and-take of round-the-table conversation. Although 
he had, for instance, an intimate knowledge of church architecture, he 
would never thrust the subject into a general conversation ; and although, 
- like nearly all men of his period, he had a store of favourite anecdotes— 
including, of course, a number of “ tacklers’ tales’ from his beloved 
Lancashire—he had none of the tiresome habits of the persistent raconteur ; 
in brief, he talked for enjoyment, and not for effect. 

In politics he was a Gladstonian Liberal, and for many years he was 
President of the Hertford Divisional Liberal Association ; less publicly, 
he was in sympathy with many aspects of socialism, and in some purely 
personal matters he was healthily tinged with conservatism. A man of 
strong courage, he was a hater of war and an active supporter of the 
League of Nations Union; and the failures of the League appeared to 
weigh more heavily upon him than any of his personal tribulations. 

He was above all things fortunate in his family life. He married in 
1919 Nesta Thomas, daughter of the Vice-Principal of the Normal College, 
Bangor. She had already proved herself a botanist of distinction, and 
had successively held a research scholarship and a Fellowship of the 
University of Wales, and a Lectureship in Botany in the Royal Holloway 
College (University of London). She continued her botanical work after 
her marriage, undertaking part-time university and other teaching and 
examining at various periods, and she obtained her London D.Sc. in 
1927—twelve years after her husband. Their only child, John, had an 
outstandingly brilliant record of achievement at school and at St. J ohn’s 
College, Cambridge, and, is now a Lecturer in Classics in the University 
of Durham. (He is, incidentally, the “John Ferguson, M.A., B.D. ’ 
who was co-author with his father of the introductory article in the 
Commemoration Number of the Philosophical Magazine, 1948.) ¢ 

Ferguson took great pride and pleasure in the academic distinctions 
and activities of his wife and son; more than this, they were a powerful 
consolation to him at a time when he most needed comforting. 

The story of the bodily afflictions of his later years, painful though 
it is, ought not to be entirely omitted; for it is also a story of high 
courage and continued activity under conditions that would have reduced 
most men to a sterile despair. His first serious troubles came in 1933 
and 1934, when he was threatened with the possibility of total blindness. 
This was averted by skilled surgery, and although the resulting condition 
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was something of a handicap in outdoor mov ement, he did not let it— 
restrict his literary and scientific activity. Early in the second world — 
war, however, there came signs of graver organic disease ; he was affected 
by a serious and progressive form of paralysis, followed by other painful | 
complications. I saw much of him in the earlier years of. this long 
illness ; we shared a room in the Cavendish Laboratory during a period 
when he was unable to lecture without the aid of powerful drugs. Never- 


theless, he insisted on continuing his work, and in spite of all he remained 


the sweetest of companions. He was not, of course, infinitely tolerant” 


of all things, and it was always possible to draw a shower of sparks from 
him by such simple devices as claiming an Englishman’s right to close 
a sentence with a preposition ; but I only once heard him grumble about — 
his own tribulations, and that was to make a purely jocular complaint _ 
that he could no longer compete on equal terms in spotting the bargains — 
on David’s bookstall. His intellect remained completely unclouded, and — 
even those who met him for the first time in this darkening period were 
deeply impressed by his personality and his mental power. For all who 
enjoyed his friendship—and they are very many—his pleasant voices, 

his nightingales, will be awake while memory lasts. 


H.R. R. 


The Directors of Taylor & Francis Ltd., proprietors of the Philosophical 
Magazine, would like to add a note telling of their indebtedness to Professor 
Allan Ferguson. He took over the Acting Editorship of this Journal on 
the death of Professor A. W. Porter in 1937. 

The war period was one of great stress for the Magazine. The difficulty 
of obtaining papers of the high standard expected of the Philosophical 
Magazine and the almost overwhelming difficulties of publication made 
Professor Ferguson’s task a very arduous one. It will be admitted that 
he overcame the various troubles extremely well, especially when it is 
remembered how badly handicapped he was with ill health during the 
last few years of his editorship. 

In addition to his work on the Philosophical Magazine his advice in 
many directions on policy was invaluable. 
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